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Abstract 

New physics searches at the LHC or Tcvatron typically look for a certain number of hard jets, 
leptons and photons. The constraints on the hadronic final state lead to large logarithms, 
which need to be summed to obtain reliable theory predictions. These constraints are 
sensitive to the strong initial-state radiation and resolve the colliding partons inside initial- 
state jets. Wc find that the initial state is properly described by "beam functions". 

We introduce an observable called "beam thrust" tb, which measures the hadronic 
radiation relative to the beam axis. By requiring re <C 1, beam thrust can be used to impose 
a central jet veto, which is needed to reduce the large background in — >■ WW — )■ iiyii? 
from tt —7- WWbb. We prove a factorization theorem for "isolated" processes, pp — t- XL 
where the hadronic final state X is restricted by <C 1 and L is non-hadronic. This 
factorization theorem enables us to sum large logarithms o^Iu^tb in the cross section 
and involves beam functions. The beam thrust spectrum allows us to study initial-state 
radiation in perturbation theory, which can be compared with experiment and Monte Carlo. 
We present results for the beam thrust cross section for Drell-Yan and Higgs production 
through gluon fusion at next-to-next-to-leading-logarithmic order. 

The beam functions depend on the momentum fraction x and the (transverse) virtuality 
t of the colliding partons. The t dependence can be calculated in perturbation theory by 
matching beam functions onto the parton distribution functions (PDFs) at an intermediate 
scale /xb ~ Vi- We calculate all the one- loop matching coefficients. Below jiB have the 
usual DGLAP evolution for the PDFs. Above /i^j the evolution of the beam function is in 
t and does not change x or the parton type. 

We introduce the event shape "N-jettiness" tat, which generalizes beam thrust to events 
with N signal jets. Requiring tat ^ 1 restricts radiation between the signal jets and vetoes 
additional undesired jets. This yields a factorization formula with inclusive beam and 
jet functions and allows us to sum the large logarithms In*" tat from the phase space 
restriction. 

Thesis Supervisor: Iain W. Stewart 
Title: Associate Professor 
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Chapter 1 



Introduction 

1.1 Effective Physics 

Over the centuries we have been studying the world around us, discovering new physics 
at new length scales. Starting from our everyday experience which ranges from scales of 
about 10^ to 10~^ meters (m), we encounter phenomena like gravity, contact forces, fluid 
dynamics and light. When we move on to the smaller scales of cells at 10~^ m or molecules 
and atoms around 10^"'^'^ m, we pass entire fields of science such as biology and chemistry. 
It is surprising that the underlying physics that leads to this rich variety of phenomena is 
simply (quantum) electrodynamics. The nucleus of the atom at scales of 10~^^ m is made 
of protons and neutrons and is held together by the nuclear force. This is a residual force 
coming from the strong force, also known as Quantum Chromodynamics (QCD). QCD is 
responsible for confining the quarks to protons and neutrons. The remaining force that is 
part of the well established Standard Model of particle physics is the weak force, which for 
example plays a role in some hadron decays. 

With the Large Hadron Collider (LHC) at CERN starting to take data, we will be prob- 
ing nature at smaller length scales than before and hope to discover new physics. In order to 
reconstruct new short distance physics from the measured leptons and QCD radiation, one 
typically looks for a signal with a certain number of jets of energetic hadrons. This thesis 
explores the effect of such restrictions on the hadronic final state in theoretical calculations, 
with an emphasis on the 0-jet case which allows us to study initial-state radiation. An 
introduction to the work in this thesis can be found in the next section. 

Our current understanding of nature as the strong, weak and electromagnetic force, 
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together with gravity, is beautiful in its simphcity as a theory. Yet its computational 
complexity makes it ill-suited to describe most phenomena in nature. In practice one 
adopts an effective description appropriate for the length scale one is looking at. We are 
already used to doing this in much simpler theories such as Newtonian gravity, where for 
our everyday purposes a constant gravitational acceleration is an appropriate description. 
Effectively we are treating the earth as an infinite plane, which is a reasonable approximation 
given that we are small compared to the size of the earth. If we instead would try to describe 
the motion of the planets in our solar system, the size of the planets and the sun becomes 
negligibly small and we can treat them as points. 

We encounter a similar situation in electrostatics, when we want to determine the electric 
field due to an object at a distance r that is large compared to the size s of the object. We 
can start by approximating the object as a point with a total charge Q and systematically 
improve this result by using the so-called multipole expansion. This is an expansion in 
s/r <C 1, where each next order contains more details of the charge distribution. The first 
correction comes from the dipole moment poi the configuration and is suppressed by a factor 
of s/r. The next correction comes from the quadrupole moment Qij and is suppressed by 
s^/r^ etc. We have thus found an effective description of the charge configuration, where 
the details of the full charge configuration are captured by only a few constants: Q, p, etc. 

This simplification to an effective description, where the relevant details of the complete 
physical picture are condensed into just a few constants, will be a theme that continues as 
we study effective descriptions of quantum field theories. In this thesis we will use Soft- 
Collinear Effective Theory (SCET) , which describes initial and final state radiation and 
provides us with a powerful tool for studying processes at the LHC. In the next section 



we will present an introduction to the work in this thesis. In Sec. 1.3 we will look at an 
effective description for a decay involving the weak force, which will allow us to discuss 
general features of effective field theories. We present an intuitive introduction to Soft- 



Collinear Effective Theory in Sec. 1.4 



1.2 Predicting Collisions at the LHC 

In this section we introduce the work described in this thesis. Starting with a schematic 
picture of proton-proton collisions, we review the topic of factorization and discuss the 
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Figure 1-1: Schematic picture of a collision at the LHC or Tevatron. It shows the evolution 
in time from left to right, starting with the proton and (anti)proton in the initial state and 
ending with jets, leptons and soft radiation in the final state. 

relevance of experimental restrictions on the hadronic final state. Our work focusses on 
how such constraints can be incorporated into factorization theorems and we give a preview 
of our results. 

The primary goal of the experiments at the LHC and Tevatron is to search for the Higgs 
particle and physics beyond the Standard Model through collisions at the energy frontier. 
The fact that the short-distance processes of interest are interlaced with QCD interactions 
complicates the search. A schematic picture of a proton-proton collision is displayed in 
Fig. |1-1[ A quark or gluon is extracted from each proton (the red circles labeled /), and 
emits strong initial-state radiation (X) prior to the hard short-distance collision (at H). 
The hard collision produces strongly interacting partons which hadronize into collimated 
jets of hadrons (Ji,2,3)) as well as non-strongly interacting particles (represented in the 
figure by the i'^l^). Finally, all the strongly interacting particles, including the spectators 
in the proton, interact with low-momentum soft gluons and can exchange perpendicular 
momentum by Glauber gluons (both indicated by the short orange lines labeled S). 

Factorization is one of the most fundamental concepts in the theoretical description of 
these collisions. It is the statement that the cross section can be computed as a product 
of probability functions, each describing a part of Fig. |1-1[ For a review of factorization 
see Ref. |77| . for examples of factorization theorems see Sec. |2.1[ Factorization allows one 
to separate the new hard physics from the array of QCD interactions in the initial and 
final states and is therefore key in the search for new physics. It is also necessary for 
controlling QCD effects. For example, the momentum distributions of the colliding partons 
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in the protons are nonperturbative, but factorization can imply that these are described by 
universal distributions which have been measured in earlier experiments. 

Monte Carlo programs are a widely used numerical method and employ notions from 
factorization and properties of QCD in the soft and coUinear limits to model the ingredients 
in the full cross section, da = H f f I I Yl^ Ji S . They have the virtue of 
providing a general tool for any observable, but have the disadvantage of making model- 
dependent assumptions to combine all the ingredients and to calculate some of them. When 
a factorization theorem for a specific observable is available this provides a better approach, 
since the various ingredients are defined and combined in a rigorous way. 

In order to probe the physics of the hard interaction, measurements often impose 
restrictions on the hadronic final state, requiring a certain number of hard leptons or 
jets [H [21 [38l I39j . For example, a typical new physics search looking for missing trans- 
verse energy may also require a minimum number of jets with pT above some threshold. 
To identify the new physics and determine the masses of new-physics particles, one has to 
reconstruct decay chains with a certain number of jets and leptons in the final state. 

However, for the majority of processes of interest at hadron colliders where one distin- 
guishes properties of the hadronic final state, so far no rigorous field-theoretic derivation of 
a factorization theorem to all orders in perturbation theory exists. The most well-known 
factorization theorem is 



da = Y, da,7* ® M^a) M^b) ■ (1.1) 



Here fi and fj are the standard parton distribution functions (PDFs) that describe the 
probability to extract a parton i,j = {g, u, u,d, . . .} out of the proton, with a fraction ^a,b 
of the proton momentum. The hard scattering of i and j is described by the partonic cross 



section daf^^^ calculated in fixed-order perturbation theory. In Eq. (1.1), the hadronic final 



state is treated as fully inclusive. Hence, in the presence of experimental restrictions that 



make a process less inclusive, Eq. (1.1) is a priori not applicable. 



Factorization theorems for processes near threshold are a well-studied case where Eq. (1.1) 



can be extended to sum large phase-space logarithms |l2l [53l |63l [66l I118| I12H 11321 1165| . 



(We discuss large logarithms in Sec. 1.3 ) However, threshold production requires the limit 



^a,b 1- The PDFs fall off steeply in this limit, so this is not as relevant for the LHC |62j . 
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Figure 1-2: Examples of events at the LHC. The denote the incoming protons. Jet a 
and b are due to initial-state radiation. 



Our goal is to study factorization for a situation where the hard interaction occurs 
between partons with momentum fractions away from the limit ^^,6 ~^ 1) and where the 
hadronic final state is measured and restricted by constraints on certain event shapes. 
These restrictions allow one to probe more details about the final state and may be used 
experimentally to control backgrounds. 



A typical event at the LHC with three high-py jets is illustrated in Fig. 1-2 (a) There 
are several complications one has to face when trying to derive a factorization theorem in 
this situation. First, experimentally the number and properties of the final-state jets are 
determined with a jet algorithm. Second, to enhance the ratio of signal over background, 
the experimental analyses have to apply kinematic selection cuts. Third, in addition to the 
jets produced by the hard interaction, there is soft radiation everywhere (which is part of 
what is sometimes called the "underlying event"). Fourth, a (large) fraction of the total 
energy in the final state is deposited near the beam axes at high rapidities. Some of this 
radiation can contribute to measurements, and when it does, it cannot be neglected in the 



factorization. In this thesis we focus on the last three items. In chapter 10 we take some 
steps towards resolving the first item. Methods for including jet algorithms in factorization 
have been studied in Refs. [32l [Ml 1173] 

To explore the implications of restrictions on the hadronic final state, we consider the 



simpler situation shown in Fig. l-2(b) where there are no hard central jets. This is already 
interesting for h — )• WW decaying non-hadronically, since it could help reduce the back- 
ground from top quarks decaying into a W plus 6-jet. We will now construct an observable 
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that we use in our analysis to impose this central jet veto. We define two hemispheres, a 
and b, orthogonal to the beam axis and two unit lightlike vectors and rib along the beam 
axis pointing into each hemisphere. Taking the beam axis along the z direction, hemisphere 
a is defined as z > with Ua = (1, 0, 0, 1), and hemisphere b as z < with = (1, 0, 0, — 1). 
We now divide the total momentum px of the hadronic final state into the contributions 
from particles in each hemisphere, px = PXa ~^ PXf We measure the components and 
Bj^ of the hemisphere momenta defined by 

Bt =na- PXa , B^ = nh- px, ■ (1-2) 

Because of the dot product with Ua or Ub, energetic particles near the beam axes only give 
small contributions to B^ or Bj^ . In particular, any contributions from particles at very 
large rapidities outside the detector reach, including the remnant of unscattered partons 
in the proton, are negligible. Demanding that i?^^ are small only allows highly energetic 
particles inside jets along the beam directions labeled "Jet a" and "Jet b" in Fig. |l-2(b)[ 



Hence, measuring and constraining i?^^ provides a theoretically clean method to control 
the remaining particles in the hadronic final state. 

In this thesis, we will prove a factorization theorem for Drell-Yan pp — )• X£^i^ where 
X is allowed to have hard jets close to the beam but no hard central jets, corresponding 
to Fig. l-2(b)[ We call this "isolated Drell-Yan". Our proof of factorization uses Soft- 



Collinear Effective Theory (SCET) [291 Ell ESI EZ] plus additional arguments to rule out 
effects of Glauber gluons based in part on Refs. |22| [78] . Our factorization theorem applies 
to processes pp — )• XL, were the lepton pair is replaced by other non-strongly interacting 
particles, such as Higgs or Z' decaying non-hadronically. This is particularly interesting for 
H —7- W^W~ — )■ I'^vt.'v where a central jet veto is required to remove the large background 
from top quarks decaying into a W plus a 6-jet. We will also write down a factorization 
formula with beam functions for processes with hard central jets. This is not as rigorous as 
the 0-jet case, in particular we have not yet shown the cancellation of Glauber gluons. 

Our main result is to show that process-independent "beam functions", Bi{t,x) with 
i = {g,u,u,d, . . .}, are required to properly describe the initial state. Generically, by 
restricting X one performs an indirect measurement of the proton prior to the hard collision. 
At this point, the proton is resolved into a colliding hard parton inside a cloud of collinear 
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and soft radiation. The proper description of this initial-state jet is given by a beam function 
in conjunction with a soft function describing the soft radiation in the event. 

One might worry that the cohision of partons inside initial-state jets rather than partons 
inside protons could drastically change the physical picture. The changes are not dramatic 
but there are important implications. The beam functions can be computed in an operator 
product expansion, giving 



Bi{t,tfiB) = S{t)fi{C,^^B) + 0[as{fiB)], (1.3) 

where /i^ is an intermediate perturbative scale and t is an invariant-mass variable closely 
related to the off-shellness of the colliding parton (and the Mandelstam variable t). Thus, 
the beam functions reduce to standard PDFs at leading order. For what we call the gluon 
beam function, this was already found in Ref. [96J , where the same matrix element of gluon 
fields appeared in their computation of 7p — ?• J/iJjX using SCET. 



Equation (1.3) implies that the momentum fractions ^a,b are determined by PDFs eval- 
uated at the scale Qi which is parametrically smaller than the scale Q of the partonic 
hard interaction. The renormalization group evolution (RGE) for the initial state now 
proceeds in two stages. For scales fi < fiB, the RGE is given by the standard PDF evo- 
lution [9l [89l 11041 1107^ I109j , which redistributes the momentum fractions in the proton to 
lower ^ values and mixes the gluon and quark PDFs. For scales > fiB, the jet-like struc- 
ture of the initial state becomes relevant and its evolution is properly described by the RGE 
of the beam function. In contrast to the PDF, the evolution of the beam function does not 
involve mixing between quarks and gluons and only changes t. In addition to the change in 
evolution, the transition from PDFs to beam functions at the scale fis also involves explicit 



dsifJ-s) corrections as indicated in Eq. (1.3). These include mixing effects, such as a gluon 
from the proton pair-producing a quark that goes on to initiate the hard interaction and an 
antiquark that is radiated into the final state. For our observables such fluctuations are not 
fully accounted for by the PDF evolution. These beam effects must be taken into account, 
which can be done by perturbative calculations. Fortunately, the standard PDFs are still 
sufficient to describe the nonperturbative information required for the initial state. 

One should ask whether the description of the initial state by beam functions, as well 
as their interplay with the soft radiation, are properly captured by current Monte Carlo 
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programs used to simulate events at the LHC and Tevatron, such as Pythia \159\ I160| and 
Herwig |23| 180] . In these event generators the corresponding effects should be described 
at leading order by the initial-state parton shower in conjunction with models for the un- 
derlying event |24| l60l I16H I162j . We will see that the initial-state parton shower is in 
fact closer to factorization with beam functions than to the inclusive factorization formula 



in Eq. (1.1). In particular, the physical picture of off-shell partons that arises from the 
factorization with beam functions has a nice correspondence with the picture adopted for 
initial-state parton showers a long time ago |45| I158j . Experimentally, measurements of 
the isolated Drell-Yan cross section provide a simple obsevable that can rigorously test the 
accuracy of the initial-state shower in Monte Carlo programs, by contrasting it with our 
analytic results. 



1.3 Weak Interactions and Effective Field Theory 

Many physical processes involve several scales and perturbative calculations typically lead 
to logarithms of the ratios of these scales, e.g. In ni/ fi2- If the scales are widely separated 
/^i S> fJ.2 the logarithms become large and perturbation theory breaks down. We will see in 
this section how effective field theories conveniently solve this problem by factoring a process 
into separate pieces, each corresponding to only one of the scales. For an introduction to 
effective field theories see for example Refs. [102^ 11371 11411 [T56] . 

As an example we will study the decay of a Z? meson to a and a tt^ meson. A 
review of effective field theories for weak decays can be found in Refs. |58t I59j . At the 
parton level we can understand this process as c — t- sud along with a spectator u quark, 
which is mediated hy a W boson of the weak force. The tree level diagram is shown in 



Fig. 1-3 (a) and given by 



(^) V;,Kd ,2 ' 2 (S^" - ^) UsI^PlUc UuluPlUd . (1.4) 

Here the g2 is the coupling of the weak interactions, Vij are elements of the CKM matrix 
from the weak interaction and Uj is the spinor for an external quark of flavor i. The quarks 
are confined to the mesons, implying that the momentum ^ myY- We can therefore 
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(a) (b) 

Figure 1-3: The tree-level process in the full theory (a) and effective theory (b). 



expand the propagator as 



A;2 - mlr (' 



m 



w 



m 



w 



m 



w 



(1.5) 



This removes the W boson as a dynamical degree of freedom and replaces the tree-level 



process in Fig. l-3(a) by the four-fermion interaction in Fig. l-3(b) The W propagator 



gets shrunk to a point because we no longer resolve the short distance physics. 

The dominant corrections to this tree-level diagram come from the strong force, which 
we now consider. The systematic approach to include these corrections in the Lagrangian 
of the effective theory is called "matching" and is based on the operator product expansion 
(OPE) |175j . First we write down all interactions in the low energy theory that are allowed 
by symmetries. The matching will occur at a scale fi ~ so we can treat the c, s, d 

and u quarks as massless and use the corresponding chiral symmetry to constrain possible 
terms. Using Fierz identities for spin and color, we are left with just two terms in our 
effective Lagrangian 



v:,Vud [Cl{^l)Oli^,) + ^2(^^)02(^)1 , 



(1.6) 



where Gp = V^92 / {^''^w) Fermi constant and j and k are color indices. Since we are 

working at leading order in k^/mw, we do not include terms with derivatives that one would 



get from higher order terms of the expansion in Eq. (1.5). To determine the coefficients 



Ci(//) we calculate some physical quantity in both the full theory and the effective theory 
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(a) (b) (c) (d) (e) (f) 



Figure 1-4: The one-loop QCD corrections in the full theory (a)-(c) and effective theory 
(d)-(f). These diagrams have symmetric counterparts. 



and match. A common choice is renormalized matrix elements, where we are free to choose 
any on-shell or off-shell states since the effective theory should completely reproduce the 



low energy limit of the full theory. From our tree-level calculation in Eqs. (1.4) and p. 5 



we know that C| (n) = 1 and C2 {p) = (the superscript (n) denotes terms of order a"). 
Moving on to the Os corrections, the one-loop diagrams that contribute to the matrix 



element with external free quarks are shown in Fig. 1-4 for the full theory and for the 
effective theory. The ultraviolet (UV) divergences we encounter in these calculations are 
regularized using dimensional regularization (DR) and renormalized in the modified minimal 
subtraction scheme (MS). Performing the matching onto the effective theory at order as 

{sud\id^^\cy^^ = {sud\iC^^\c) 



-effl \(1) 

i=l,2 



results in 



Ci(^) = l + A^ln"^^ 



C2(M) = -3"#ln"t. (1. 



Here A^c = 3 is the number of colors in QCD. It is worth pointing out that the full theory 
and effective theory results depend on the momenta of the external quarks (which we 
took to all be equal), but that this dependence cancels in the matching. This is an explicit 
check that effective theory reproduces the full theory in the low energy limit. 
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We now take a closer look at the full theory calculation. Schematically it has the 
following structure 



{sud\iC^;^\c) = 1 + ttsL + a;L^ + a^L^ + ... 



as + a^L + a^L^ + . . ■ 

+ a^L + ... , (1.9) 

where 1 denotes the tree level result and L = ln(— m^/p^). For our D meson decay the 
momentum <^ fnw, which implies that L is large and logarithmic terms such as UsL are 
more important than the finite (non- logarithmic) Us corrections. Assuming that L ~ l/c^s, 
the terms on the first row would all be of 0(1) and thus need to be summed to get the 
correct treelevel amplitude. This is called leading logarithmic (LL) resummation. The 
finite as correction (first term on the second row) is of the same order as the logarithmic 
terms a^'^^L"' on the second row. Summing these terms corresponds to next-to- leading 
logarithmic (NLL) resummation. We will see how this resummation can be achieved in the 
effective field theory. 

Large logarithms are a generic feature of theories with disparate scales. Our effective 
field theory approach allows us to separate these scales and the corresponding physics: 

2 

{sud\iC'^'\c) = 1 + a. In ^ + . . . 

2 2 

1 + a. In ^ + . . . ) (l + a. In ^ + . . . ) = {sud\iC^^\c) . (1.10) 



^cud\Oi{ii) \ b) 



The physics associated with mw can only enter in the Wilson coefficients. The low energy 
physics remains as dynamical degrees of freedom in the effective theory and drops out in the 
matching. The large logarithms can now be avoided by evaluating the matching coefficients 



in Eq. (1.8) at a scale fx ~ mw and the matrix elements at a scale ~ vP^- 

We can use the renormalization group evolution (RGE) |61l IIOH \nO\ \171\ I174j to run 
the Wilson coefficients Ci{^) from a scale ~ my/ down to /u ~ The renormalization 

scale can roughly be interpreted as the scale at which the UV divergences are cut off. 
Lowering ^ thus corresponds to reducing the degrees of freedom in the effective theory and 
absorbing the corresponding change in the Wilson coefficients Cj. With this picture in mind 
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matching 


running 
anom. dim. 


r 

/3(a,) 


LO 


0-loop 


_ 




NLO 


1-loop 






NNLO 


2-loop 






LL 


0-loop 


1-loop 


1-loop 


NLL 


1-loop 


2-loop 


2-loop 


NNLL 


2-loop 


3-loop 


3-loop 



Table 1.1: Order counting for fixed and resummed perturbation theory with single loga- 
rithms. 



it makes sense that we should match at ;U ~ m\Y and then run down to /x ~ The 
full theory calculation contained no UV divergences in the final answer, because the mass 
of the W acted as a UV cutoff. In our effective theory we are expanding about the limit 
mw — )• oo and hence In — m^/p^ leads to UV divergences, which is what allows us to factor 
and resum these large logarithms. 

We will now see explicitly how the RGE resums the large logarithms in our example. 
The UV divergences from the one- loop effective theory diagrams in Fig. |l-4| ^d)-(f) lead to 
the following RGE 



^ C'i (/^) = Hj (m) Cj (/^) , lij ilJ) 



asip) I ^/Nc -3 
-3 3/iVc 



2tt 



(1.11) 



This can be solved by diagonalizing 



C± = Ci±C2, C±if,) = C±imw) 



7±/(2/3o) Tvr -r 1 

. = (1.12) 



where /3o = (HC'yi — ^Tpnf)/?) is the lowest order coefficient in the running of as- As 
a consequence, even though C2{mw) = at tree level a non-zero C2(/u) gets generated 
through the RGE. By expanding we explicitly see how the RGE sums the LL series 



as{mw) 



7±/(2/3o) 



exp 



'y±— — In 



47r 



~ l + asL + a'iL'^ + ... . 



(1.13) 



Resummed perturbation theory is powerful because it captures the most important pieces 
of higher loop diagrams, without having to calculate them. The above resummed result 
contains for example a^L^ which in a full theory calculation would only show up at two 
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loops. The appropriate resummation when we match at order as is NLL resummation, 
for which the running is determined from the UV divergences of two-loop diagrams. The 



After matching and running, the final step to calculate the D — )• K~^tt decay consists 



decays Heavy Quark Effective Theory (HQET) ^TJ 11031 1108j allows you to exploit some 
further symmetries, but at the end of the day this is a non-perturbative matrix element that 
needs to be determined using either lattice QCD, modeling or experiments. The good news 
is that in effective field theories the non-perturbative matrix elements are often universal: 
they can be extracted from some experiment and used in a different one. The PDFs that 
we encountered in the previous section are an example of this. 

The weak decay we have studied captures the generic situation in effective field theories: 
a large separation of scales allows us to expand onto an effective theory but also leads to 
large logarithms. These large logarithms are factored in the effective theory. The Wilson 
coefficients Cj(^) contain the effects of physics at the high scale and the operators Oj(^) 
describe the dynamics of the low energy physics. The renormalization scale /x can be thought 
of as the scale that separates what goes into Ci{fi) and what is described by Oi{fj,). The 
large logarithms are resummed by calculating Ci{fi) at the high scale ("matching") and 
using the RGE to run it down to the low scale ("running"). The relevant matrix elements 
of the operator Oi{fi) at the low scale may be non-perturbative, in which case they should 
be extracted from experiment or may be determined using lattice QCD. 

We conclude this section with a somewhat technical point related to our choice of renor- 
malization scheme. For this we will look at the running of the electromagnetic coupling in 
MS and the off-shell momentum subtraction (ms) scheme. The one-loop contribution from 
the electron with mass m leads to 



running is generally at one higher loop that the matching, as is shown in Table |1.1[ 



of determining the matrix elements 





(1.14) 



For /i > m we have ^""(/n) = (fi), but 



e3(/i) ^2 



for /i ^ m . 



(1.15) 
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So below the electron mass the effect of the electron decouples from the runnmg, with a 
smooth transition through fi = m. By contrast, the electron doesn't decouple in the MS 
scheme {P^^ doesn't depend on m at all). As this example illustrates, adopting a so- 
called physical renormalization scheme causes high energy degrees of freedom to decouple 
naturally. However, the expressions become much more involved (compare /S™*^ with f3^^), 
which is why we prefer to use the MS scheme and manually perform the decoupling by 
matching and running. 



1.4 Soft-Collinear Effective Theory 

In a typical event at the LHC the final state contains photons, leptons and hadrons (see 



Fig. l-2(a)). The energetic hadrons generically come in collimated jets and there are also 
much less energetic hadrons everywhere. These two types are referred to as "coUinear" and 
"soft" radiation and are enhanced due to the presence of collinear and soft infrared(IR) 
singularities in QCD. SCET is an effective theory of QCD describing these modes |29l 
[3H [35l [37] . In this section we present an intuitive introduction to SCET. We start by 
introducing light-cone coordinates and describing the degrees of freedom in SCET. We will 
discuss the matching onto SCET, reparametrization invariance (RPI), the factorization of 
soft and collinear degrees of freedom and the resummation of large logarithms. 

It is convenient to use light-cone coordinates to describe jets of energetic hadrons. 
Choose a light-cone vector n'^ vector pointing in the direction of one of the energetic jets 
and an additional light-cone vector such that n ■ n = 2. For n^^ = (1,^) a common 
choice is = (1, —n), where of course = 1. We can then decompose any four- vector 
ptM — (p+,p^,p^) as 

= p'^ — + — + p'^ , p^=n-p, p~=n-p. (1-16) 
The momenta of energetic hadrons in the jet along the direction are parametrically 

p" = ip'^,p',p'i) ^ QiX^l,\) (collinear), (1.17) 
where Q is the scale of the hard collision and X ^ p±/p^ ^ 1 characterizes the size of the 



jet. The power counting in Eq. (1.17) corresponds to a particle with virtuality p^ ~ X^Q^ 
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H 





Figure 1-5: Schematic diagram of a proton- (anti) proton collision up to the hard interaction. 
A quark or gluon is taken out of each of the protons / and propagates while emitting initial- 
state radiation I until it enters the hard interaction H. The two shades of green correspond 
to the two coUinear modes. For our observables the proton remnant is much more collinear 
and goes straight down the beam pipe. 



that is boosted in the n'^ direction (we consider light particles whose mass we can ignore). 



For the isolated Drell-Yan process shown in Fig. l-2(b) , we have two jets from initial- 
state radiation corresponding to the directions Ua = (1,0,0,1) and = (1,0,0,-1). We 
thus have two collinear modes 

Pa = {pt,Pa:Pa±) ~ Qi^^^ 1' ^) (^a collinear) , 

Pfe = (P^T^Pj^^Pfcx) ~ A) (nf, collinear) , (1.18) 

where we take = and = fif, for convenience. By convention the will always 
correspond to the large light-cone component. Each of the two initial-state jets originates 
from a single quark or gluon coming out of the proton, which propagates and radiates until 



it enters the hard interaction, as shown schematically in Fig. 1-5 

Since the final state hadrons are colorless, there must be some degrees of freedom con- 
necting the jets that carry color charge. The collinear modes cannot interact directly because 
that would lead to a momentum 

P^^+K~Q(A',l,A) + g(l,A2,A) = Q(l,l,A), {pa+Pbf^Q^. (1.19) 

In SCET we have integrated out the hard interaction and so the corresponding modes 
with virtuality are no longer dynamical degrees of freedom. The interactions between 
the collinear modes is provided by additional "soft" degrees of freedom which have the 
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^2 



Figure 1-6: Feynman rule for the interaction between an n-collinear quark (dashed hne) 
and a soft ghion, at leading order in the power counting. 



following power counting^ 

= {pt,Ps,Ps±) ~ A^, (soft) . (1.20) 

They describe much less energetic radiation E ^ }?Q without a preferred direction and 
correspond to long range fluctuations since their virtuality ~ QA^ is smaller. At a 
more technical level we expand the collinear momenta = -\- into a discrete large 
momentum p^ = {0,p~,pj_) ~ Q(0,1,A) and a small residual momentum pr ~ QX^ that 
can be exchanged with the soft radiation. The interaction between a soft and a collinear 



mode at leading order in the power counting is shown in Fig. 1-6 

SCET differs in several ways from effective theories where a heavy particle has been 
integrated out (such as the W boson in the previous section). The expansion parameter 
A is not set by the mass of a heavy particle and there are different modes corresponding 
to the same particles, e.g. we have both collinear gluons and soft gluons. We thus have to 
be careful not to double count degrees of freedom, which can be removed through so-called 
zero-bin subtractions |140j . 

One might wonder if there are any other degrees of freedom that we should include in 
our effective theory. The answer to this question depends on the observable you want to 
study. For the weak decay in the previous section the small external momenta told us that 
we did not need to treat the W boson as a dynamical degree of freedom in our effective 
theory. For the isolated Drell-Yan process the invariant mass of the leptons sets the scale 
of the hard interaction Q. Measuring the and Bj^ variables introduced in Eq. (1.2) 



and restricting B^^^ <^ Q implies that energetic radiation is restricted to be along the beam 



axis. This is what leads to the two collinear modes in Eq. (1.18) with A ~ B'^^/Q. The 



'^In the literature these are sometimes referred to as ultrasoft modes, where soft is reserved for modes 



with power counting (p ,p ~ Q(A, A, A). As discussed below Eq. (1.21|, we do not need these modes. 

30 



contribution to and B'^ from the soft radiation is of the same parametric size as that 
of the colUnear radiation, so we need to include the soft degrees of freedom. 
In principle you could imagine including modes with a scaling such as 

(a) p~Q(A,A,A), 

(6) p~Q(A^A^A^), 

(c) p~Q(A^l,A2), 

(d) J9~Q(A^A^A). (1.21) 

In (a) we are considering a mode with virtuality p"^ ~ X^Q^ that, like the soft radiation, has 
no preferred direction. This mode would lead to contributions of order \Q to B^^^. However, 
our measurement demands i?^^ ~ X^Q and so (a) is ruled out in final state radiation. We 
now address virtual contributions from (a), by arguing that these modes cannot couple to 
the collinear and soft degrees of freedom. Attempting to couple (a) to collinear modes leads 
to a virtuality p'^ ~ AQ^, which is too large and has been integrated out. The momentum 



components of (a) are all larger than those of the soft modes in Eq. (1.20), and so these 
short range fluctuations are not resolved by our soft modes. There is no coupling between 
(a) and the soft degrees of freedom in the leading order Lagrangian. 

The degrees of freedom in (b) correspond to a soft mode with virtuality p"^ ~ X^Q'^ and 
is already contained in our soft degrees of freedom. Alternatively, the contribution of (b) to 
B^ and Bj^ is power suppressed, so there is no need to separately include them. The same 



is true for (c), which is a collinear mode with a smaller virtuality than in Eq. (1.17) and is 
therefore already contained in our existing collinear degrees of freedom. 

Gluons with a scaling as in (d) are called Glauber gluons. They cannot enter in the 
final state since they are offshell modes, but they can mediate perpendicular momentum 



transfer between the two different collinear sectors in Eq. (1.18). In particular they could 
connect the active parton from one proton with a quark or gluon in the remnant of the other 
proton. As a consequence the proton remnant no longer goes straight down the beam pipe 
and contributes to B'^^^, see Fig. 



1-7 



We will discuss below how to factorize the soft and 
the various collinear degrees of freedom. If Glauber gluons need to be taken into account, 
this factorization breaks down. In Sec. 17.3.41 we show that the contribution from Glauber 
gluons cancels out for our factorization theorem. 
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Figure 1-7: Schematic diagram of a proton- (anti)proton collision up to the hard interaction. 
Glauber gluons (double line) can couple an active parton with a spectator from the other 
proton. The spectator was very collinear and would have gone straight down the beam pipe 
but the Glauber gluon kicks it into active region, where it contributes to -B^^^. 



To summarize, identifying the relevant modes is really the art of effective field theories. 
Including any unnecessary degrees of freedom is permitted and will not affect the outcome, 
as long as one avoids double counting by appropriate zero-bin subtractions. However, it 
will dramatically increase the computational effort required. 

We now discuss the matching onto SCET in which we integrate out the hard inter- 
action [3D]. In the previous section our effective Lagrangian only contained two terms 

~ J2i=i 2 CiOi, because the weak interaction is insensitive to the details of our external 
states. In SCET the collinear degrees of freedom have large p~ ^ Q components, which 
enter in the hard interaction. Matching onto SCET thus yields Wilson coefficients and 
operators that depend on both the direction and size p~ of these large components, 

^eff^^^ y y dp^Ci(na,p;;;,n5,p^)C'i(na,Pa ,nfc,p^) . (1.22) 

i na,ni, 



As an example we show the matching for Higgs production through gluon fusion in Fig. 1- 
The collinear directions are fixed by the directions of the protons and p^^^ are related to the 
virtuality and rapidity of the Higgs by momentum conservation. 

Collinear gluons have the scaling An = {A^ , A~ , A'^j_) ~ (A^,1,A), which e.g. follows 
from the homogenous power counting of the covariant derivative. This implies that we can 
add arbitrary many collinear n ■ An fields in the matching at the same order in the power 
counting, see Fig. |l-8[ c) and (d). Fortunately gauge invariance restricts these operators, 
organizing these gluons into collinear Wilson lines Wn |37] . 
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(a) (b) 




Figure 1-8: The matching for Higgs production through gluon fusion. The full theory 
diagram (a) gets matched onto the effective interaction (b). Springs with a line through 
them are collinear gluons. The different shades of green correspond to the two collinear 
modes and the solid black lines with the hard degrees of freedom that get integrated out. 
In (c) we consider the emission of an additional n-A which is at the same order in the power 
counting and gets matched onto (d). Including arbitrary many emissions of n-^ gluons 
generates a collinear Wilson line Wn in SCET. 



There was some arbitrariness in our choice of and n^, since the light-cone vector 
only needs to point roughly in the direction of the jet and we can choose any that 
satisfies n-n = 2. The precise choice of these parameters should not affect our calculations, 
providing us with additional symmetries known as reparametrization invariance (RPI) [139j . 
We will often use RPI-III that transforms — )• e"n'^ and — )• e~°'n^. As an example, an 
RPI-III invariant function f{p'^,p^) can only depend on the invariant combination p^p~ . 

We now discuss the factorization of the soft and various collinear degrees of freedom. It 
turns out that at leading order in the power counting the coupling of soft degrees of freedom 
to any collinear sector organizes itself into a Wilson line Y of soft gluons [35j , as shown in 



Fig. 1-9 This is due to the simple coupling between soft and collinear degrees of freedom 



in Fig. 1-6 and their relative power counting in Eqs. (1.17) and (1.20). By factoring out 



this soft Wilson line from the collinear quark and gluon field and An 



AU^) = Y{x)A':^'\x)Y\x), 



(1.23) 
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(a) (b) 

Figure 1-9: The soft radiation in (a) organizes itself into emissions from soft Wilson lines 
shown in (b). The dashed lines correspond to collinear quarks, and the springs with (with- 
out) a solid line through them correspond to collinear (soft) gluons. 

we get collinear fields \ An^^^ that are decoupled from the soft degrees of freedom at 
leading order in the power counting. This is known as the BPS field redefinition. 

We conclude this section by discussing the resummation of large logarithms. Our mea- 
surement sets the expansion parameter ~ B^/Q <^ 1 and SCET will allow us to resum 
the corresponding large logarithms L = In B^/Q. However, rather than single logarithms 
UsL we now encounter double logarithms a^L^ which leads to a different structure than in 



Eq. (1.9) 



cr ~ 1 + asL^ + a^L^ + a^L^ + ... 

asL + aj.L^ + a^L^ + . . . 
as + alL^ + a^L^ + ... . (1.24) 

If we assume L ~ l/os, getting the correct tree level result requires us to sum the first 
row, the second row, everything on the third row except for the first term etc. This looks 
complicated, but once these logarithms are factored in SCET the story is very similar to 



that in Sec. 1.3 Each factor only depends on one physical scale and by evaluating it at 



this scale we avoid large logarithms. We then use the RGE to evolve all the factors to a 



common scale fi, which resums the logarithms in Eq. (1.24). The evolution of a factor has 
the following parametric structure 

?7(Aio,/i) ~exp [^(a,L)^L + (a,L)*^ + a,(a,L)*^ + ...] . (1.25) 

A; 



The first term, Y^^{asL) L, was not present for the weak decays in Eq. (1.13). It is due 



to the cusp anomalous dimension which comes from the cusp of Wilson lines. This term is 
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matching 


non-cusp 


running 
cusp 




LO 


0-loop 








NLO 


1-loop 








NNLO 


2-loop 








LL 


0-loop 


0-loop 


1-loop 


1-loop 


NLL 


0-loop 


1-loop 


2-loop 


2-loop 


NNLL 


1-loop 


2-loop 


3-loop 


3-loop 



Table 1.2: Order counting for fixed and resummed perturbation theory with double loga- 
rithms. 



related to the different structure of the logarithms in Eq. (1.24), as is clear from expanding 
U{fiQ,fi). Assuming L ~ l/a^, we need the first two terms of the exponent in Eq. (1.25) 
to get the correct tree level result, which requires the two- loop cusp anomalous dimension. 
The non-cusp piece of the anomalous dimension does not contribute to the first term of 



the exponent in Eq. (1.25) and is only needed at one-loop. The general counting for fixed 



order and resummed perturbation theory is shown in Table 1.2 We will not consider LL 
resummation because this does not include the '^^i'^sL)^ term in the exponent in Eq. ( 1.25 ), 
which parametrically contributes at 0(1). 



1.5 Outline 

Chapter [2] provides an introduction to beam functions and factorization for isolated pro- 
cesses, in which we discuss our main results. The goal of this chapter is to give a thorough 
discussion of the physical picture behind our results which is accessible to non-expert read- 
ers. We define the "beam thrust" event shape that we use to impose the central vet jeto in 
isolated Drell-Yan, and we present the factorization theorem. The renormalization of the 
beam functions and the relationship between the beam functions and the PDFs is discussed, 
and we make the comparison with the initial-state parton shower. We compare our factor- 
ization theorem for isolated Drell-Yan to the fixed-order calculation and discuss its RGE 
structure. For convenience we have included a quick overview of our notation and symbols 
in App. |A} In App. [B] we list the plus distribution definitions and identities that we use, as 
well as identities for discontinuties needed in the calculations in Ch. |4] and Ch. [5l 

In chapter [3] we discuss several formal aspects of the quark and gluon beam functions, 
including their definition in terms of matrix elements of operators in SCET, their renor- 
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malization, their analytic structure, and the operator product expansion relating the beam 
functions to PDFs. In particular, we show that the beam functions have the same RGE 
as the jet functions to all orders in perturbation theory. (Part of the proof is relegated to 
App.[C}) 

We perform the one-loop matching of the quark beam function onto quark and gluon 
PDFs in Ch. |4} Using an offshellness IR regulator we give explicit details of the calcula- 
tions for the quark beam function and PDF. We verify explicitly that the IR singularities 
cancel in the matching and extract results for the next-to-leading order (NLO) matching 
coefficients. In App. [D] we repeat the matching calculation for the quark beam function 
in pure dimensional regularization. The calculation of the one-loop matching coefficients 
for the gluon beam function is performed in Ch. [5] In chapter [6j we show and discuss 
plots for the quark and gluon beam function at NLO in fixed-order perturbation theory and 
next-to-next-to-leading logarithmic (NNLL) order in resummed perturbation theory. 

In chapter[7j we derive in detail the factorization theorem for isolated processes pp — t- XL 
using SCET, and apply it to the case of Drell-Yan. 

Plots for the beam thrust cross section for Drell-Yan at NLO and NNLL are shown 
and discussed in Ch. |8j In Ch. [9] we apply our factorization theorem to Higgs production 
through gluon fusion and show results for the corresponding beam thrust cross section. The 
necessary ingredients for evaluating these cross sections at NNLL are collected in App. [El 



In Ch. 10 we extend our work to final state jets. For a signal with N jets we define a 
global event shape called "N-jettiness" to veto any unwanted additional jets. We describe 
theoretical and experimental benefits of our approach and present a factorization formula 



for the corresponding A'^-jet cross section. We conclude in Ch. 11 
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Chapter 2 



Introduction to Beam Functions 
and Isolated Drell-Yan 



This chapter provides an extensive discussion of how factorization with beam functions 
works, including the necessary kinematic definitions for the variables that constrain the 
hadronic final state. In the interest of avoiding technical details, we only discuss the physics 
contained in the factorization theorems. Readers interested in the field-theoretic definitions 
for the beam functions are referred to Ch. |3] while those interested in the derivation of the 
factorization theorem in SCET and explicit definitions for all its ingredients are referred to 
Ch. [7j The work in this chapter was first presented in Ref. |166| . 



In Sec. 2.1, we review the factorization theorems for inclusive Drell-Yan and threshold 
Drell-Yan, and then explain the factorization theorem for our isolated Drell-Yan process. 
We use a simple setup where measurements on the final-state hadrons use hemispheres 



orthogonal to the beam. These observables are generalized in Sec. 2.2 to uniformly account 
for measurements that sample over a wide variety of boosts between the hadronic and 
partonic center-of-mass frames. We explain the relation between beam functions and parton 



distribution functions in Sec. 2.3 We compare the beam- function renormalization group 



evolution to initial-state parton showers in Sec. 2.4 In Sec. 2.5, we show how the various 



pieces in the factorization theorem arise from the point of view of a fixed-order calculation. 



In Sec. 2.6 we compare the structure of large logarithms and their resummation for the 
different factorization theorems. This yields an independent argument for the necessity of 
beam functions. 
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2.1 Drell-Yan Factorization Theorems 



To describe the Drell-Yan process pp — )■ Xl^i or pp — )■ Xi^i , we take 



(2.1) 



where P^^ are the incoming (anti)proton momenta, Ecm = \/ {Pa + Pb)^ is the total center- 
of-mass energy, and is the total momentum of the i^i" pair. We also define 



' 



y = iin^, 

2 Pa-q 



Xb 



(2.2) 



where y is the total rapidity of the leptons with respect to the beam axis, and Xa and Xb 
are in one-to-one correspondence with r and Y. Their kinematic limits are 



< r < 1 



2|y|<-lnr, T<Xa<l: 



T < Xfo < 1 



(2.3) 



The invariant mass of the hadronic final state is bounded by 



(2.4) 



In Drell-Yan 



Q = > Aqcd 



(2.5) 



plays the role of the hard interaction scale. In general, for factorization to be valid at some 
leading level of approximation with a perturbative computation of the hard scattering, the 
measured observable must be infrared safe and insensitive to the details of the hadronic 
final state. 



For inclusive Drell-Yan, illustrated in Fig. 2-1 (a) , one sums over all hadronic final states 



X allowed by Eq. (2.4) without imposing any cuts. Hence, the measurement is insensitive to 



any details of X because one sums over all possibilities. In this situation there is a rigorous 
derivation of the classic factorization theorem |52^ [75l [78] 



1 do- 



(To dgMy 



, (2.6) 
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Jet 2 Soft Jet 2 

(d) Dijet near threshold. (e) Isolated dijet production. 



Figure 2-1: Different final-state configurations for pp collisions. The top row corresponds to 
Drell-Yan factorization theorems for the (a) inclusive, (b) threshold, and (c) isolated cases. 
The bottom row shows the corresponding pictures with the lepton pair replaced by dijets. 



where ctq = ^'^ot'eni/ {"^-^cEl^^q^), and the integration limits are < < 1 and a^f, < < 1. 
The sum is over partons i, j = {g, n, u,d, . . .}, and fi{S,a) is the parton distribution function 
for finding parton i inside the proton with light-cone momentum fraction along the proton 
direction. Note that (^a,b are partonic variables, whereas Xa^b are leptonic, and the two are 
only equal at tree level. The inclusive hard function -ff™'^^ can be computed in fixed-order 
perturbative QCD as the partonic cross section to scatter partons i and j [corresponding 

part 



to dcrP^" in Eq. and is known to two loops [S O ESI [HOl fTT3]. 



For threshold Drell-Yan, one imposes strong restrictions to only allow soft hadronic final 



states with mx ^ Q, as illustrated in Fig. 2-l(b), Using Eq. (2.4), this can be ensured 
by forcing (1 — \/t)^ ^ ''") so that one is close to the threshold r — t- 1. In this case, 
there are large double logarithms that are not accounted for by the parton distributions. 
Furthermore, since 

1 > L,b > Xa,b > r ^ 1 , (2.7) 
a single parton in each proton carries almost all of the energy, ^a,b 1- The partonic analog 
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of r is the variable 



< 1. 



(2. 



and r — )• 1 implies the partonic threshold limit z — t- 1. As Eq. (2.3) forces y — t- for r — t- 1, 
it is convenient to integrate over Y and consider the r — t- 1 limit for dcj/df;^. The relevant 



factorization theorem in this limit is 1631 1165 



1 dcr 
do dg2 



xQSti,, Qil 



(2.9) 



where we view Eq. (2.9) as a hadronic factorization theorem in its own right, rather than 



simply a refactorization of Hlj^^^ in Eq. (2.6). This Drell-Yan threshold limit has been 



studied extensively 021 SI [Ml HIS IMl EMI ESQ]. Factorization theorems of this type 
are the basis for the resummation of large logarithms in near-threshold situations. In 



contrast to Eq. (2.6), the sum in Eq. (2.9) only includes the dominant qq terms for various 



flavors, ij = {uu,uu,dd, . . .}. Other combinations are power-suppressed and only appear 
at 0(1 — r) or higher. The threshold hard function Hij ~ |CjC'j'| is given by the square 
of Wilson coefficients in SCET, and can be computed from the timelike quark form factor. 
The threshold Drell-Yan soft function Sthr is defined by a matrix element of Wilson lines 
and contains both perturbative and nonperturbative physics. If it is treated purely in 
perturbation theory at the soft scale Q(l — r), there are in principle additional power 



corrections of 0[Aqcd/Q{1 - t)] in Eq. [T25]. 

Our goal is to describe the isolated Drell-Yan process shown in Fig. 2-l(c)[ Here, the 
colliding partons in the hard interaction are far from threshold as in the inclusive case, 
but we impose a constraint that does not allow central jets. Soft radiation still occurs 
everywhere, including the central region. Away from threshold, the hard interaction only 
carries away a fraction of the total energy in the collision. The majority of the remaining 
energy stays near the beam. The colliding partons emit collinear radiation along the beams 
that can be observed in the final state, shown by the green lines labeled "Jet a" and "Jet 



6" in Fig. 2-l(c), This radiation cannot be neglected in the factorization theorem and 



necessitates the beam functions. In the threshold case, these jets are not allowed by the 
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limit r — )■ 1, which forces all available energy into the leptons and leaves only soft hadronic 
radiation.^ In the inclusive case there are no restrictions on additional hard emissions, in 
which case initial-state radiation is included in the partonic cross section in Hf^'^^. 



Also shown in Fig. 2- 1(c) is the fact that the leptons in isolated Drell-Yan need not be 



back to back, though they are still back to back in the transverse plane [see Sec. 7.2 . In 
this regard, isolated Drell-Yan is in-between the threshold case, where the leptons are fully 
back to back with Y ^ 0, and the inclusive case, where they are unrestricted. 



In Figs. 2-1 (d) and 2-1 (e) we show analogs of threshold Drell-Yan and isolated Drell-Yan 



with the leptons replaced by final-state jets. We will discuss the extension to jets in Ch. 10 



To formulate isolated Drell-Yan we must first discuss how to veto hard emissions in 
the central region. For this purpose, it is important to use an observable that covers the 
full phase space. Jet algorithms are good tools to identify jets, but not necessarily to veto 
them. Imagine we use a jet algorithm and require that it does not find any jets in the central 
region. Although this procedure covers the full phase space, the restrictions it imposes on 
the final state depend in detail on the algorithm and its criteria to decide if something 
is considered a jet or not. It is very hard to incorporate such restrictions into explicit 
theoretical calculations, and in particular into a rigorous factorization theorem. Even if 
possible in principle, the resulting beam and soft functions would be very complicated 
objects, and it would be difficult to systematically resum the large logarithms arising at 
higher orders from the phase-space restrictions. Therefore, to achieve the best theoretical 
precision, it is important to implement the central jet veto using an inclusive kinematic 
variable. This allows us to derive a factorization theorem with analytically manageable 
ingredients, which can then be used to sum large phase-space logarithms. For a clean 
theoretical description, this observable must be chosen carefully such that it is infrared safe 
and sensitive to emissions everywhere in phase space. Observables satisfying these criteria 
for hadron colliders have been classified and studied in Refs. \27\ I28| . and are referred to 
as global event shapes. (Issues related to non-global observables have been discussed for 
example in Refs. [201 HHl [H21 EZ]-) 



^Note that the proof of factorization for the partonic cross section in the partonic threshold hmit z ^ 1 
is not sufficient to estabhsh the factorization of the hadronic cross section, unless one takes the limit r — )■ 1. 
The hadronic factorization theorem assumes that all real radiation is soft with only virtual hard radiation in 
the hard function. The weaker limit 2 — >■ 1 still allows the incoming partons to emit energetic real radiation 
that cannot be described by the threshold soft function. Only the r — >■ 1 limit forces the radiation to be 
soft. This point is not related to whether or not the threshold terms happen to dominate numerically away 
from T — >■ 1 due to the shape of the PDFs or other reasons. 
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hemisphere b | hemisphere a 




Figure 2-2: Definition of hemispheres and kinematic variables for isolated Drell-Yan. 



We will consider a simple kinematic variable that fulfills the above criteria, leaving the 
discussion of more sophisticated generalizations to the next section. The key variables for 



the isolated Drell-Yan process are shown in Fig. 2-2, The proton momenta Pa and Pj^ are 
used to define lightlike vectors Ua and n^, 



~2 



Pa'=^<, K=^<^ (2.10) 



where the protons are massless and = 0, = 0, and ria-rif, = 2. Using the beam axis, 
we define two hemispheres a and b opposite to the incoming protons. We then divide up 
the total hadronic momentum as 

p^^ = Bi^ + B^, (2.11) 

where Ba = Px and B^ = p^^ are the total final-state hadronic momenta in hemispheres 
a and b. Of these, we consider the components 

B+ = na-Ba = ^0(1 + tanhy,,) e'^?^" , B+ = Ub-Bb = B^{1 + tanhyb) e-^^' , (2.12) 

where B^^ are the energy components and ^ are the total rapidities of B^^ with respect 
to the forward direction for each hemisphere. Here, liuiy^ooi^ + tanhy) = 2 and 
1 + tanhy > 1.8 for y > 1, so i?^^ scale exponentially with the rapidities ya,b- 
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In terms of the measured particle momenta in hemisphere o, 



Bt = Yl ""--P^ = + tanh%)e-2'?'= . (2.13) 

Here, Ek and r/^ are the experimentahy measured energy and pseudorapidity with respect to 
fia, and we neglect the masses of final-state hadrons. An analogous formula applies for . 
Hence, B'^ and B^ receive large contributions from energetic particles in the central region, 
while contributions from particles in the forward region are suppressed. Thus, requiring 
small B'^^ ^ (5 is an effective way to restrict the energetic radiation in each hemisphere 
as a smooth function of rapidity, allowing forward jets and disallowing central jets. At the 
same time, soft radiation with energies ^ Q is measured, but not tightly constrained. 

As an example, consider the cut 

5+6 < Qe-^J^-t _ (2.14) 

This constraint vetoes any events with a combined energy deposit of more than Q/2 per 
hemisphere in the central rapidity region |y| < ycut- In the smaller region |y| < ycut — 1, the 



energy allowed by Eq. (2.14) is reduced by a factor of e ~ 7, essentially vetoing any jets 
there. In the larger region \y\ < ycut+l, it is increased by the same factor, so beyond ycut + 1 
the hadronic final state is essentially unconstrained. Thus, a typical experimental value 
might be ycut = 2, which vetoes energetic jets in the central region \y\ < 1. The precise 
value of the cut on i?^^ will of course depend on the requirements of the experimental 
analyses. 

Note that the variable B^ is similar to the total transverse energy in hemisphere a, 
defined as 

^Ta = E = E ^^(1 + tanhr?,)e-'"= . (2.15) 

B^ has two advantages over -Eto- First, the exponential sensitivity to rapidity is much 
stronger for B^, which means it provides a stronger restriction on jets in the central region 
and at the same time is less sensitive to jets in the forward region. Second, since B^ is a 
specific four-momentum component and linear in four-momentum, (pi +^2)+ = pf ~^ P2 ^ 
it is much simpler to work with and to incorporate into the factorization theorem. It is 
clear that the isolated Drell-Yan factorization theorem discussed here can be extended to 



43 



observables with other exponents, e , much hke the angularity event shapes in e~^e |47j . 

One should ask, down to what values can B'^j^ be reliably measured experimentally? In 
principle, particles at any rapidity contribute to -B^^, but the detectors only have coverage 



up to a maximum pseudorapidity r/detj as indicated in Fig. 2-2 For the hadron calorimeters 
at the LHC r/^et — 5 and at the Tevatron j^jet — 4. In the hadronic center-of-mass frame, 
the unscattered partons inside the proton have plus components of O {Aq^jj^ / Ecm) , so any 
contributions from the unmeasured proton remnants are always negligible. The question 
then is, what is the maximal contribution to B^^^ from initial-state radiation that is missed 



as it is outside the detector? In the extreme scenario where all proton energy is deposited 
right outside r/det) we would have B^^^ = 14TeVe~^° = 0.6 GeV at the LHC and B^f^ = 
2TeVe~^ = 0.7 GeV at the Tevatron. In more realistic scenarios, the contribution from 
such radiation is suppressed by at least another factor of 10 or more. Therefore, the finite 
detector range is clearly not an issue for measuring values B^^ > 2 GeV, and the relevant 
limitation will be the experimental resolution in B^f^. 

The factorization theorem for isolated Drell-Yan, which we prove in Ch. [7j reads 



c^ dq2dyd^+d^+ = 5^^^^^'^''^) / <iktdk+q''B,[u;a{B+ - kt),Xa,fi] 

AqCD ^afi^tb 



X Bj [u}b{Bl -k^),Xb, fJ^]Sa,^^i{kZ , k^,fi){l + 



Q ' 



(2.16) 



The physical interpretation of Eq. (2.16) is that we take partons i and j out of the initial- 
state jets Bi, Bj and hard-scatter them to final state particles with Hij, while including 
'S'ihemi to describe the accompanying soft radiation. The hard function Hij is identical to 



the one in the threshold factorization theorem in Eq. (2.9), and the sum in Eq. (2.16) is 
again only over ij = {uu,uu,dd, . . .}. The quark and antiquark beam functions Bg and 
Bq describe the effects of the incoming jets and have replaced the PDFs. The variables 
^^a,b = XafiEcm- The hard partons are taken from initial-state jets rather than protons, so 
unlike in the threshold case the gluon PDF now contributes via the beam functions. We will 
see how this works in more detail in Sec. |2.3[ Finally, S'ihemi is the initial-state hemisphere 
soft function. 



The kinematic variables in Eq. ( 2.16 ) are displayed in Fig. 2-2 The soft function depends 



on the momenta k'^ = Ua- ka and k^ = Ub- kj, of soft particles in hemispheres a and b, 
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respectively. Much like PDFs, the beam functions Bi{ta, Xa, and Bj{tf,, X),, n) depend on 
the momentum fractions Xa and Xh of the active partons i and j participating in the hard 
collision. In addition, they depend on invariant-mass variables 



ta = UJabt>0, tb = UJbb+>0, (2.17) 

where ujafi = Xa^bEcm are the hard momentum components and = ria-ba- The momentum 
ba is defined as the total momentum of the energetic particles radiated into hemisphere a, 



as shown in Fig. 2-2, and similarly for 6^. (The kinematics are shown in more detail in 



Fig. 7-1 ) Before the hard interaction, the momentum of the active quark can be written as 



-af -6;:^. (2.18) 

The first term is its hard momentum along the proton direction, and the last two terms 
are from the momentum it lost to radiation, where = —b'^j^ contains the transverse 
components. The quark's spacelike invariant mass is —ujab'^ — b'^rp = —ta — ^^y- The beam 
function Bi for hemisphere a depends on ta = (j-^aba = XaEcmba, which is the negative of the 
quark's transverse virtuality. (When the distinction is unimportant we will usually refer to t 
simply as the quark's virtuality.) By momentum conservation b^ = B^ — , leading to the 



convolution of the beam and soft functions as shown in Eq. (2.16). Physically, the reason 
we have to subtract the soft momentum from B^ is that the beam function only properly 
describes the collinear radiation, while the soft radiation must be described by the soft 
function. An analogous discussion applies to Bj and for hemisphere b. The convolutions 
in the factorization theorem thus encode the cross talk between the soft radiation and 
energetic collinear radiation from the beams. 

By measuring and constraining B^ we essentially measure the virtuality of the hard 
parton in the initial state. As the proton cannot contain partons with virtualities larger 
than Aqqj-,, the initial state at that point must be described as an incoming jet containing 
the hard off-shell parton. This is the reason why beam functions describing these initial- 



state jets must appear in Eq. (2.16). It also follows that since t ^> ^qcd calculate 



the beam functions perturbatively in terms of PDFs, which we discuss further in Sec. 2.3 

It is convenient to consider a cumulant cross section, including all events with B^f^ up 
to some specified value, as in Eq. (2.14). Integrating Eq. (2.16) over < i?^^ < -B^ax 
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obtain 



1 da 



AqCD Wa^foS+ax 



(2.19) 



where the soft function 5ihemi is the same as in Eq. (2.16), and we defined the integrated 
beam function 

/x)= [dtBi{t,X,^i)9{tme..-t). (2.20) 



The cut S+j, < B+^^ imphes the hmit 6+^ < B+^^ - and ta,b < uJa,b{B^ 



k 



leading to the convolutions in Eq. (2.19). 



The factorization theorem Eq. (2.16) and its integrated version Eq. (2.19) are valid in 
the limit ta^b/Q"^ — B'^^/Q = ^ 1, and receive power corrections of 0{\^). Thus, for 
-^max — Qe"^^'^"* with Ucut = 1, we expect the power corrections not to exceed ~ 10%. 



This is not a fundamental limitation, because the power corrections can be computed in 
SCET if necessary. If the soft function is treated purely perturbatively, there are additional 
power corrections of C'(Aqcd/-B^^), which account for soft singularities as i?^^ — >• 0. The 
variables B^^ are infrared safe with respect to collinear splittings |164j . 

The hard function receives perturbative corrections at the hard scale ~ Q, the 
beam functions have corrections at the intermediate beam scale /i| ~ 

^max — QB^aax^ 

and the soft function at the soft scale ns — ^max- For example, for Q ^ 1 TeV and 
ycut = 2 we have /Ub — 140 GeV and jj-s — 20GeV. Even with a very small Q ~ 100 GeV, 
perhaps for Higgs production, //^ ~ 14 GeV and fis — GeV are still perturbative (although 
at this point nonperturbative contributions ~ Aqcd//^S' to the soft function might no 
longer be small and may be incorporated with the methods in Refs. jll7l I135j ). In fixed- 
order perturbation theory, the cross section contains large single and double logarithms, 
ln{B^^^/ Q) ~ —4 and In^ {B^^^/ Q) ~ 16, invalidating a fixed-order perturbative expansion. 
The factorization theorem allows us to systematically resum these logarithms to all orders 



in perturbation theory, which is discussed in more detail in Sec. 2.6 



The factorization theorem Eq. (2.16) also applies to other non-hadronic final states 
such as Z' — )• l^£~, or Higgs production with — )• 77 or i7 — )• ZZ* — )• 4:1. In each 
case, and Y are the total non-hadronic invariant mass and rapidity, and central jets are 



46 



vetoed with a cut on -B^^- The only dependence on the process is in the hard function, 
which must be replaced appropriately and can be taken directly from the corresponding 
threshold factorization theorem. One may also consider W production with W — )• iu, 
with an appropriate replacement of and Y with the charged lepton's rapidity. For a 
light Higgs with Q ~ mn, the isolated Drell-Yan factorization theorem applies to Higgs 
production through gluon fusion gg ^ H and Higgs-strahlung qq — )• VH, which are the 
dominant production channels at the LHC and Tevatron, respectively.^ For a generic process 
pp — )• XL, the sum over ij = {gg, uu, uu, dd, . . .} includes a gluon-gluon contribution, but 
still no cross terms between different parton types, and there will be two independent soft 
functions Sf^^^^ and Sf^^^^. [As shown in Ch. ^ only the qq soft function contributes 



to isolated Drell-Yan, so the labels were omitted in Eq. (2.16).] Indeed, the gluon-gluon 
contribution involving the gluon beam and soft functions, Bg and S^^j^j, gives the dominant 
contribution in the case of Higgs production. 

With the above physical picture, we can understand why the gluon beam function 
appeared in jp — )• J/ipX in the analysis of Ref. |96j in the limit where -Bj/^ — s- E.y. 
Taking px as the total momentum of final-state hadrons other than the J/ip, one has 
n-px ~ ^cm(l — Ej/^/E^), where n is the proton direction. For Ej/^ close to E^, energetic 
radiation in the final state is restricted to a jet close to the n direction. Just as for our 
B^f^, the measurement of -Ej/^ probes the radiation emitted by the colliding gluon in the 
initial state. Thus, the proton is broken apart prior to the hard collision, and the gluon 
beam function is required to describe the initial state. 



2.2 Generalized Observables 



The factorization theorem in Eq. (2.16) applies for ta ^ q and ^ q . This includes 
the situation where in the hadronic center-of-mass frame there is a numerically significant 
asymmetry uJa = XaEcm > = XhEcm- This means that the boost between the hadronic 
and partonic center-of-mass frames, given by the leptonic y = In y^uia/^b = In \/xa/xt,, is 
significantly different from zero. We explore the implications of this here. 

If there is no hierarchy, oja ^ i^b ^ \J^a^b = Q-, corresponding to F 0, we can define 



^In vector-boson fusion and associated production gg — >■ ttH, the situation is more complicated and one 
has to explicitly consider the process pp — >■ XjjH with two forward (top) jets. 
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a simple variable to constrain both hemispheres simultaneously. 



From Eq. (2.16), this gives 
1 da 



(To dq'^dYdB 



B= " ' ^ . (2.21) 



^Hij{q'^,^l) / dtadtbBi{ta,Xa,^i) Bj{tb,Xb,IJ.) 

xQSb(qB-^-^,^, (2.22) 



where the soft function is defined as 

SB{k+,i^) = I dktdk+ 5ihemi(A:^ , fc+ /i) dik+ -kt- k+) . (2.23) 

The advantage of using B is that the soft function now only depends on the single variable 
k^ = k'^ -\- k^ , much like the soft function for thrust in e~^e~ collisions. 

If we have a hierarchy cjf, < Q < uja, the final state has a substantial boost in the Ua 



direction, as shown in Fig. 2-3 In this case, the energetic radiation will generically be much 
closer to the beam axis in hemisphere a than in hemisphere b. To take this into account, it 
is natural to impose different cuts on B^ and i?^. Using the boost-invariant combinations 
LOaB^ /q"^ and UbB'^ /q^ to define the cut, we obtain 

= ^ < 6-2?^-* , = ^< , (2.24) 

q^ Ub q^ iOa 

so B'^ has a tighter constraint than B^ , as desired. If we simply replace B by B^ jujb + 
Bb l^ai the soft function analogous to Sb in Eq. (2.23) will depend on {ujaka +^bkb)/Q'^- 



However, we should also adjust the hemispheres themselves to take into account the 
significant boost of the partonic center-of-mass frame. We therefore define a generalized 
hemisphere o as y > 1" and hemisphere 6 as y < y, as shown in Fig. |2-3[ The corresponding 
total hemisphere momenta are denoted as B^^{Y) and the soft hemisphere momenta as 
k'^f^iY). The original definitions in Fig. 2-2 correspond to -^^^^(O) = i?^^ and A;^{,(0) = k'^^. 



The generalization of B is given by the boost-invariant combination 



u:aBt{Y)+LObBtiY) 
TB = ° ^ \ ^ ^ ^ . (2.25) 
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hemisphere b Y ,' hemisphere a 




Figure 2-3: Generalized definition of liemisplieres. Tlie total rapidity of the leptons is Y, 
Kb = na,b-ba,b, and A;+^(y) = Ua^b-kafiiX)- 



With the generalized definition of the hemispheres, i5^^(y) and uja^b transform under a 
boost by y in the Ua direction as 



Bt{Y)^B+'{Y + y) = e-yBtiY)^ 



B+{Y)^B+'{Y + y) = eyB+{Y), 



(2.26) 



Thus, boosting by y = —Y from the hadronic to the partonic center-of-mass frame gives 



UJ 



TB 



iBfiO) + u^B^'m _ BfiO) + i?+'(0) 



Q 



(2.27) 



In the partonic center-of-mass frame we have = = Q, so there is no hierarchy. 
Correspondingly, the generalized hemispheres in this frame are again perpendicular to the 



beam axis, so Eq. (2.27) has the same form as B. 



Note that for e~^e~ — t- jets, one can use the thrust axis to define two hemispheres with 
Ua^b analogous to our case. In the 2-jet limit, thrust is then given by 1 — T = {Qna-pXa + 
Q nb-pxi,) I'^Q'^ ■ Hence, we can think of tb as the analog of thrust for incoming jets. For 
this reason we will call tb the "beam thrust". 



In analogy to Eqs. (2.14) and (2.24), we define the cutoff on tb by 



TB < & 



(2.28) 
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For — )• or equivalently — )• oo the jets along the beam axes become pencil-Hke, while 
for generic y^* we allow energetic particles up to rapidities y < yg'* (with y measured in 
the partonic center-of-mass frame). 

The beam functions are boost-invariant along the beam axis, so the different hemi- 
sphere definitions do not affect them. The soft function is boost-invariant up to the 
hemisphere definition, which defines its arguments /c^^^. Hence, boosting by —Y we have 
5'ihemi[e^A:+,e~'^/i;j;'';y] = Si\,crm[K '^^'\ = Si\,cmi{K ^K)^ where the third argument de- 
notes the definition of the hemispheres. This implies that the soft function for tb is the 



same as in Eq. (2.23). The factorization theorem for tb following from Eq. (2.16) is 



1 da 

an dg^dydrR 



Q 

Integrating over < < exp(— 2?/^*) we obtain 



xQSb(Qtb-^^,/x). (2.29) 



do- , /•exp(-2s/g't) J 

ivT) = / dTB , ■ (2-30) 



dg2dy'^^ ' Jo dq^dYdTB 
We will use Eqs. (2.29) and (2.30) to show plots of our results in chapters [s] and [9l 



2.3 Relating Beam Functions and PDFs 



The beam functions can be related to the PDFs by performing an operator product expan- 
sion, because tafi ^ ^qcd- This yields the factorization formula 



/"^ d^ / X \ 



(2.31) 



where we sum over partons j = {g,u,u,d, . . .}, lij are perturbatively calculable Wilson 
coefficients, and fj is the standard PDF for parton j. The 0(AQ(^p/t) power corrections 



in Eq. (2.31) involve proton structure functions at subleading twist. Further mathemat- 



ical details on Eq. (2.31) are discussed in Ch. [3| whereas here we focus on the physical 
ramifications. 



The interpretation of Eq. (2.31) is illustrated in Fig. 2-4 At a hadronic scale fiA 
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fJLx changing x fl^ changing t /jLh 



Figure 2-4: Evolution of the initial state. Starting from the low scale //Aj the incoming 
proton is described by the x-dependent evolution of the PDFs, which redistributes the total 
momentum of the proton between its constituents. At the scale fiB, the proton is probed 
by measuring the radiation in the final state and breaks apart. This is the scale where the 
PDFs are evaluated and the x-dependent evolution stops. Above /j^b, the proton has ceased 
to exist, and the initial state behaves like an incoming jet, whose evolution is governed by 
the virtuality t of the off-shell spacelike parton that eventually enters the hard interaction 
at the scale /j^h- 



1 GeV, the initial conditions for the PDFs fj can be specified, and one has the standard 
DGLAP evolution up to the scale fiB, 



d 



m,^) = J2 /y (2-32) 



The anomalous dimensions Pjj' are the standard QCD splitting functions for quarks, an- 



tiquarks, and gluons (including the color factors and coupling constant). Equation (2.31) 
applies at the scale ^ = since this is the scale at which a measurement on the proton is 
performed by observing the soft and collinear radiation contributing to B^^^. At this scale, 
a parton j with momentum fraction ^ is taken out of the incoming proton according to 
the probability distribution /j (^, /u). As the parton continues to propagate and evolve with 
/i > /x^, it is modified by virtual radiation and by the emission of real radiation, which forms 
a jet. The evolution in this region no longer depends on ^, but instead on the virtuality t. 
This evolution occurs with fixed x and fixed parton type i, via the beam function RGE 

d f 
fi—B,{t,x,f^) = J dt'fB{t-t',^,)Bi{t',x,^l). (2.33) 

This result for initial-state jet evolution has the same structure as the evolution for final- 
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state jets. In fact, the anomalous dimension 7^ is identical to that for the quark jet function 
to all orders in perturbation theory. We further discuss the renormalization of the beam 
function and prove this correspondence with the jet function in Sec. |3.2[ 

The effect of initial-state real and virtual radiation is described by the perturbatively 
calculable Wilson coefficients x/^, fj,) at the scale fi = fiB- They encode several physical 
effects. The virtual loop corrections contribute to the la and modify the effective strength 
of the various partons. If the radiation is real, it has physical timelike momentum. Hence, 
it pushes the active parton in the jet off shell with spacelike virtuality — t < and reduces 
its light-cone momentum fraction from ^ to x. 

In addition, the real radiation can change the identity of the colliding parton, giving rise 



to the sum over j in Eq. (2.31). For example, an incoming quark can radiate an energetic 
gluon which enters the hard interaction, while the quark itself goes into the final state. This 
gives a contribution of the quark PDF to the gluon beam function through Igq. Similarly, an 
incoming gluon can pair-produce, with the quark participating in the hard interaction and 
the antiquark going into the final state. This gives a contribution of the gluon PDF to the 
quark beam function through Iqg. There are also of course real radiation contributions to 
the diagonal terms, Iqg and Igg, where the parton in the PDF and the parton participating 
in the hard interaction have the same identity. 

At lowest order in perturbation theory, the parton taken out of the proton directly enters 
the hard interaction without emitting radiation, 

lT''{t,^,f^)=6ij6{t)6(l-^). (2.34) 

Thus at tree level, the beam function reduces to the PDF 

Bl''^%t,x,^i) = 6{t)fiix,^,). (2.35) 

Beyond tree level, Xij{t,x/£^, n) can be determined perturbatively as discussed in more detail 
in Ch.[3j where we give precise field-theoretic definitions of the beam functions. We calculate 
the one-loop coefficients lij for the quark beam function in Ch. [4] and for the gluon beam 
function in Ch. [5l 



Interestingly, in the threshold factorization theorem Eq. (2.9), cross terms between quark 



and gluon PDFs are power suppressed, so the gluon PDF does not contribute at leading 
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order. In the inclusive case Eq. (2.6), such cross terms are leading order in the power 
counting. For isolated Drell-Yan, there are no cross terms between quark and gluon beam 
functions, but there are leading-order cross terms between different PDFs, which appear 



via the contributions of different PDFs to a given beam function in Eq. (2.31). Thus, the 
isolated case is again in-between the inclusive and threshold cases. 



2.4 Comparison with Initial-State Parton Shower 

The physical situation associated with the beam evolution has an interesting correspondence 
with that of initial-state parton showers. As pictured in the region between and in 



Fig. 2-4, the parton in the beam function evolves forward in time while emitting a shower of 



radiation into the final state governed by the anomalous dimension 7^(t— i', fi) in Eq. (2.33). 
This equation has no parton mixing. Each emission by the radiating parton increases the 
magnitude of its spacelike virtuality —t < 0, pushing it further off-shell in a spacelike 
direction. At the time the parton is annihilated in the hard collision, it has evolved to some 
t with |t| <C g^, so the large momentum transfer q'^ guarantees that no partons in the final 
state are spacelike. This description agrees quite well with the physical picture associated 
with the evolution of the primary parton in an initial-state parton shower, as summarized 
in Ref. fT59]. 

Differences in the description arise when one considers the initial-state parton shower in 
more detail (for simplicity we focus on the so-called longitudinal evolution). The shower is 



based on the evolution equation for the PDFs in Eq. (2.32 ). An evolution forward in time is 
not practical because of the lack of prior knowledge of the scale of the hard interaction, so 
the shower uses backward evolution starting at a given partonic hard scale Q |158| . Knowing 
the identity of the final parton i, the shower evolves based on the probability dVi/dt that 
parton i is unresolved into parton j via the splitting j — )• ij' at an earlier (lower) scale t. 
The evolution equation is |159j 



dt 



D ^^ fjix/z,t) 

fi{x,t) 



\Vi{x,t^^^,t) , (2.36) 



where 'Pi(x, tmax, t) is the shower Sudakov exponential, which is interpreted as the probabil- 
ity for no emissions to occur between the initial value tmax and t. The evolution variable i, 
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which determines the scale of the sphtting, is usually chosen as the virtuality or transverse 
momentum of the parton. 

The mixing of partons in the PDF evolution influences the shower. In particular, the 
evolution kernel depends on the PDF fj{x/z,t), which determines the number density of 
partons of type j at the scale t, and inversely on the PDF fi{x,t). Thus, unlike in the 



beam evolution in Eq. (2.33), the shower evolution in Eq. (2.36) still knows the identity of 
the initial-state hadron. Double logarithms in the initial-state parton shower are generated 
m q ^ qg and g ^ gg splittings because of the soft-gluon singularity ~ 1/(1 — 2;) in the 
splitting functions. This singularity is regulated [159] by the upper cutoff -Zmax = x/{x + Xe), 
where x^ provides a lower cutoff on the gluon energy in the rest frame of the hard scattering, 
Eg > Xe7i?cm/2 — 2GeV (where 7 is the boost factor of the hard scattering). Hence, one 
logarithm, Inx^, is generated by the z integration, and one logarithm, Int, by the collinear 
1/t singularity. In contrast, the beam function contains double logarithms In^t similar to a 
final-state parton shower, where the z integration yields a kernel ~ (lnt)/t that produces a 
double logarithm In^ t via the t evolution. 

The above comparison is very rough. For example, the influence of soft radiation on 
both the shower and on the isolated factorization theorem was not compared and is likely 
to be important. Furthermore, the goal of the shower is to provide a universal method for 
populating fully exclusive final states, while the beam function applies for a more inclu- 
sive situation with a particular measurement. Note that just the presence of mixing in the 
initial-state parton shower and absence of mixing in the beam-function evolution does not 
imply an inconsistency. For example, it is well known that the final-state parton shower 
reproduces the correct double logarithms for e~^e~ event shapes |64] . even though there is 
no parton mixing in the evolution of the corresponding hard, jet, and soft functions. In the 
future it would be interesting to test in detail the correspondence between the double loga- 
rithms generated by the initial-state parton shower and those predicted by our factorization 
theorem for the isolated Drell-Yan process. 



2.5 Relation to Fixed-Order Calculation 



The factorization theorem for the cross section in Eq. (2.16) and the factorization for the 



beam function in Eq. (2.31 ) together allow us to describe in more detail how various Feyn- 
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Figure 2-5: Factorization for isolated Drell-Yan in pictures. The left-hand side of each 
equality are graphs in QCD, while the right-hand side shows the sum of the corresponding 
SCET diagrams. Dashed lines are collinear quarks, and springs with a line through them 
are collinear gluons. The double lines denote soft Wilson lines, and the gluons attached to 
them are soft. 



man diagrams that would appear in a fixed-order calculation contribute to the cross section 



in our kinematic region. Various examples are shown in Fig. 2-5 
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we have the tree-level qq annihilation producing a 7 or Z, which involves 
the tree-level 0{a^) hard function, beam functions, and soft function, denoted by a super- 
script (0) in the figure. In Fig. 2-5(b)[ initial-state gluons couple to a quark loop (e.g. a top 



quark), which subsequently annihilates into a 7, Z, or Higgs. The quarks in this loop are 
far off shell, so they can be integrated out and appear as one-loop corrections, Hg]i\ to the 
hard coefficient in the factorization theorem. Other possibilities for this graph are power 
suppressed. 



The situation for the vertex correction in Fig. 2-5 (c) is more involved. If the gluon 



in the loop is hard, all particles in the loop are far off shell and can be integrated out, 
giving the one-loop hard function Hj^^^ shown as the first term on the right-hand side. In 
the second term, the gluon is collinear to the incoming quark beam and gives a virtual 
one-loop contribution to the quark beam function, Bq^\ The third term is the analog of 
the second, but now with the gluon collinear to the incoming antiquark. Finally in the 
fourth term, the gluon is soft, communicating between the incoming collinear beams. Here, 
the eikonal approximation holds for describing the quark propagators. The generalization 
of this to all orders in as leads to the fact that the soft function is a matrix element of 
Wilson lines. Although a single loop graph contributes in several different places in the 
factorization theorem, all of these contributions have a precise separation in SCET. We will 
use this separation in Ch. [7] to prove the isolated Drell-Yan factorization theorem. 



An interesting contribution occurs in Fig. 2-5(d), where a gluon is radiated into the 
final state. Because of the kinematic restrictions in isolated Drell-Yan, this gluon can only 
be collinear to the incoming quark, collinear to the incoming antiquark, or soft, and these 
three possibilities are represented by the diagrams on the right-hand side of the equality. 
In the first case, we have a real-emission correction to the quark beam function, Bq^\ In 
the second case, the intermediate quark is far off shell and can be integrated out, and the 
gluon collinear to the antiquark arises from a collinear Wilson line contribution in B^^\ 
The third case gives a real-emission correction to the soft function, S^^ . The full-theory 



graph in Fig. 2-5 (d) has a t-channel singularity. An important fact about the isolated Drell- 
Yan factorization theorem is that it fully captures the dominant parts of this singularity, 
and allows a simple framework for a resummation of higher order as corrections enhanced 
by large double logarithms due to this singularity. For threshold Drell-Yan, the kinematic 
restrictions are stronger and only allow the third graph with soft initial-state radiation. In 
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inclusive Drell-Yan, the gluon is treated as hard, and the graph in Fig. 2-5(d) only corrects 



H^^^, without providing a framework for summing the large double logarithms that appear 
when we make a global measurement of the radiation in each hemisphere defined by the 
beams. 



The situation is a bit simpler for Figs. 2-5 (e) and 2-5 (f) In Fig. 2-5 (e), the incoming 
collinear gluon from the PDF pair-produces a quark and antiquark both collinear to this 
beam direction, and the quark enters the hard interaction. Therefore, this is a one-loop 
correction to the quark beam function, Bq^\ proportional to the gluon PDF fg. The 
beam functions again allow us to resum the possibly large logarithms due to this t-channel 
singularity. Other possibilities for the final-state antiquark in Fig. 2-5(e)| lead to power- 



suppressed contributions. Similarly, the s-channel graph in Fig. 2-5(f), which has the same 



initial and final states as Fig. 2-5 (e) , has no leading-power contribution and only contributes 



to Eq. (2.16) in the power-suppressed terms. The same is also true for Drell-Yan in the 
threshold region. Only inclusive Drell-Yan receives a leading-order hard contribution from 
the s-channel graph, which is then treated as of the same size as the t-channel graphs. 



2.6 Renormalization Group Evolution 

In this section, we discuss and compare the structure of large logarithms in the cross sections 
for inclusive, threshold, and isolated Drell-Yan. These large logarithms may be summed 
using the renormalization group evolution of the individual functions appearing in the fac- 
torization theorems. In fact, the structure of large logarithms in the differential cross 
section allows us to infer the necessity of the beam functions in the isolated factorization 
theorem. This procedure provides a method of determining whether beam functions enter 
for other observables or processes than those studied here. The consistency of the RGE was 
used to provide a similar consistency check in Ref. when deriving a new factorization 
theorem for the invariant-mass distribution of jets initiated by a massive quark in e~^e~ col- 
lisions. In that case, the RGE consistency provided important constraints on the structure 
of the factorization theorem at scales below the heavy-quark mass. 

In inclusive Drell-Yan, the hard functions Hfj'^^ are sensitive to the scale hh — Q oi 
the hard interaction, and the proton mass defines a low scale /iA — 1 GeV > Aqcd (which 
is still large enough so perturbation theory can be applied for the PDF evolution). The 
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Figure 2-6: RGE running for different Drell-Yan scenarios. Case (a) corresponds to the 
inclusive case. Case (b) corresponds to the threshold case, where the kinematics forces all 
hadrons in the final state to be soft. Case (c) corresponds to the isolated case. Here, the 
PDFs freeze out at the intermediate beam scale fi b , above which they are replaced by beam 
functions. 



measurement of and Y in this case does not introduce additional scales, and thus does 
not influence the structure of the logarithms. Thus, we have the hierarchy <^ fj,H, 
and the large logarithms are L = \d.{ii\/ ^u). Here, only single-logarithmic series, [agL)^, 
are generated at higher orders in perturbation theory. The logarithms are factorized as 



ln{ii/iiH) + lii(/iA//^) in the factorization theorem in Eq. (2.6) and may then be resummed. 
The general form of the running is pictured in Fig. |2-6[ a). The logarithms ln(/iA//i) are 
summed by evolving the PDFs fi{^a^^J') and fj{(,b,lJ') from ^ua up to the common scale 
/X. The inclusive hard function, W'^^^{xa/Ca,xi,/Cb,Q^, fJ-), is evolved from fin down to /i, 
summing the logarithms ln(/i//i//). The choice of /i is arbitrary. Taking fi ~ fin corresponds 
to only running the PDFs up, while for fi fi/<^ only _f/'™^' runs down. The equivalence of 
these two choices implies that must be convoluted with the two PDFs and exhibit a 

factorized structure for logarithms in the a and b variables. 

For threshold Drell-Yan, the kinematic restrictions only allow soft radiation in the final 
state. This induces additional large logarithms ln(l — r). These can be written in terms 
of a ratio of scales ln(fj,s/ f^n), where the soft scale fis — (3(1 — t) is another important 
scale in the analysis. The logarithms L = \n{^s / I^h) appear as double- logarithmic series 



{agL ) in the cross section. In the threshold factorization theorem in Eq. (2.9), these 



double logarithms can be summed by evolving the PDFs and the threshold soft and hard 
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functions, S'thr and H, to a common scale /x, as shown in Fig. 2-6 'b). Since ^a,b — ^ 1; 
the logarithms ln(l — ^q) and ln(l — S^h) are also large. The RGE for the PDFs must be 
expanded, and the result sums a double- logarithmic series of ln^(l — ^) terms. The threshold 
soft function sums double logarithms ln'^{fi/^s) between fis and fi, while the threshold hard 
function sums double logarithms \v?{^/ ^h) between and /i. The RG equations are 

= (2.37) 

The consistency of the RGE at the scale /i shown in Fig. |2-6[ b) implies that the double 
logarithms in /j, fj, and Sthr combine in such a way that the RGE of the convolution 
fifj ^ Sthr is identical to that of H, and hence only depends on q^. 

For isolated Drell-Yan, the kinematic restrictions allow both soft and coUinear initial- 
state radiation, and induce an invariant-mass scale for each beam function, fi^ ~ XaEcmB^ 
and ~ XbEcmB^ , and a soft scale fis ~ B^b- simplicity, we use a common scale 
^jLB for both beam functions in our discussion here. (Since the evolution of the two beam 
functions is independent, one can just as easily implement two independent beam scales.) As 



we saw in Sec. 2.1 , at partonic center-of-mass energies of a hundred GeV to a few TeV there 
is a large hierarchy between the different scales, ^ /^s ^ /^s ^ ^/f, and correspondingly 
large double and single logarithms of the ratios of these scales. The RGE running for this 
case is shown in Fig. |2-6| ^c). Here, the PDFs are not restricted to their endpoints, so their 
evolution is given by Eq. ( |2.32 ), which involves the unexpanded and nondiagonal Pij{S,/^') 



and sums single logarithms, (a^L) . For each fj this evolution joins at ^ = fiB with the 



Wilson coefficients lij in the beam function factorization Bi = lij (8) fj of Eq. (2.31). The 
lij cancel the ^-dependent evolution of fj, and turn it into the t-dependent evolution of 
Bi, which sums a double-logarithmic series. The objects meeting at the common scale /u in 
Fig. |2-6|^c) are the hard function (same as threshold case) and the beam and soft functions, 



d f 
tJ'-^Bi{t,x,fi) = J dt'7B(i - t',/^) Bi{t',x,n), 

d/i 



f^-^Sii,emi{kt,k^,f^) = I dfc^,dA;[,5ihemi(^a.^b.^)75ihcmi(^^ - K^k^ - k'b,fJ-) ■ (2.38) 
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The consistency of the RGE at fi now imphes that the double-logarithmic running in the 
different variables for Bi, Bj, and Sihemi cancels such that the convolution BiBj (gi 5ihcmi 
has an RGE identical to H, which only depends on q^. (A detailed discussion of this 
consistency can be found in Ref. [94] for the analogous case of two jet functions and the 
final-state hemisphere soft function, J J S'hemi-) It is important that this cancellation 
would not be possible if we tried to replace Bi by fi in the isolated factorization theorem. 
Given the type of double logarithms in the cross section, the single logarithms summed by 
the PDFs at generic x cannot combine with the double logarithms in S'ihemi to give a result 
in agreement with the double logarithms in H. Thus, the structure of double logarithms 
necessitates the presence of beam functions in the isolated factorization theorem. 

By the same argument we can conclude that for all processes involving a threshold-type 
hard function H with double logarithms, and with Xa^b away from one, the description 
of the initial-state radiation will require beam functions B^. This includes all situations 
where H is the square of Wilson coefficients of SCET operators, H = X^fcl^'fcP (for example 
when the energetic partons in the hard collision all have distinct collinear directions). In 
particular, the theoretical description of any threshold process with x — )• 1 can be extended 
to a factorization theorem for the respective isolated case with x away from one. This is 
achieved by adding variables -B^^, replacing the PDFs by beam functions, and replacing 
the threshold soft function by an appropriate soft function for the isolated case. 

Thus, beam functions are quite prevalent for cross sections that one may wish to study 
at the LHC. In situations where the hadronic final state is constrained with variables that 
are more complicated than B^^, one generically expects to find different beam functions 
and different soft functions encoding these constraints. This extension is analogous to how 
the choice of jet algorithm modifies the definition of the jet and soft functions for central 
jets produced by the hard collision [32]. Even with this generalization, the beam and soft 
functions will both sum double-logarithmic series, and we expect that the factorization 
relating the beam function to the PDFs will carry through, just with different Wilson 
coefficients la. 
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Chapter 3 

Beam Functions 



We discuss several formal aspects of the quark and gluon beam functions in this chapter, 

as 



which were first reported in Ref. |168j . The beam functions are defined in Sec. 3.1 



matrix elements of operators in SCET. We include a brief review of the necessary SCET 



ingredients. In Sec. 3.2 we work out the renormalization of the beam functions and show 
they have the same RGE as jet functions to all-orders in perturbation theory. We relate 



beam functions to PDFs in Sec. 3.3 by performing an OPE, and calculate the tree-level 



matching coefficients in Sec. 3.4 Finally, in Sec. 3.5 we discuss the analytic structure of the 
beam functions and derive a relationship with matrix elements of time-ordered products of 
fields. This allows us to use the usual Feynman rules when we calculate the beam functions 
in the next chapters. 



3.1 Definition 

In this section we discuss the definition of the quark and gluon beam functions in terms of 
matrix elements of operators in SCET, and compare them to the corresponding definition 
of the PDF. The operator language will be convenient to elucidate the renormalization 
structure and relation to jet functions in the following section. 

We first discuss some SCET ingredients that are relevant later on. We introduce light- 
cone vectors n'^ and n'^ with = = and n • n = 2 that are used to decompose 
four- vectors into light-cone coordinates = {p'^ ,P~ tP±), where = n ■ p, p~ = n ■ p and 

contains the components perpendicular to and n^. 

In SCET, the momentum p^ of energetic collinear particles moving close to the n direc- 
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tion is separated into large and small parts 



=P^ +p'^. =n-p—+p^'^^+pi; . (3.1) 

The large part = (0,p~,pn±) has components p~ = n ■ p and pn± ~ Ap~, and the 
small residual piece pr = {pt -^PV iPta) ~ P~ i ■> with A ^ 1. The corresponding 
n-collinear quark and gluon fields are multipole expanded (with expansion parameter A). 
This means particles with different large components are described by separate quantum 
fields, in,p{v) and An.p{y), which are distinguished by explicit momentum labels on the 
fields (in addition to the n label specifying the collinear direction). We use y to denote the 
position of the fields in the operators to reserve x for the parton momentum fractions. 

Interactions between collinear fields cannot change the direction n but change the mo- 
mentum labels to satisfy label momentum conservation. Since the momentum labels are 
changed by interactions, it is convenient to use the short-hand notations 

Uv) = E ^r.M > A'niy) = E <dy) ■ (3.2) 

The sum over label momenta explicitly excludes the case = to avoid double-counting 
the soft degrees of freedom (described by separate soft quark and gluon fields). In practice 
when calculating matrix elements, this is implemented using zero-bin subtractions |140j or 
alternatively by dividing out matrix elements of Wilson lines |120| I134j . The dependence on 
the label momentum is obtained using label momentum operators Vn or V^j_ which return 
the sum of the minus or perpendicular label components of all n-collinear fields on which 
they act. 

The decomposition into label and residual momenta is not unique. Although the ex- 
plicit dependence on the vectors and breaks Lorentz invariance, the theory must 
still be invariant under changes to and which preserve the power counting of the 
different momentum components and the defining relations = v? = 0, n ■ n = 2. This 
reparametrization invariance (RPI) \Q9\ I139j can be divided into three types. RPI-I and 
RPI-II transformations correspond to rotations of n and n. We will mainly use RPI-III 
under which and transform as 



(3.3) 
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which impHes that the vector components transform as p"*" — )• e"p^ and p 



. In 



this way, the vector p^ stays invariant and Lorentz symmetry is restored within a cone about 



the direction of n^. Since Eq. (3.3) only acts in the n-colhnear sector, it is not equivalent 



to a spacetime boost of the whole physical system. 
We now define the following bare operators 



j^bare 



q {y 
T'%y-,u^) 



e-^^^y-/'Xn{y-'^)^[6{u:-Vn)Xn{0)] 



-ip+y-/2 



tr<^Xn{y-'^)[6{u;-V. 



-ip+y /2^c 



'n±^,[y 



_ n 



n)Xn(O)]}, 



[<5(.. - P„)e^l(0)] . 



(3.4) 



Their renormalization will be discussed in the next section. The corresponding renormalized 



operators are denoted as Oi{y ,u),n) and are defined in Eq. (3.21) below. Here, p~^ is the 



momentum operator of the residual plus momentum and acts on everything to its right. 
The overall phase is included such that the Fourier-conjugate variable to corresponds 



to the plus momentum of the initial-state radiation, see Eq. (3.9) below. The operator 



6{uj — Vn) only acts inside the square brackets and forces the total sum of the minus labels 
of all fields in Xn(0) and Bn±{^) to be equal to lo. The color indices of the quark fields are 
suppressed and summed over, c is an adjoint color index that is summed over, and the trace 
in Oq is over spin. The operators are RPI-III invariant, because the transformation of the 
5{ijJ — Vn) is compensated by that of the ^ in Oq^q and the overall uj in Og. 



The fields 



Xn{y) = wl{y)Uy) 



(3.5) 



with ii?^_|_ = + are composite SCET fields of n-collinear quarks and gluons. In 



Eq. (3.4) they are at the positions y^ = y and y'^ = 0. The Wilson lines 



Wn{y) 



perms 



n-An{y) 



(3.6) 



are required to make Xn{y) and B'^_^[y) gauge invariant with respect to collinear gauge 
transformations [3T[ |37j. They are Wilson lines in label momentum space consisting of 
n-An{y) collinear gluon fields. They sum up arbitrary emissions of n-collinear gluons from 
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an n-collinear quark or gluon, which are 0(1) in the SCET power counting. Since Wniy) is 
locahzed with respect to the residual position y, Xniu) and 13^j_iy) are local operators for 



soft interactions. The fields in Eqs. (3.4) and (3.5) are those after the field redefinition [35 



decoupling soft gluons from collinear particles. Thus at leading order in the power counting 
these collinear fields do not interact with soft gluons through their Lagrangian and no 



longer transform under soft gauge transformations. Hence, the operators in Eq. (3.4) are 
gauge invariant under both soft and collinear gauge transformations. The soft interactions 
with collinear particles are factorized into a soft function, which is a matrix element of soft 



Wilson lines (see Sec. 7.3.3). 



Note that our fields in Eq. (3.4) have continuous labels and hence are not the standard 
SCET fields with discrete labels. They only depend on the minus coordinate, y~ , corre- 
sponding to the residual plus momentum, p+, and not a full four- vector y^. As discussed in 
detail in the derivation of the factorization theorem in Ch. [7| it is convenient to absorb the 
residual minus and perpendicular components into the label momenta which then become 
continuous variables. For example, for the minus momentum (suppressing the perpendicular 
dependence) 

E e-'^'^^/'xn,p-(y-,y+) = E/ e-'^^"+"^")^^/\n,p-(y-) 

p- p- 

'dp-e-'p-y^/^xn,p-iy-)- (3.7) 

In this case, Wn{y^n/2) can also be written in position space where all gluon fields sit 
at the same residual minus coordinate, y~, and are path ordered in the plus coordinate 
(corresponding to the label minus momentum) from y~^ to infinity. 

Next, we introduce the Fourier-transformed operators 

0^-(u;6+,u;) = ^ e'^''^'/^ a^-(y-,a;) . (3.8) 

zvr J z\uj\ 

For example, for the quark operator 

O^^^^^icob+,u;) = ^ |^e'(^+-^^)^-/2(ei^-^^'/2x„(0)e-'^^+^-/2)|[5(.;-P„)Xn(0)] 

= XniO) Siu;b+ - ojp+) I [6{u; - Vn)Xn{0)] . (3.9) 
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In the first step we used residual momentum conservation to shift the position of the field. 



Here we see that the overall phase in Eq. (3.4) allows us to write the dependence in terms 



of 6{ujb'^ — ojp'^), which means that measures the plus momentum of any intermediate 
state that is inserted between the fields. For the corresponding renormalized operators 



Oi{ujb'^ fi) see Eq. (3.22) below 



We divide by |a;| in Eq. (3.8) to make the integration measure of the Fourier transform 



RPI-III invariant. (Taking the absolute value ensures that the definition of the Fourier 
transform does not depend on the sign of uj.) As a result, the Fourier-transformed operators 
are still RPI-III invariant and only depend on 6^ through the RPI-III invariant combination 
t = ub~^ . The beam functions are defined as the proton matrix elements of the renormalized 
operators Oi{t,uj, fi), 



(3.10) 



The matrix elements are always averaged over proton spins, which we suppress in our 



notation. Note that part of the definition in Eq. (3.10) is the choice of the direction n such 



that the proton states have no perpendicular momentum, = P~n^/2, which is why we 
write \pn{P^))- By RPI-III invariance, the beam functions can then only depend on the 
RPI-III invariant variables t = ojb'^ and x = uj/P^ . The restriction 9{u}) on the right-hand 



side of Eq. (3.10) is included to enforce that the XniO), Xn(0), or Bn±iO) fields annihilate a 



quark, antiquark, or gluon out of the proton, as we discuss further at the start of Sec. 3.5 



The definition of the beam functions can be compared with that of the quark and gluon 
PDFs. In SCET, the PDFs are defined [30] in terms of the RPI-III invariant operators 

Q^-^(a;') = eiu;') Xn(0)| [5(^' - P„)xn(0)] , 
Qf'^u;') = e{J) tr{|xn(0) [5{J - P.)xn(O)] } , 



-ijj 



(3.11) 



as the proton matrix elements of the corresponding renormalized operators Qi{uj' , //) defined 



in Eq. (3.14) below, 



Moo'/p-,fl) = {pn{P-)\Q^{i0',^l)\pn{P-)). 



(3.12) 
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By RPI-III invariance, the PDFs can only depend on the momentum fraction ^ = w'/P . 
Beyond tree level ^ or uj' are not the same as a; or a; , which is why we denote them differently. 



Without the additional ^(w') in the operators in Eq. (3.11) the quark and anti-quark PDFs 
would combine into one function, with the quark PDF corresponding to a; > and the 
antiquark PDF to a; < 0. We explicitly separate these pieces to keep analogous definitions 
for the PDFs and beam functions. 



It is important to note that the coUinear fields in Eq. (3.11) do not require zero-bin 



subtractions, because as is well-known, the soft region does not contribute to the PDFs. 
If one makes the field redefinitions ^„ — t- Y^n and — >• YAnY^ to decouple soft gluons. 



then the soft Wilson lines Y cancel in Eq. (3.11). Equivalently, if the fields in Eq. (3.11) 



include zero-bin subtractions then the subtractions will cancel in the sum of all diagrams, 
just like the soft gluons. (This is easy to see by formulating the zero-bin subtraction as a 
field redefinition [T34J analogous to the soft one but with Wilson lines in a different light- 



cone direction.) In contrast, the collinear fields in the beam function operator in Eq. (3.4) 



must include zero-bin subtractions. We will see this explicitly at one loop in our PDF and 
beam function calculations in Ch. 01 

The SCET definitions of the PDFs are equivalent to the standard definition in terms of 
full QCD quark fields ■0 in position space. For example, the quark PDF in QCD is defined 
as [71] 



dy+ 
in 



p„(p-)|[v;(y+^)iy(y+^,0)|V'(0)] \pn{P 



(3.13) 



The square brackets denote the renormalized operator. Here, the fields are separated along 
the n direction and the lightlike Wilson line W{y^n/2,0) is required to render the product 
of the fields gauge invariant. The relation to the SCET definition is that the SCET fields 



in Eq. (3.11) (without zero-bin subtractions) involve a Fourier transform of "0 in U'^ to give 



the conjugate variable lo' . The corresponding Wilson lines in Eq. (3.12) are precisely the 
Wn contained in the definitions of Xn and Hence, the QCD and SCET definitions of 

the PDF are equivalent (provided of course that one uses the same renormalization scheme, 
which we do). 



Comparing to Eq. (3.4), the difference between the beam functions and PDFs is that for 



the beam functions the fields are additionally separated along the n light-cone, with a large 
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separation ^ y+ corresponding to the small momentum <^ uj. Thus, formulating 
equivalent definitions of the beam functions directly in QCD would be more challenging, as 
it would require QCD fields that are simultaneously separated in the n and n directions. 
For this case, it is not clear a priori how to obtain an unambiguous gauge-invariant defi- 
nition, because Wilson lines connecting the fields along different paths are not equivalent. 
This ambiguity is resolved in SCET, where the multipole expansion distinguishes the dif- 
ferent scales and divides the possible gauge transformations into global, collinear, and soft 
transformations, allowing one to treat the separations along the two orthogonal light-cones 
independently. The large y~ separation corresponds to soft Wilson lines and soft gauge 
transformations that are independent from collinear gauge transformations corresponding 



to the small dependence. As already mentioned, the operators in Eq. (3.4) are separately 
gauge invariant under both types of gauge transformations. 



3.2 Renormalization and RGE 

The beam functions and PDFs are defined as the matrix elements of renormalized operators. 
The renormalization of the operators immediately yields that of the functions defined by 
their matrix elements. In this section we derive the RG equations and show that the anoma- 
lous dimensions of the beam and jet functions are the same to all orders in perturbation 
theory. 

We start by considering the known renormalization of the PDF, but in the SCET op- 
erator language. The renormalized PDF operators are given in terms of the bare operators 



in Eq. ( |3.11[ ) as 

'do;' 



j 

In general, operators with different i and w can (and will) mix into each other, so the 
renormalization constant z(-{uj/ijj' , ^) is a matrix in i,j and (jJ,Ljj' . RPI-III invariance then 
restricts the integration measure to be dw'/w' and zfj{uj/oj', /u) to only depend on the ratio 



u/uj' . Hence, the form of Eq. (3.14) is completely specified by the SCET symmetries. The 
independence of the bare operators Q^^'^'^(a;) yields an RGE for the renormalized operators 
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in MS 



j 



(3.15) 



where the inverse (Z-^) .^^(z, ^) is defined as 

•dz' 



z). 



(3.16) 



Taking the proton matrix element of Eq. (3.15) yields the RGE for the PDFs 



d/x 



m). 



(3.17) 



The solution of this RGE can be written in terms of an evolution function which acts 
on the initial PDF /;(■?', /^o) and takes it to fiiCfJ-), 



/.(e,M)= /yC^i(|,;",/^o)/,(e',M). 



(3.18) 



Prom Eq. (3.17) we can identify the anomalous dimensions 'jfj{z) in terms of the QCD 



splitting functions. For example, in dimensional regularization in the MS scheme, the one- 
loop anomalous dimensions for the quark PDF are the standard ones 



TT 



TT 



(3.19) 



with the q ^ qg and g ^ qq splitting functions 



P^g{z) = Co{l-z){l + z') + -5{l-z) 



2n 



9(1 - z) 



1 + z 
1 - z 



P^^(z) = eii-z)[{i-zf + z'] 



(3.20) 



The plus distribution Co{x) = [9{x)/x]+ is defined as usual, see Eq. (B.2). For later 



convenience we do not include the overall color factors in the definitions in Eq. (3.20) 



We now go through an analogous discussion for the beam functions. The renormalized 
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operators Oi{y ,uj,fi) are given in terms of the bare operators in Eq. (3.4) by 



O^- iy-,uj) = Z},(y-,fAo,{y-,u,f,) 



\2uj' 



(3.21) 



where Z^(y^/2aj, /x) is the position-space renormahzation constant. In App. [c| we give 
an expHcit proof that the beam function renormahzation is multiphcative in this way to 
all orders in perturbation theory. The underlying reason is that the renormahzation of the 
theory should preserve locality, so renormalizing the nonlocal beam function operator should 
not affect the y~ separation between the fields. For example, mixing between operators with 
different y~ would destroy locality at distance scales within the validity range of the effective 
theory. RPI-III invariance then implies that can only depend on the ratio /2uj (the 
factor of 1/2 is for convenience). In principle, one might think there could also be mixing 



between operators with different z or w in Eq. (3.21) [as was the case for the PDFs in 



Eq. (3.14)]. Our derivation in App. [C] shows that this is not the case. 



Taking the Fourier transform of Eq. (3.21) according to Eq. (3.8), we find 



©^■-^(t, a;) = / dt' ZUt - t', fi) 0,{t', oj, f,) , 



1 



2tt / 2 a; 



(3.22) 



Since the bare operator is fi independent, taking the derivative with respect to //, we find 
the RGE for the renormalized operator 



/i— Oi(i, w, /x) = / dt'^sit - t', fi) Oi{t', uj, /i) 



7B(i,M) 



dt'iZh)-Ht-t',f,)i,^ZW,t,) 



(3.23) 



where the inverse of Z'^Q{t,fi) is defined as usual. 



dt'{Z'B)-\t-t\i,)ZUt\li) = 5{t) 



(3.24) 



Taking the proton matrix element of Eq. (3.23) we obtain the corresponding RGE for 
the beam function, 



H—Bi{t,x,n) = I dt' -i'B{t-t\n)Bi{t' ,x,ijl) . 



(3.25) 
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As discussed in App. [Cj to all orders in perturbation theory the anomalous dimension has 
the form 

YB{t,fi) = -2r^„,p(a,)-^£o(^) +jh{o^s)S{t), (3.26) 



where Co{x) = [9{x)/x]^ is defined in Eq. ( |B.2| ), r*^gp(as) is the cusp anomalous dimension 
for quarks/antiquarks (i = q) or gluons {i = g), and 7^(0^) denotes the non-cusp part. Since 
there is no mixing between operators Oi{t,oj,^) with different i or uj^ the beam function 
RGE only changes the virtuality t but not the momentum fraction x and does not mix 
quark and gluon beam functions. By rescaling the plus distribution, 



(3.27) 



we can see that 7^(t, ^u) has logarithmic /^-dependence, which means that the RGE sums 
Sudakov double logarithms. 



The solution of the RGE in Eq. (3.25) with the form of the anomalous dimension in 



Eq. (3.26) is known [26l [Ml EM] • It takes the form 



Bi{t,x,fj,) = J dt' Bi{t - t',x,no) UB{t',fio,n) 



(3.28) 



where the evolution kernel can be written as [135] 



W+vh) 



t 



Mo ^Mo 



(3.29) 



The distribution C^{x) is defined in Eq. (B.2), and the RGE functions K^^ = K^(/iO)M) 



and = 77^(^x0) m) ^"^^ given in Eq. (E.8). 



The SCET quark, antiquark, and gluon jet functions are given by [351 195] 



(2vr)2 [dy 



(3.30) 



(2vr)2 rdy- 



Nr J 2 a; 



Xn (y' ^) I [5(^ + Vn)S\uj± + Vn±)Xn{0)] 



, 
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(2^)2 fdy- 



- 1 ./ 2\lo 



where the notation [. . again denotes the renormahzed operators. Here, we used the 
same conventions as for the beam functions where the large label momenta oj and a;_L are 
continuous, so the only position dependence of the fields is in the minus component. RPI 
invariance requires that the jet function only depends on the total invariant mass of the jet, 
= bjp'^ + w^. When the jet function appears in a factorization theorem, the direction of 
the jet is either measured (e.g. by measuring the thrust axis in e^e~ — >• 2 jets) or fixed by 
kinematics (e.g. in i? — )• Xg'^ the jet direction is fixed by the direction of the photon) and 
n is chosen along the jet direction, so one typically has (jJ±_ = 0. Taking the vacuum matrix 
element of 0^^^'^{t,u:), we get 

^'^''^'^-{0\oY'%-t,-uj)\{)) 



= J (fuj±J^'''%t - CSl) = J^^'%t) = J^^'^t) . (3.31) 

In the last step we used that the quark and antiquark jet functions are the same. The 
analogous relation holds for the antiquark operator, Oq{t, lv, /x). The uj± integral is bounded 
and does not lead to new UV divergences, because the jet function only has support for 
nonnegative argument, < < t, and t is fixed. Similarly, for the gluon operator we have 



(2^) 



'l^{0\O^^-^'{-t, -a;)|0> = I d^^x Jb-<=(t - (^i, ^) = jj-^(i) . (3.32) 



The renormalization of J^^^^{t) does not depend on the choice of u± in Eq. (3.30). Since 
Jr'%t) is simply an average over different choices for uj± it has the same renormalization. 
Hence Ji{t,^) and Ji{t,fj,) have the same anomalous dimension, 

fi-^Jiit, fi) = j d^u^ds-i'j{t - s,n) Ji{s,ij) = Jdt'j){t-t',n) J d'^u}± Ji{t' - ul, fi) 
''dt'J){t-t',^l)Ji{t',^i). (3.33) 
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On the other, Eq. (3.31) says that Jr'^it) is a matrix element of the same operator 



as the beam function, whose general renormalization was given in Eq. (3.22). We thus 



conclude that the beam and jet function anomalous dimensions are identical to all orders 
in perturbation theory, 

fj{t,fi)=jUt,f^)- (3.34) 



For the cusp part this result already follows from our explicit one-loop calculation, since 



cusp 



is universal and its coefficients are the same at one loop. Our one-loop result provides 



a cross check for the identity of the one-loop non-cusp part of the anomalous dimension, 
which agree. Furthermore, jj{as) and hence 7^(0^) can be obtained to three loops from 
Refs. \147\ I148j . and for completeness the result is given in App. [E| 



3.3 Operator Product Expansion 



The difference between the beam function operators in Eq. (3.4) and the PDF operators 



in Eq. (3.11) is the additional separation in the coordinate between the fields. Hence, 
by performing an operator product expansion about the limit y~ — ?■ we can expand the 
renormalized operators Oi{y~ ,uj, fi) in terms of a sum over Qi{uj',fj,), 



(3.35) 



For completeness we included the identity operator on the right-hand side. The form of the 
matching coefficients lij and Jj is again constrained by RPI-III invariance so the structure 
of the OPE is completely determined by the SCET symmetries. 



Fourier transforming both sides of Eq. (3.35) with respect to y we get 



Oi{t,uj,fi) = Ji{-t,n)l + 



j 



(3.36) 



Taking the vacuum matrix element of both sides, and using (0|Qj|0) = 0, we just get 



the coefficient of the identity operator on the right-hand side, which from Eqs. (3.31) and 



(3.32) is thus given by Ji{—t,fi). Taking the proton matrix element of Eq. (3.36) with 



a; > according to Eq. (3.10), this first term drops out, because the jet functions only have 



support for —u > (or alternatively because the corresponding diagrams are disconnected) , 
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and we obtain the OPE for the beam function 



Bi{t,x,fi) = ^I^I,j{t,'^,fi)f,{tfi) l + 0[-^) . (3.37) 



For Bg this equation was already derived in Ref. [96] using a moment-space OPE for the 
matrix element (modulo missing the mixing contribution from the quark PDF). The higher- 



order power corrections in Eq. (3.37) must scale like 1/t and are therefore of 0{AQQY)/t) 
where Aq^j-, is the typical invariant mass of the partons in the proton. (Equivalently, one 
can think of the scaling as (Aqqj-j/cj) where Aqq^/cj is the typical plus momentum of the 
parton in the proton.) They are given in terms of higher-twist proton structure functions. 
Since Eq. (3.36) is valid for t ^ Aqqj-,, this also means that we can calculate the matching 
coefficients in perturbation theory at the beam scale /i^ ~ t. This matching calculation is 
carried out in the usual way by computing convenient matrix elements of the operators on 



both sides of Eq. (3.36) and extracting the matching coefficients from the difference. This 
is carried out at tree level in the next section, while the full one-loop matching calculation 
for the quark beam function is given in Ch. |4] and for the gluon beam function in Ch. [5j 
On the other hand, for t ~ Aqqj-, the beam functions are nonperturbative and the OPE 
would require an infinite set of higher-twist proton structure functions. In this case, the 
beam functions essentially become nonperturbative ^'''-dependent PDFs. 



The physical interpretation of the beam function OPE in Eq. (3.37) leads exactly to the 



physical picture shown in Fig. 2-4 and discussed in the introduction. At the beam scale 
/is — ^1 the PDFs are evaluated and a parton j with momentum fraction ^ is taken out 
of the proton. It then undergoes further collinear interactions, which are described by the 
perturbative Wilson coefficients Iij{t, z, fj,). By emitting collinear radiation it looses some 
of its momentum, and the final momentum fraction going into the hard interaction is x < 
In addition, the sum on j indicates that there is a mixing effect from terms without large 
logarithms, e.g. the quark beam function gets contributions from the quark, gluon, and 
antiquark PDFs. For example, when an incoming gluon from the proton pair-produces, 
with the quark participating in the hard interaction and the antiquark going into the beam 
remnant, then this is a mixing of the gluon PDF into the quark beam function. These are 
the physical effects that would usually be described by the PDF evolution. The difference is 
that once we are above the beam scale these effects only cause non-logarithmic perturbative 
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corrections, which means the parton mixing and x-reshuffling now appears in the matching, 
while the RG evolution of the beam function only changes t, as we saw above. In Ch. [6] 
we will see that these matching corrections are still important numerically and must be 
taken into account. For example, since the gluon PDF at small is very large compared 
to the quark and antiquark PDFs, it still gives an important contribution to the quark and 
antiquark beam functions. 

The consistency of the RGE requires that the /i dependence of the Wilson coefficients 
Iij{t,z,fj,) turns the RG running of the PDFs into the proper RG running of the beam 



functions. Taking the fi derivative of Eq. (3.37) we find the evolution equation for the 
Wilson coefficients 



d f d^' 



(3.38) 



The solution to this RGE can be easily obtained in terms of the evolution factors for the 



PDF and beam function in Eqs. (3.18) and (3.29), 



^'t' ^'^ik{t - t', ^,Moj Usit' , f^o, fJ-) ulj{z',no,fJ') ■ (3.39) 



Hence as expected, the RGE running of Iij{t,z,ii) cancels the running of the PDFs and 
adds in the running of the beam function. 



3.4 Tree-level Matching onto PDFs 

To illustrate the application of the OPE, we will calculate the Wilson coefficients lij at 
tree level, starting with Iqq. We can use any external states for the computation of the 
Wilson coefficient as long as they have nonzero overlap with our operator. Thus, we pick 
the simplest choice, n-collinear quark and gluon states, \qn{p)) and \gnip)), with momentum 
= {p^,p^,0) where p~ > is the large momentum. In the following section we will use 
a small < as an IR regulator, but otherwise is set to zero. The tree-level diagrams 



with an external quark for the quark PDF and beam function are shown in Figs. 3-l(a) and 
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(b) 



Figure 3-1: Tree-level diagram for the quark PDF (a) and the quark beam function (b). For 
the latter, the y~ coordinate separation in the operator is indicated by drawing separated 
vertices for each field. 



3- 1(b) They give 



(9n(p)|Qg(w',/i)|g„(p))^°^ = 9{uj')Un{p)6{uj' -p ) ^Un(p) = e{uj') 6{1 - Uj' /p ) , 



{qn{p)\Og{t,U, fl)\qn{p)f^ = Un{p) 5{t) 5{uJ - p") '^Un{p) = 6{t) 6{1 - Uj/p-) . (3.40) 



Here and in the following the superscript {i) indicates the 0{a\) contribution. Note that 



the results in Eq. (3.40) are the same whether we use a state with fixed spin and color 



or whether we average over spin and color. Taking the matrix element of both sides of 



Eq. (3.36) and using Eq. (3.40), we can read off the tree- level matching coefficient 



X(°)(t, z, /x) = Xg-)(t, z, /i) = 5{t) 6{l-z). 



(3.41) 



Similarly, the tree-level results for the gluon PDF and beam function are 



{9n{p)\Qa{i^\^^)\gn{p)f^ = 9{UJ') 6{l - J/p 



{9n{p)\Og(t,u,ij)\gn{p)f^ = -oje*-e6{t)6{u-p-) = 6{t)5{l - u/p-) 



(3.42) 



leading to 

I^fit,z,^,) = 6{t)6{l-z). (3.43) 
Since at tree level the quark (gluon) matrix elements of the gluon (quark) operators vanish, 



9n{p)\Qg{^^' , fJ')\gn{p)y^^ = (^n (p) | Qg(w', ^) | (p) 



0, 



{9n{p)\Oq{t,UJ, ^i)\gn{p))^°^ = {qn{p)\Og{t, UJ , ^i)\qn{p) 



,(0) 



(3.44) 
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we obtain 

xg)(f,z,M)=x(°)(t,z,/x) = 0. (3.45) 
To summarize, the complete tree-level results are 

I^^\t,z,fi) = 5ijS{t)S{l-z), B^''\t,x,fi) = d{t)fi{x,^i). (3.46) 

The interpretation is simply that at tree level the parton taken out of the proton goes 
straight into the hard interaction. However, even at tree level the OPE already provides 
nontrivial information. Prom our general discussion we know that the matching should be 
performed at the beam scale /^(^ ~ f to avoid large logarithms in the 0{as) terms, and this 
determines the scale at which the PDFs must be evaluated to be ji = hb- 



3.5 Analytic Structure and Time- Ordered Products 

In this section we discuss the analytic structure of the beam functions. For the OPE 
matching calculation we want to calculate partonic matrix elements of Oq{t, u, (j). Por this 
purpose it is convenient to relate the matrix elements of the products of fields in Oq{t, uj, n) 
to discontinuities of matrix elements of time-ordered products of fields, since the latter are 
easily evaluated using standard Peynman rules. Por notational simplicity we only consider 
the quark operator Oq{t, uj) and suppress the spin indices and ^ dependence. The discussion 
for the antiquark and gluon operators are analogous. 

We are interested in the forward matrix element of Oq{t,u}) between some n-collinear 
state \pn) = \Pnip~^ ,P~)) with large momentum p~ and small residual momentum p'^. 
Inserting a complete set of states J2x\-^){'^\i S^^ 

{pn\Oq{t,u)\pn) = ^(pn|xn(0)| 6(1 - U;p+)\x){X\ [5{u; - ^n)Xn(O)] \pn) (3.47) 
X 

= Y.S{t- Ujp+) 6{U -p- +p^)(^Pn\xn{0)^\x){X\xn{0)\pn) ■ 
X 

The S{uj — Vn) by definition only acts on the field inside the square bracket, returning 
its minus momentum, which by momentum conservation must be equal to the difference 
of the minus momenta of the external states. Since uj = p~ — p'^, requiring a; > implies 
p'^ < p^ . This means that the action of the field reduces the momentum of the initial state 
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so it effectively annihilates a parton in the initial state . Similarly, for a; < we would 
have > P~ ^iid the field would effectively create an antiquark in {X\. Also, since \X) 
are physical states, we have > so w < and t = wpj and to have the same sign. 

Hence, for the beam function, where = \pn{P~)) is the proton state, the restriction 



to a; > in its definition, Eq. (3.10), enforces that we indeed take a quark out of the proton. 



(Note that uj < does not correspond to the anti-quark beam function.) Taking the states 
\X) to be a complete set of physical hadronic states, the beam function has the physical 
support 

< X < 1 - < 1 , t> > , (3.48) 



where Pxmin ^ ^ smallest possible hadronic momenta (which are strictly positive 

because with an incoming proton \X) can neither be massless nor the vacuum state). For 
the jet function the external state is the vacuum \p„) = |0) yielding 5{uj+p'^), so the matrix 



element in Eq. (3.47) vanishes for > 0. 



Next, consider the following time-ordered analog of (^pn\Oq{t,uj)\pn), 



(p„|Tg(a;6+,a;)|pn/ 



T{xn{y-\)\[S('^ -'Pn)Xnm] 



Pr. 



(3.49) 



Writing out the time-ordering, 



n 



0{y~)Xn {^^ - ^n)Xn(O)] - ) - T'OXnW] Xn (v" |) \ , (3.50) 



using 



9{±y-) = — Uk- , 

^ ^ ^ 27r / K + iO ' 



(3.51) 



inserting a complete set of states, and translating the fields to spacetime position zero. 



{pn\Tq{u:b'^ ,uj)\p, 

i f dy~ dn 



(27r)2 J 2\oj 



K + iO V 



Ab+-p+-K)y /2 



P +Px) 



)(pn|xn(0)||x)(X|xn(0)|p„> 



+ ei(b++4+Mj/-/2^(^+^-_^^)<p^|^^(0)|X>(x|xn(0)||p„) 
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ilm 6+ 



-Px min 



Px min 



Re 6+ 



Figure 3-2: Cuts in the complex 5+ plane for the time-ordered product in Eq. (3.52) 



27r|a; 



E 

X 



Pn 



Xn{0)^ Xj{X\xn{0)\pn) 

i 



0+ -I- p J — lO ' ' \ I 



Pr. 



(3.52) 



The first term creates a cut in the complex 5^ plane for 6+ > p\^^-a- This cut is shown as 



the dark red line in Fig. 3-2 The second term produces a cut at Ir < —Pxmin^ shown as 



the light blue line in Fig. 3-2 



The beam function matrix element in Eq. (3.47) can be identified as precisely the dis 



continuity of the first term in Eq. (3.52) with respect to b^. Thus, for the beam function 
we have 

BqiLob+,uj) = DisCb+>o {pn{P')\e{cj)Tq{ujb+ ,uj)\pn{P-)) . (3.53) 
Taking the discontinuity only for 6"*" > ensures that we only pick out the cut due to the 



first term in Eq. (3.52). Here, the discontinuity of a function g{x) for x > xq is defined as 



DisCi:>a:o g{x) = lim 9{x - xq) \gix + i/3) - g{x - i/3)l 



(3.54) 



and we used Eq. (B.7) to take the discontinuity of 1/(6^ ~ Px) 



Disc6+ 



>o 



2ttU 6^ 



Px 



6{b+-p+) = S{Lob+-cop+). 



(3.55) 



Since we explicitly specify how to take the discontinuity, we can drop the iO prescription 
in the denominators. (Alternatively, we could multiply by i and take the imaginary part 
using the iO prescription.) Since we are eventually only interested in the case a; > anyway. 
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we may as well restrict the matrix elements to a; > from the very beginning and instead 
simply take the discontinuity in t for t = ujb'^ > 0, so 

{pn\0(uj)Og{t,Uj)\pn) = DisCf>o {pn\0{uj)Tq{t, tj)\pn) ■ (3.56) 

For the matching calculation is a partonic quark or gluon state. For any contri- 
butions with real radiation in the intermediate state, i.e. diagrams where the two Xn or 
Bn± fields in the operator Oi are joined by a series of propagators and vertices, we can 
use the standard Feynman rules to evaluate the time-ordered matrix element of Tq{t,u}). 
However, with partonic external states, we can also have the vacuum state as an interme- 
diate state, because the fields in the operator are spacetime separated. For such purely 
virtual contributions it is simpler to directly start from Oq{t,uj), insert the vacuum state 
between the fields, and then use standard Feynman rules to separately compute the two 
pieces (pn|xn(0)^/2|o) and (o|xn(0)|pn)- In fact, this is exactly what we already did in 
our tree-level calculation in Sec. |3.4[ and we will see another example in Ch. |4| Thus, we 
will obtain the total partonic matrix element as 

{pn\9iuj)Og{t,Uj)\pn) = {pn\0{Lj)Oqit,u)\pn)^.^^^^y+ {pn\0iuj)Og{t,lj)\pn)^^^.^^.^^ 



5{t)5{uj-p~)(pn Xn(0)| O) (0|xn(0)|p^ 
\ LI connected 



' connected 



+ Disct>o (p„|0(u;)rg(t,aj)|p„>^^^^^^^^^ . (3.57) 

The virtual contribution must be kept, since it only looks superficially disconnected because 
the operator itself is spacetime separated. As always, we still disregard genuinely discon- 
nected diagrams, e.g. diagrams involving vacuum bubbles, when calculating the matrix 
elements in the second line. 
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Chapter 4 



The Quark Beam Function 



In this chapter we compute the matching coefficients Igqi^t, z, fi) and Igg{t, z, fi) in the OPE 
for the quark beam function to next-to- leading order in as{fj,), ffist reported in Ref. |168) . 



As explained in Sec. 3.3 and Sec. 3.4, the matching can be done by computing the 



partonic matrix elements of both sides of Eq. (3.36) to NLO. We use the same n-collinear 



quark and gluon states, \qn) = \qn{p)) and \gn) = \gn{p))i as in the tree-level matching 



in Sec. 3.4 with momentum = {p^,p~,0). Since only Iqq{t, z, fj,) is nonzero at leading 
order, we will only need the NLO matrix elements of the quark operators, Oq{t,Ljj, fi) and 
Qq{u;,fi). We will write all the matrix elements in terms of the RPI-III invariant variables 

t = h^uj , t' = —p'^oj = —zp'^p~ , z = — . (4.1) 

p 

Here, z is the partonic momentum fraction of the quark annihilated by the operator relative 
to the momentum of the incoming quark or gluon, and will coincide with the argument of 
the Wilson coefficients Xjj (t, z, /x). 

To regulate the UV we use dimensional regularization with d = 4 — 2e dimensions and 
renormalize using the MS scheme. Since the matching coefficients in the OPE must be 
IR finite, the matrix elements of Oq and Qq must have the same IR divergences, i.e., the 
beam function must contain the same IR divergences as the PDF. To explicitly check that 
this is the case, we separate the UV and IR divergences by regulating the IR with a small 
p'^ < 0. This forces the external states to have a small offshellness p^p^ < 0, and since 
p'^p~ = —t'/z the IR divergences will appear as Int'. This also allows us to directly obtain 
the one-loop renormalization constants and anomalous dimensions for Oq and Qq from their 
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Figure 4-1: Nonzero one-loop diagrams for the quark PDF. The minus momentum w enters 
the vertex through its outgoing fermion hne and leaves through its incoming fermion line. 
Diagram (c) represents the inclusion of the wave-function renormalization constant for the 
renormalized fields together with the corresponding residue factor in the LSZ formula for 
the S-matrix. Diagrams (b) and (c) have symmetric counterparts which are included in 
their computation. 



one-loop matrix elements. 

We first compute the renormalized one- loop matrix elements of the quark PDF operator 



Qq in Sec. 4.1 , This calculation of the PDF for general x using the SCET operator definition 



and with an offshellness IR regulator is quite instructive, both by itself and in comparison 

we 



to the beam function calculation, which is why we give it in some detail. In Sec. 4.2 



compute the renormalized one-loop matrix elements of the quark beam function operator 
Oq. We use these results to extract expressions for Iqq^t, z, fi) and Iqg{t, z, fi) valid to NLO 
in Sec. [Ol 

Assuming that the IR divergences in the beam function and PDF will cancel, the match- 
ing calculation can be performed more easily using dimensional regularization for both UV 
and IR. We do this as an illustrative exercise in App. |Dj which, as it should, yields the same 
result for the matching coefficients. 



4.1 Quark PDF with Offshellness IR Regulator 



We start by calculating the bare S'-matrix elements 



Q^-^(a;)|(z„(p)) , (<7n(p)|Q^'=(^)|5n(p)) , 



(4.2) 



using Feynman gauge to compute the gauge-invariant sum of all diagrams. The relevant 



one- loop diagrams are shown in Fig. 4-1 Since Qq is a local SCET operator, we can use the 



usual time-ordered Feynman rules in SCET (without any of the complications discussed in 
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Sec. 3.5 for Oq). The collinear qnqnQn vertex factor is 



ig T'^V'^ip, i) \ with £) = n'^ + ^ + - ^n^^ , (4.3) 

where and are the label momenta of the outgoing and incoming quark lines. (Because 
we have a single collinear direction the computation can also be done with QCD Feynman 
rules, still accounting for zero-bin subtractions, with the only difference being the Dirac 
algebra in the numerator of the loop integral. We checked that the final results for each 
diagram are indeed the same either way.) 



The diagram in Fig. 4-l(a)| is 



where g = g{^) is the renormalized MS coupling. The Dirac algebra for the numerator gives 

Un{p)Vn^{p,l)Vnf,{e,p)^Un{p){rf = UniPhlhhlXf.^UniP) = P~ (d " 2)^1 . (4.5) 

To compute the loop integral we write d'^i = di'^d£~d'^^'^£±/2, where i± is Euclidean, so 
£^ = — The integral is done by contour integration as follows. For £~ < all poles 
are above the axis and for £^ > all poles are below the axis, so both cases give zero. 
Hence, the i~ integration range is restricted to < < P -, where there is a double pole 
below the axis from the 1/(^^+10)^ and a single pole above the axis from the l/[(^— p)^-|-iO]. 
Taking the single pole above amounts to replacing the second denominator by 2it\/ [d.^ — p^) 
and setting = — ^/ (p~ — £~) everywhere else. After performing the contour integral 
the iO have served their purpose and can be set to zero everywhere. The £~ integral is trivial 
using the 5{l~ — uj) and turns the ^~ limits into an overall 0{uj)0[p^ — u). The remaining £±_ 
integration is done in d — 2 = 2(1 — e) Euclidean dimensions as usual. Putting everything 
together, we obtain 

= —r— a Cf 9(uj)9(p - uj) / - — -t4t 
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"-('''^-«(.Wi-.)r(.)(5^y(i-^)'-(i-.)^ 



27r 

asifJ-)CF 



t' 



1 



9{z)e{l-z){l-z) 
zvr e 



where in the last hne we expanded in e. 



In 



7 



ln(l - z)-2 



(4.6) 



In the diagram in Fig. 4- 1(b) , the gluon is annihilated by the Wilson line inside one of 
the Xn fields. The contraction with the one in Xn is oc 6{i~ — u) and the contraction with 



the one in Xn is oc S{p~ — uj). The 1/Vn in the Wilson lines [see Eq. (3.6)] contributes a 
factor l/{£~ — p~) with a relative minus sign between the two contractions. (There is also 
a diagram where the gluon connects both Wilson lines which vanishes because the Wilson 



lines only contain n-A gluons and we use Feynman gauge.) Adding Fig. 4-l(b) and its 
mirror graph, which gives an identical contribution, we get 



{qn\Qr"i^)\Qn) 



47r 
as{\i)Cp 

TT 
TT 



{2TTy {£- - p- ) + iO) [{£ - p)2 + iO] 

{i-/p-y-^ 



t' 



(l-z)i+^ 



il-£-/p-y+' 
r(2 - 6)r(-e) 
r(2 - 2e) 



[6{r - u) - 5{p- - u;)] 
(4.7) 



In the first step we used n^l^C = 2 and Un{p)^Un{p) = 2p , performed the integral by 
contours and did the £± integral as usual. The integral has the same pole structure as in 



Fig. 4- 1(a) (except that the double pole at = is now a single pole), which restricts the 



£ integral to the finite range < £ < p . Expanding Eq. (4.7) in e, using the distribution 



identity in Eq. (B.3), we get 



TT 



t' 







n 


e{z) 



(5(1 -z) + £o(l - z)z - e£i(l - z)z 



1 / vr^ 



TT 



e{z 



+ 5{l-z) 
|(^-ln-^)[£o(l-^)^ + 5(l 



£i(l - z)z 



+ 6{l-z)(2 



IT 



(4.^ 



where Cnix) = [9{x){ln"' x)/x]+ are the usual plus distributions defined in Eq. (B.2) 
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In the last step in Eq. (4.7), the t~ mtegral produces an additional 1/e pole in each 
of the two terms corresponding to real and virtual radiation from the two different Wilson 
line contractions. It comes from the singularity at l~ = , where the gluon in the loop 
becomes soft. (This soft IR divergence appears as a pole in e because the offshellness only 
regulates the collinear IR divergence here.) The soft IR divergences cancel in the sum of 



the virtual and real contributions, as can be seen explicitly in the first line of Eq. (4.8) 



where the 1/e poles in curly brackets cancel between the two terms. One can already see 



this in the I integral in Eq. (4.7), because for I = p the two 6 functions cancel so 



there is no soft divergence in the total integral. Thus, in agreement with our discussion in 



Sec. 3.1, we explicitly see that contributions from the soft region drop out in the PDF. As 
a consequence, the PDF only contains a single 1/e pole and correspondingly its RGE will 
sum single logarithms associated with this purely collinear IR divergence. 

Since the gluon in the loop is supposed to be collinear, the soft gluon region must be 



explicitly removed from the collinear loop integral, which is the condition p ^ in Eq. (3.2 ). 
For continuous loop momenta this is achieved by a zero-bin subtraction. However, since the 
soft region does not contribute to the PDF, it also does not require zero-bin subtractions 
in SCET. (If we were to include separate zero-bin subtractions for the virtual and real 
contributions, they would simply cancel each other.) We will see shortly that the situation 
for the beam function is quite different. 



The last diagram with external quarks. Fig. 4- 1(c), is 



{<ln\Q^r(^)\qnf^ = S{1 - z){Z^ - 1) = -^^^ S{1 - .)|i - In^ + l} . (4.9) 

Here we used the result for the one-loop on-shell wave-function renormalization with an 
offshellness IR regulator, which is the same in SCET and QCD. 



Adding up the results in Eqs. (4.6), (4.8), and (4.9) we obtain for the bare one-loop 
quark matrix element 

{Qn\Q^ri^)\Qnf^ = e{z)l^{l - In ^)P,,{Z) - A(l - z){l + z') 

+ 6{l - z)(^ - ^) - 9{l - z)2(l -z)\, (4.10) 
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where 



P^^{z)=C,{l-z){l + z') + -5{l-z) 



21 



[l-z 



l + z 
1 - z 



(4.11) 



is the q ^ qg splitting function, see Eq. (3.20). 



Next, we consider the matrix element of Qq between gluon states = \gn{p))- The 



only relevant diagram is shown in Fig. 4- 1(d) 



(5n|Qf'^^(o;)|<7. 



in 



Ad) 



d'^^ {-ele.)ir[V^V::%{rf{i--p-) 



(£2 + i0)2[(^-p)2 + i0] 



(4.12) 



Here e = e{p), Vn = Vn{£ — p,i) and = V^{£,i — p). Since the physical polarization 
vector is perpendicular, n • e{p) = n ■ e{p) = 0, we only need the perpendicular parts of the 
collinear vertices. The numerator then becomes 



tr 



' n ' n 



-tr 



{rfir-p-) 



+ 



l~ — p~ 



.llg^^-+SrtJ\ = 2g^^'p-[- 



in\r-p-) 



Az 
d-2 
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(4.13) 



In the last step we used that under the integral we can replace i~ = uj = zp~ 
I'^^g^ / {d — 2) . The remaining loop integral is exactly the same as in Fig. 4-l(a) 
one-loop gluon matrix element becomes 



and tj""^ = 
so the bare 



{gn\Qt-{u:)\g^r-'''^'^^'^'^^^^^ .^w.r^^V^^. ...... 



27r 
2tt 



e{z)e{l-z)T{e) 



t' 



\l - zy{l - 2z + 2z'^ - e) 



9{z) 



- - In 



Le 



7 



ln(l - z) Pqg{z) -e{l-z)\. (4.14) 



Here 



Pqg{z) = e{l- z){l-2z + 2z'^ 



(4.15) 



is the g ^ qq splitting function from Eq. (3.20) 



Note that the diagram analogous to Fig. 4-l(d) with the two gluons crossed can be 
obtained from Fig. 4-l(d) by taking p^^ — )• —p^, which takes z — )• —z. The limits resulting 



from the integral are then —1 < z < or —p < cj < 0, and since we require a; > 
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for Qq, this diagram vanishes. The diagram involving the SCET vertex with two cohinear 
gluons vanishes because the integral does not have poles on both sides of the axis. 



From the bare matrix elements in Eqs. (4.10) and (4.14) we can obtain the renormal- 
ization of Qq. Taking parton matrix elements of Eq. (3.14) and expanding to NLO, 



(<!„|cS"°(<^)k„>'" 

"E/"^ [zr(^./.)<''"l«'('^''')l''">'°' + <'°'(5-'')<'"|e'("''''W''' 

= Zfq^^\z,n) + {qn\Qq{u:,li)\qn)^^\ 



Z4«(z,/x) + (5n|Q,(a;,/i)|5n>^'\ 



(4.16) 



where we used the tree- level matrix elements in Eqs. (3.40) and (3.44) and zl-^\z,^) 



5ij (5(1 — z). The MS counter terms required to cancel the 1/e poles in the bare PDF matrix 
elements are then 



Zlq{z) = 5{l-z) + \ ^^iM^ e{z)Pqq{z) , Z^z) = \ ^^iM^ e{z)Pq,{z) . (4.17) 



Expanding Eq. (3.15) to NLO, the one-loop anomalous dimensions are obtained by 



(4.18) 



which with Eq. (4.17) yields the anomalous dimension for the quark PDF in Eq. (3.19), 



Jfqqiz,^^) = '^^^e{z)Pqqiz), ^{^{z^^.) = 9 {z) Pq,[z) . 



(4.19) 



Finally, the renormalized NLO PDF matrix elements, which we will need for the matching 



computation in Sec. 4.3 below, are 



Q,(a;, ^)|g„>^'^ = -^iM^ 0{z)[pqq{z) In ^ + £i(l - z){l + z^) 



27r 



5{l-z)[ 



7_^ 
2 3 



+ e{i - z)2{i - z) 
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(ffn|Qg(^^,/U)|5'ri 



27r 



e{z] 



t' 

In ^ + ln(l 



+ 9il-z) 



(4.20) 



4.2 Quark Beam Function with OfFshellness IR Regulator 

Next, we calculate the bare beam function 5- matrix elements, 

{qn{p)\e{u)0^^'''%t,u)\qn{p)) , {gn{p)\e{uj)0^^'''%t,u;)\gn{p)) , (4.21) 



to NLO. The corresponding one-loop diagrams are shown in Fig. 4-2, The matrix elements 



are calculated as explained in Sec. 3.5 in Eq. (3.57): For the virtual diagrams with vacuum 



intermediate state we explicitly insert the vacuum state, while for the real-emission diagrams 



we use Eq. (3.56). In the latter case, we first take the Disc, then expand in e to extract 
the UV divergences, and at last take the t' — )• limit to isolate the IR divergences into Int' 
terms. Some helpful formulas for calculating the discontinuity and taking the limit — t- 
are given in App. [B| 

For the beam function calculation the p'^ < actually plays a dual role: For the UV 



divergent piece we can treat the calculation as in SCETj, and so p+ ~ 6^ ~ X^p , which 



allows us to explicitly check the structure of the convolution in Eq. (3.22 ). The renormalized 
result contributes to the matching onto PDFs, matching from SCETj onto SCETn. In the 
matching, —p~^ <^ plays the role of the IR regulator, since we are required to use the 
same states as in the PDF calculation. We will see that the IR divergences Int' match up 
with those present in the PDF calculation, and hence drop out in the coefficients lij. 

The diagrams in Fig. |4-2| have the same Dirac and propagator structure and overall 



factors as the corresponding PDF diagrams in Fig. 4-1, so we can reuse those parts from 
the previous section. The difference compared to the PDF calculation is that for the real- 
emission diagrams, instead of doing the integral by contours, is fixed by the additional 
6 function in , and since we use time-ordered perturbation theory we must now take the 
discontinuity. This also alters the structure of the remaining i± integral, for which we 
now use Feynman parameters to combine the denominators. After carrying out the £± 
integration, we will need the following two Feynman parameter integrals 

hiA,B,e)= f\a [(1 - a)A - aB]'^-^ = ' ~ f ' , (4.22) 

Jo €{A + B) 



y-/uj,-b+ 



(X) (X) 



- i 

(a) 



y-/uj,-b+ 




P 



\P 



(d) 



y-/^?-^^ 
(X) ^ 




(b) 



y-/uj,-b+ 
I 

i-p 

(e) 



y-/uj,-b+ 
6^ ^ 




(c) 




(f) 



Figure 4-2: One-loop diagrams for the quark beam function. The minus momentum u is 
incoming at the vertex and the 6"*" momentum is outgoing. Diagram (d) denotes the wave- 
function contribution. Diagrams (b), (c), and (d) have symmetric counterparts which are 
equal to the ones shown and included in the computation. Diagram (f) and the diagram 
with the gluon connecting both vertices vanish. 



l2{A,B,e)= / da {I - a)[{l - a)A - aB] 
Jo 



-l-e 



A- 



e(l-e)(A + S)2 e{A + B) 



The first diagram, Fig. 4-2(a) , has real radiation in the final state, so we use Eq. (3.56) 



(a) 



p-{d-2)il 



9 Cf \i — - Disct>o / da 



47r 



(27r)2z 



1 ^ d^"¥x {l-a)P^ 
J (2^)'^-2 [p^ + (1 _ a)A - aS]3 



27r z 



r(l + e){e^^^?Y{l - ef [-^DisCi>o h{A, B, e) 



(4.23) 



where we abbreviated 



I — z t 

A = t + t' , B = t, A + B = - + t' ^ 

z z 



(4.24) 



Since t' > and z > 0, the only discontinuity in l2{A,B,e) for t > arises from {—B). 
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Using Eq. (B.7) to take the Disc, we obtain 



i i (—B)^'^ 

~Bisc,,ol2iA,B,e) = -Disc,>o _ ^^^^ ^ 



0{t) 



sinvre 9{B)B 



l-e 



7re(l-e) {A + B)^ ' 
sinvre [{1 - z)t]^-^z^+' 



7re(l-e) {t + zt') 



i\2 



(4.25) 



Note that there is only a discontinuity for i? > 0, so taking the discontinuity for t > 
requires (1 — z)/ z > 0, and since z > we obtain the expected limit z < 1. Since there are 



no UV divergences, we can let e — )• 0, and Eq. (4.23) becomes 



2tt 



{t + zt'Y 



(4.26) 



The above result has a collinear IR singularity for t — )• which is regulated by the nonzero 
t' . We can isolate the IR singularity using Eqs. ( [R2| ) and (IRSJ) by letting p = zt'/n^ 
while holding i = t + zt' fixed^ , 



lim 



e{t)t 



lim 



e{i-zt') 9{t-zt')zt' 



<^(t)(ln^+l) . (4.27) 



The final result for Fig. 4-2 (a) is thus 



qn\eiu;)0^g^'%t,u;)\qnf^ = eiz)9{l - z)il - z){ 



27r 



2'-'0\ 2 



J(f)(ln^ + l)}. 

(4.28) 



Next, we consider the real-emission diagram in Fig. 4-2 (b) It corresponds to the 



-OJ 



term in Eq. (4.7). Together with its mirror graph, giving an identical contribution, we obtain 

{qn\9{u:)0^^-'%t,io)\qnf^ 

^./e^^H^Y 2^ G{u}) ^. f d'^i 2p-e- 6{r-uj)6{e++b+-p+) 
2i( — : I g Cf — — Discf>o 



(4.29) 



as{n)CF 



TT 



e{z)T{l + e] 



{2ttY {i--p-){£^ + io)[{e-py + io] 

r(l + e){e"'^fiy [-;^Disct>o h{A, B, e) 

1 — Z L ZTT J 

e^-E^2 mvre e{t) 6(1 - z)z^+' 



ne t + zt' (1 - z) 



where in the second step we performed the loop integral as before, and in the last step we 



^We keep the dependence on i in our calculation as it will be useful for checking the structure of the 
renormalization in the following section. 
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used Eq. (B.7) to take the discontinuity. As for Fig. 4-2(a), the loop integral produces no 



UV divergence. However, as in the PDF calculation for Fig. 4-1 (b) , there is a soft gluon IR 
divergence at z — )• 1 or ^~ — )• p~ producing a 6{1 — z)/e IR pole when expanding the last 



factor using Eq. (B.3). In contrast to the PDF calculation, the soft gluon region must now 



be explicitly excluded from the collinear loop integral. In dimensional regularization with 
an offshellness IR regulator the relevant zero-bin integral is scaleless and vanishes. Thus, 
including the zero-bin subtraction removes the 1/e IR divergence and replaces it by an equal 
1/e UV divergence such that all 1/e poles in the final result are UV divergences. Expanding 
in e, we have 



{qn\9{u;)0^g-'%t,u)\q„ 



(b) asip.)C, 



(4.30) 



and taking the same limit as in Eq. (4.27) to isolate the IR divergences 



t'^O t + Zt' 2t'/^2-S>0 t 



^Co[^] -6{i)ln 



0{t) 



t 



Oii-zt')^ i-zt' 1 
m — 



zt^ 



t'-s>o t + zt' /x^ zt'/ti'2^o i 



6{t) - ln2 — + 



(4.31) 

^2 ' e)' 



the final result for Fig. 4-2 (b) is 



{qn\e{u;)0^g^'%t,u;)\qn) 



IT 



+ 



12^ 



6{t) In 



5(l-z)+£o(l 

2, 



z)z 



vr 

2 + "6 



6(1 



(4.32) 



For the diagram in Fig. 4-2 (c) (and its mirror diagram) we insert the vacuum interme- 



diate state between the fields in Oq as in Eq. (3.57), resulting in a one-loop virtual diagram 
involving a single field. The calculation is exactly the same as for the S{p~ — uj) term in 



Eq. (4.7) times an overall S{t), 



(g„|0(u;)O,^-^(t,a;)|g„ 



-2i 



4tt 
as{n)CF 



g^CF6{t)5{p--u) 



d'^l 2p-i- 
(27r)rf (£--p-)(£2 + i0)[(^-p)2 + iO] 

r(2 - e)r(-6) 
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IT 



6{i)5{l - z) 



1 1/ t' , ^ 
+ - 1 - In ^ + ;t In 



1.2*' 



, t' TT 

ln^ + 2 



12 



(4.33) 



In the last step we expanded in e and took the IR hmit. To be consistent we have to use the 
same IR hmit in the virtual diagrams as in the real-emission diagrams above, which simply 
turns the overall 6{t) into a 5{t), 



lim 6(t)= lim 5(i - zt') = 5(i) . 



(4.34) 



As in the PDF calculation, the UV divergence in the loop produces a r(e) and the soft 
IR divergence a r(— e). The latter is converted by the zero-bin subtraction into a UV 
divergence, producing the 1/e^ pole. The 1/e^ poles do not cancel anymore between Figs. 



4-2(b) and 4-2(c) as they did for the PDF in Fig. 4-l(b) , because the phase space of the real 
emission in Fig. 4-2(b) is now restricted by the measurement of 6"*" via the 5{l'^ + 6^ ~P^)- 



For the same reason Fig. 4-2 (a) has no UV divergence anymore, while Fig. 4- 1(a) did. 



The (l/e)lnt' terms in Eqs. (4.32) and (4.33), which are a product of UV and collinear 



IR divergences, still cancel between the real and virtual diagrams, ensuring that the UV 
renormalization is independent of the IR, as should be the case. 



The final one-loop contribution to the quark matrix element. Fig. 4-2(d) and its mirror 
diagram, comes from wave-function renormalization. 



Qs(^)Cf 



11 



47r 



6{t)5{l-z) 



, t' 



(4.35) 



Adding up the results in Eqs. (4.28), (4.32), (4.33), and (4.35), we obtain the bare beam 



function quark matrix element at one loop. 



m 



2 1 

2e/ e fi 



2^0V „2 



6{l-z) 



- S{i) [P,,(z) In ^ - 6{1 - z) - ^) + ^(1 -z){l- z) 



We now consider the beam function matrix element with external gluons. The corre- 
sponding diagrams are shown in Figs. 4-2(e) and 4-2(f)[ For Fig. 4-2(e) , which is analogous 
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to Fig. 4-1 (d) , we find 



{gn\e{u)0\-^%t,uj)\gn) 



1 

1 - z 



4z \ f d'^i ?[d{£--uj)6{£++b+-p+) 

— jDisQ>oy (^2 + io)2[(£-p)2 + iO] 



27r z V 1 — z 



2z 



27r 



Disct>o/2(A5,e) 



2tt 



{t + zt 



/\2 ■ 



(4.37) 



The loop integral and discontinuity are exactly the same as for Fig. 4-2 (a) The diagram 
in Fig. 4-2(f)| does not contribute to the quark beam function. It can be obtained from 
Eq. (4.37) by replacing p^^ — )• —p^, which takes t' — )■ —t' and z — )• —z. Doing so, the only 
contribution to the discontinuity is still from B = —(1 -|- z)t/z for B > 0, which for t > 



requires — 1 < z < 0. Hence, Fig. 4-2(f) does not contribute. Using Eq. (4.27) to take 



t' — )• in Eq. (4.37), we get the final result for the bare one-loop gluon matrix element 



{gn\eioj)0^^^^%t,u;)\gnf^ = '^A^ e{z)P,g{z)\ -C,[-^]-5{t)[\n~- + l]\. (4.38) 



zt[_ 



As for Fig. 4-2(a) , it has no UV divergences because of the measurement of , which means 



that the renormalization does not mix Oq and Og. 



4.3 Renormalization and Matching 



Using the bare matrix elements calculated in the previous section, we can extract the 
renormalization of Oq. We first take t = t + zt' ^ t \n the bare matrix elements. Then, 



expanding the quark matrix element of Eq. (3.22) to one- loop order 



{qn\OY'%t,u:)\qn) 



(1) 



dt' 



{t,n) 5(1 - z) + {qn\Oq{t,UJ,^l)\qn 



(4.39) 
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we can then read off the MS renormahzation constant from Eq. (4.36) 



2tt 



2 1 



2eJ e fj.' 



Co 



(4.40) 



The fact that the gkion matrix element is UV finite and the UV divergences in the quark 
matrix element are proportional to (5(1 — z) confirms at one loop our general result that 
the renormahzation of the beam function does not mix quarks and gluons or change the 
momentum fraction. 



In Eqs. (4.39) and (4.40) we used that we already know the structure of the renormal- 



ization from our general arguments in Sec. 3.2 i.e. that only depends on the difference 
t — t' . Alternatively, we can also use the dependence on z and the finite dependence on t' 
via t to explicitly check the structure of the renormahzation. In this case, we must use the 



same IR limit also for the tree-level result in Eq. (3.40), which using Eq. (4.34) becomes 



(g„,|0(6^)Og(t,w,;u)|(7n>^°^ = Urn 5(t)<5(l - z) = 5(i)<5(l - z) 



(4.41) 



Taking Z'L(t, t' ,uj/u;' , /j,) to be a general function of t, t' and uj/u' , we now get for Eq. (4.39 ) 



(g„|0,b-^(t,a;)|g„> 



(1) 



dt"^ zfHt, t", ^, /.) 6{t" + z't') <5f 1 - ^ 

UJ' \ OJ' / \ p~ 

Zf \t, -t',z) + {qn\Og{t,io,fi)\qny^^ . 



[qn\Oq{t,UJ,n)\qny^^ 



(4.42) 



In the first step we used Eq. and zf^\t,t',z) = 6{t - t')6{l - z). From Eq. ( [436 ) 



we now find 



ZUt,t',z,^i) = {d{t-t') + 



27r 



/X/ 2 3 \ 2 1 



Co 



t-t' 



5{l-z), (4.43) 



thus explicitly confirming at one loop that Z'j^{t,t',z,fi) = Z|(t-f,/i)(5(l 



The one- loop anomalous dimension for the quark beam function follows from Eq. (4.40), 



(4.44) 



It is identical to the one-loop anomalous dimension of the quark jet function. The coefficient 
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of £o(*//^^)/^^ can be identified as the one-loop expression for — 2rcusp- Thus, Eq. (4.44) 



exphcitly confirms the general results in Eqs. (3.26) and (3.34) at one loop 



Taking the bare matrix elements in Eqs. (4.36) and (4.38) and subtracting the UV 



divergences using Eqs. (4.39) and (4.40) gives the renormalized one-loop beam function 
matrix elements, 



{qn\0{uj)Og{t,UJ, fl)\qn] 



zf 

Pqg{z) In 



'7 vr^ 



(1) = ^^,(,)|^^,fi^>(i_,)+J^^„fi^)^,(i_,)(i+,2, 
-m 

(1) _ as{lx)TF 



5{l-z)[---]+e{l-z){l-z) 



27r 



^(^)^..(-){;;^^o(^) -<^(t)(ln^ + l)}. (4.45) 



For the matching onto the PDFs, we must take t' — )• and have therefore set t = t 
everywhere, only keeping t' in the IR divergent Int' terms. 



Expanding the OPE for the quark beam function, Eq. (3.37), to one loop, we have 



{qn\9{uj)Oq{t,U,fi)\qn)^^'^ 

= X«(t,z,^) + <5(t)(g„|Qg(a;,/x)|g„>^') , 



(1) 



(1) 



(4.46) 



Thus, the one-loop matching coefficients, T^^^ {t, z, /i), are obtained by subtracting the renor- 



qi 



malized PDF matrix elements in Eq. ( 4.20[ ) from those in Eq. (4.45). Doing so, we see that 



the Int' IR divergences in Eqs. (4.20) and (4.45) precisely cancel, as they must, such that the 



matching coefficients are independent of the IR regulator and only involve large logarithms 
that are minimized at the scale /x^ ~ t. 
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The final result for the NLO matching coefficients is given by 



Tqq{t, z, fi) = d{t) d{l - z) 



(4.47) 



+ 



27r 



+ d{t) 



3 

d 

C,{l-z){l + z^)-'^S{l-z) + e{l-z)(l-z-^:^hiz) |, 
{t, z, n) = e{z) I ^Co[^) Pqg{z) + S{t) [P,,(z) (in 1^ - l) + ^(1 - z)j I . 
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Chapter 5 

The Gluon Beam Function 



In this chapter we compute the matching coefficients for the gluon beam function at one 
loop. We determine the renormalization group equation for the gluon beam function and 
verify that its anomalous dimension equals that of the gluon jet function. These calculations 
will appear in Ref. |48j . 

The matrix element which defines the gluon beam function first appeared in Ref. [96] and 
was named the gluon beam function in Ref. |166j . In Ref. |142j the impact-parameter gluon 



beam function was introduced, which differs from the gluon beam function in Eqs. (3.4) and 
(3.10) because the Lorentz indices on the two B'^j_ fields are not contracted and the fields 



have an additional separation in the residual perpendicular coordinate y±. In Refs. j96j and 
|142j only the matching onto gluon PDFs was considered and there is a discrepancy in one 
of the terms. Our results agree with those presented in |142j and we also determine the 
matching onto quark PDFs. 

In our calculation we will regularize both the UV and the IR using dimensional regular- 
ization and employ the MS renormalization scheme. To perform the renormalization of the 
beam function and calculate the matching coefficients we replace the external proton states 



by free quarks or gluons, as was explained in Sec. 3.2 for the tree-level matching. In our 
partonic calculations we will abbreviate the matrix element notation used in the calculation 
of the quark beam function in Ch. |4]by writing and z, /u) to indicate that 

the proton is replaced by a parton of type j in the PDF or beam functions, i.e. 



An 



jn{P~ 



n 



Bi^Ay--]5{u:-Vn)BZ_{Q) 



jn{P~ 
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(5.1) 



We start in Sec. 5.1 with recalling the known results for the partonic PDF for our choice 



of regulator and scheme. In Sec. 5.2 we list the necessary Feynman rules, for which we 
introduce some shorthand notation. We calculate the one-loop gluon beam function in 
Sec. 15.31 and extract the results in Sec. 15.41 



5.1 The Gluon PDF at One Loop 



At tree level the partonic PDF is normalized to 



4°)(z,m) = 5.,<5(1-z). 



(5.2) 



The renormalization of the PDF in the MS scheme is given by 

J 3 

where we sum over parton species j = q, q, g, and the entries in the matrix Z-^ are a series 
of 1/e poles with coefficients in terms of the renormalized coupling Q!s(//). At one- loop, the 
results for the standard gluon PDF are [9], 



Zf^iz, /.) = S{1 -z) + \^e{z) [CAPggiz) + 6{l - z) 



(5.4) 



where /3o = [I'^Ca — 4?T'/77)/3, is the lowest order coefficient of the QCD /3 function. The 
splitting functions are given by 



2Co{l- z)z + 29{l- z 

z 



1 - z 



+ z(l-z] 



29{1 - z 

e{i-z) 



(i-z)+ 

1 + (1 - z)2 



1 — z 
+ + z{l-z) 



(5.5) 



We will extract the one-loop partonic gluon PDF from the above results. With our 
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p,a,oc 

(a) 



p,a,a. 



q,b,/3 



(b) 



Figure 5-1: The SCET Feynman diagrams for the gluon field strength at 0{g^) (a) and 



0{g) (b) given by the Feynman rules in Eq. (5.9) and Eq. (5.10), respectively, 



choice of dimensional regularization as IR regulator, the only Lorentz invariant quantity 
they can depend on is z. This is dimensionless, so all diagrams contributing to the PDF 
calculation vanish. Hence, the one-loop PDF has IR divergences with signs opposite to the 



UV divergences. The UV divergences follow from Eq. ( 5.4 ) by looking at the renormalization 
of the partonic gluon PDF at one-loop 



bare (1) 



g/g 



gj 



"(i.)/j;](.',M)+z5»>(i. )/«(.',.) 



/.bare (1) / n 



(5.6) 



This implies that the UV divergent part of the partonic f^^j{z,fj,) is just Z-j- {z,ij,). We 



thus obtain the renormalized PDF 



■^^^e{z)[CAPgg{z) + l(3o6{l-z) 



1 as{fJ-)CF 
e 2tt 



e{z)p,g{z). 



(5.7) 



5.2 SCET Feynman Rules 



We now list the SCET Feynman rules necessary for the calculations. The Feynman rules 
for the gluon field strength, 



(5.8) 
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are illustrated in Fig. 5-1 The Feynman rule at order 0{g^) is 



Fig.[5^= <5-(5r - ^ S^'^B^^l'-'^ip) 



n-p 



(5.9) 



and at 0{g) it is given by 



Fig.[5^T(b)]=i5/ 



cab 



+ 



n-p n-q \n-q n-pJn-{p + q) 



^gr'l3^^l'''^{p,q). 

(5.10) 



We abbreviate the collinear quark-gluon vertex as 



n-p n-q n-pn-q /I I 



(5.11) 



where p and q are the momenta of the outgoing and incoming quark lines, so the gluon 
carries incoming momentum p — q. Finally, we abbreviate the triple gluon vertex by 



gr'"'yrivx.v2.p^) = gr^w'^kvY - p2y + g'''{P2 - psr + g'^ips - piT] , (5.12) 



eabcrnu 



where all momenta are incoming and momentum conservation holds, Pi + P2 + P3 = 0. 



5.3 The Gluon Beam Function at One Loop 



We now turn to the calculation of the partonic gluon beam function. The beam function 



was defined in Eq. (3.10) such that at tree level it is normalized to 



Bl'j]{t,z,fi) = 6ij6{t)6{l-z). 



(5.13) 



The one-loop Feynman diagrams for the partonic calculation are shown in Fig. 5-2 Since 



we regulate both the UV and IR with dimensional regularization, the virtual diagrams 

vanish because there is only the p~ momentum of the external 



in Figs. 5-2(d) and 



5-2 e 



gluon flowing into the loop, which is insufficient to give a nonzero Lorentz-invariant result 
for the loop integral. This means that the UV divergences cancel the IR divergences. 
Performing the wave function renormalization in the on shell scheme, both the wave function 
renormalization as well as the residue that enters in LSZ are equal to one. (A different 
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y-/uj,-b+ 



y-/uj,-b+ 




(d) 



(e) 




Figure 5-2: One-loop diagrams for the gluon beam function. The minus momentum uj is 
incoming at the vertex and the momentum is outgoing. Graphs (b), (d), and (e) have 
symmetric counterparts which are equal to the ones shown and included in the computation. 
The diagrams with the external lines crossed vanish and are not shown. 



scheme would lead to contributions to and that cancel each other in the final result.) 
To compute the one-loop real radiation diagrams for the gluon beam function, we use 



the relationship with the matrix element of the time-ordered product in Eq. (3.56). We 



always average over the polarizations of the external gluons, and since their momentum is 
= p-fiP- 12 this gives us the rule 



d-2 

pol 



(5.14) 



where d = 4— 2e. We also average over the color of the external gluons (which is trivial). This 
choice for the polarizations and colors of the external states is just a matter of calculational 
convenience, in order to determine the renormalization and matching all we need is states 
with non-zero overlap with the operator. 



Let us start with the diagram in Fig. 5-2 (a) To make our expressions more palatable, 
we first perform the color algebra, which simply yields 



srab 



jaec jbde ^fc ^fd _ _(J^ 



(5.15) 
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The diagram is then given by 



Fig. |5-2(a) 



(5.16) 



-11—^) g^CA9iu;)Bisct>o 



9i 



^ ^ yaau^p^ ^ _ _e)Vsp,A-P, e,p-i) B'^l^ii) Ci^'^C 



?(o) ^J.p/ 



d-2 



in J (27r)2 



Disct>o / da (1 — a) 





(2vr) 



d-2 



(^1 + A)3 



27r 



2)^(1 - z)2 



z 



with A = f(l — a/ z). To derive the second equahty we integrated and /~ using the delta 
functions and applied a Feynman parametrization. The last line follows from performing 
the i±_ integral and using the discontinuities listed in App. [BI By expanding in e we find 



Fig. |5-2(a) 

~ 2-K 



(5.17) 



{z)e{i - z)2 



[1 — z ,^ , z 
— + .(l-.) + - 



m 



1 ^ l-z 
h In 

e z 



1 

7' 



Co 



where the plus distribution Co{x) is defined in Eq. (B.2). The color structure for the diagram 

gab 



in Fig. 5-2 (b) yields 



m-1 



jaec gfcjfeb ^ 



Including a factor of two for its mirror graph, we get 



47r 



(27r)^ (£2 + iO) [{£ - pY + iO] 



d-2 



1 



^' ^) (2^ —z i ^« J (2^ 



e{z){i + z)z 



1+e 



[l-z 



27r ' ' '' ' (1 - z) 

using the relations listed in App. |B| Expanding in e yields 



^r(l + 6)(6^-/x2)^ 



sinvre Q{t) 



(5.18) 



(5.19) 
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Fig. |5-2(b) 

asiiJ,)CA 



(5.20) 



+ 



2tt 
2 



1 r/ N 1 ^ / i 



£i( 4jV(1 - ^) + + ^) fA(l - 2) - £0(1 -z)lnz- ^(5(1 - z) 



12 



The diagram in Fig. 5-2 (c) is zero. Because the external gluons have perpendicular polar- 



ization, contracting the two B^]_ in Eq. (5.10) leads to n • n = in the numerator. 

Crossing the external lines of these diagrams corresponds to p — ?■ —p (the polarization 
and color are symmetric and hence unaffected). By changing i ^ —i this yields the same 
as the original diagram but with 6{£~ — uj)5{i'^ + 6+) — )• (5(^~ -|- u;)S(i^ — 6+) and thus 



A 



1 



a 



1 + - It. 

z . 



(5.21) 



Since a, z > the discontinuity for t > vanishes for these graphs. Other one-loop 
real radiation diagrams, like the one involving a four gluon vertex, vanish because their 
discontinuity vanishes. All cuts of these diagrams involve cutting one of the external lines, 
which sets A = t in the loop integral and leads to a vanishing discontinuity for t > 0. 

Adding up the two non- vanishing diagrams in Figs. 5-2(a)| and 5-2(b) , we find 



2-K 



e{z) 



2 1 



t 



e fi^ \fj, 

+ i?^'(?)'<i-'-> + i?^»(?)^"<^' 

g(l-z)lnz] 2[l-z(l-z)]^ TT^ 
z 6 



6{l-z)--5{t)Pgg{z) (5.22) 



<5(t){[A(l 



1 



6(1 



'})■ 



This expression agrees with the expression found in Ref. |142j . Transforming to moment 
space, we find agreement with the expression quoted in Ref. [96j up to an extra term of 
-asCA6{l - z)tt^/8. 

We now determine the mixing contribution of the quark PDF to the gluon beam function, 
which was not calculated in Refs. [9611142] . The relevant matrix element has external quarks 
instead of gluons, and the corresponding diagram is shown Fig. 5-2(f)[ Using fixed spin and 
color or averaging over them yields again identical results. Here, we average over spins. 
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using 5 ^spins Un{p)un{p) = \i>- The color average gives {\/Nc)ti[T°-T^]5°''^ = Cp. The 
diagram is thus given by 



Fig.F2(f) 



(5.23) 



II rCF^(w) Disct>o ' 



{iTT^ (£2+i0)2[(p-^)2 + i0] 



(2^ 



in 



d-2 



0{z)(- h^)Disct>o/ da(l-a) 

\ 1 — z z / In 



2tt 



e{z)0{l - z)T{l + e)ie'^^ fi^y 



i + a-z] 



— ez 



sinvre 6{t) 



vre VI 



This is the same loop integral and discontinuity as in Fig. 5-2 (a) For the crossed graph, 
where the external lines are interchanged, we change i ^ —i to find the same expression as 



inEq. (5.23) with S{i- - uj)5{£+ + b+) ^ 5{£- + uj)6{£+ - b+). This leads to Eq. (5.21) and 



a vanishing discontinuity. Expanding in e, we obtain for the bare one- loop quark matrix 
element 



a/i 



(5.24) 



as(/^)Cj 
27r 



- +ln^^ — -] +9(l-z)z 
e z ' 



The matrix element with external antiquarks gives the same result, so the mixing contribu- 
tions from quarks and antiquarks are identical. 



5.4 Renormalization, Running and Matching 

We know that the anomalous dimension of the gluon beam function equals that of the gluon 
jet function, and that the IR divergences of the gluon beam function cancel those of the 
gluon PDF in the matching. Since we regulate both the UV and IR using dimensional 
regularization, we cannot separate the UV and IR divergences and we can only verify one 
of these statements in our calculation. We will assume the cancellation of IR divergences, 
which then allows us to verify that the anomalous dimension equals that of the gluon jet 
function. Equivalently, we could have taken the anomalous dimension of the gluon jet 
function for granted and have checked that the IR divergences cancel in the matching. 
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We denote the counterterm that renormahzes the beam function by Z^{t — t' , fi) 



Bf^^^it, x,fi) = J dt' Z^it - t', fi)Bgit', X, /.) 



(5.25) 



In order to extract the renormahzation for the beam function from our partonic calculation, 
we expand this to one-loop order for gluon states 



B 



bare(l) 
9/9 



{t,z)= dt' 



zf\t,^,)5{l-z) + B^^l{t,z,n). 



9/9^ 



(5.26) 



Similarly expanding Eq. (3.37) for the matching onto PDFs 



9/9^ 



l(l\t,z,fi) + 6{t)fl'/^{z,f,). 



(5.27) 



We know that the IR divergences of the PDF cancel those of the beam function in 



the matching Eq. (5.27). Hence, after subtracting the IR divergences of the gluon PDF in 



Eq. (5.7) from the gluon beam function in Eqs. (5.22) and (5.24), all remaining divergences 



are ultraviolet and enter into the renormahzation. There is no UV divergence in the mixing 
graph and hence the RG evolution does not mix the quark and gluon beam functions. Using 



Eq. (5.26) we conclude that 



2tt 



11 / 1 

— 2-^0 1 ~^ 



+ ^/3o'^(t) 







vr 



+ "4^5(t) . (5.28) 



We see that our one-loop calculation agrees with the anomalous dimension of the gluon jet 



function, as it should according to Eq. (3.26). 



Since the PDFs in Eq. (5.7) contain no finite pieces, we infer from Eq. (5.27) that the 
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remaining terms are simply the matching coefficients 



as{iJ,)CA 

9(1- z)lnzi2[l- z{l- z''^^ ^2 



I,,{t, z, = + ^^^^^ e{z) ( ^£1 (l^)s{l-z) + ^Co (^)P99{Z) 



i^^it,z,i,) = ^^^l^e{z){p,,{z) 



) + 5{t) In 



1 - z' 



+ e{l- z)z6{t)) . (5.29) 
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Chapter 6 

Beam Function Plots 



We will now compare results for the beam functions at LO and NLO in fixed-order per- 
turbation theory as well as at NLL and NNLL in resummed perturbation theory, which 
were first reported in Refs. \4:8 \ 1168] . Our conventions for the Os loop counting are given 



in Table 6.1 To evaluate the required convolutions of plus distributions at NNLL we use 
the identities from App. B of Ref. [135]. We always use the MSTW2008 [HH] parton distri- 
butions at NLO for as{mz) = 0.117 and with two-loop, five-flavor running for as{f^)- The 
uncertainty bands in the plots show the perturbative uncertainties, which are estimated by 
varying the appropriate scales as explained in each case. They do not include the additional 
uncertainties from the PDFs and as{mz)- 

The order of the running of as{fj.) deserves some comment. Working consistently to 
NLO in the matching corrections requires us to use NLO PDFs, for which the two-loop 
running of was used in Ref. jl43j . On the other hand, the double- logarithmic running 
of the beam functions at NNLL requires the three-loop running of a^, which poses a slight 
dilemma. Ideally, we would need NLO PDFs using three- loop running for as(^), which as 
far as we know are not available. The numerical difference between Og run at two and three 
loops is very small, at most 2%. Hence, we use the following compromise. To be consistent 





matching 




r 

cusp 




LO 


0-loop 








NLO 


1-loop 








NLL 


0-loop 


1-loop 


2-loop 


2-loop 


NNLL 


1-loop 


2-loop 


3-loop 


3-loop 



Table 6.1: Order counting in fixed-order and resummed perturbation theory. 
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with our PDF set, we use the above as{mz) and two-loop, five-flavor running to obtain the 
numerical value of at some required scale, and to be consistent with the RGE, we use 
the two- and three-loop expression for the QCD /3 function in the RGE solutions at NLL 
and NNLL. (For simplicity we use the same NLO PDFs and also at NLL.) 

To illustrate the importance of the various contributions to the quark and gluon beam 
functions, we also consider the beam functions in the threshold limit and without the mixing 



contribution. In the threshold limit we only keep the terms in Eq. (4.47) and Eq. (5.29) 
which are singular as z 



1, 



X*^-^^t,z,M) = Smi - z) + ^fi!^e{z)l^^C^{^)5{l -z) + ^Co{^)Coil - z) 



27r 



+ 5{t) 2Ci{l-z)-—5{l-z) 





^^t, z, = - ^) + ^^^^ ^(z) I (^) 5(1 - .) + -^£o ) /:o(l - ^) 



+ 5{t) 2Ci{l - z) - —5{1 - z) 



X*^'-^t,z,M)=Xj^-^^t,z,^) = 0. 



9Q 



(6.1) 



The mixing terms Xqg and Xgq contain no threshold terms (which reflects the fact that in 
threshold Drell-Yan the gluon PDF does not contribute). For the results without mixing 
contributions we keep the full Xa but set Xij to zero for i j, which corresponds to adding 
the remaining non-threshold terms in Xa to the threshold result. In the plots below, the 
results in the threshold limit are shown by a dotted line and are labeled "x — )• 1", and the 
results without the mixing contribution are shown by a dashed line and are labeled "no 
for the quark beam functions and "no g" for the gluon beam function. The full results are 
shown by a solid line. Hence, the size of the non-threshold terms in Xa, and therefore the 
applicability of the threshold limit, is seen by the shift from the dotted to the dashed line, 
and the effect of the mixing is given by the shift from the dashed to the solid line. 

To be able to plot the beam functions as a function of the momentum fraction x including 
the virtual terms proportional to 5{t), we integrate over t up to some maximum tmax. 



Bi{t^s.x,X,n) = j dtBi{t,X, fl)9{trai,^ - t) 



(6.2) 



where Bi{t,x,fj,) is given by Eqs. (3.28) and (3.37). In the plots, we always choose t 
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Figure 6-1: The u (upper left), u (upper right) and g (lower left) beam functions at the hard 
scale HH = X 7TeV at LO, NLO, NLL, and NNLL, integrated up to W = {x e"^ TTeV)^. 
The bands show the perturbative uncertainties estimated by varying for the fixed-order 
results and the matching scale ^\ ~ tmax for the resummed results, explained in the text. 



{xe-^TTeVf, which one should think of as tmax = (e ^'^^ X-Ecm)^- Hence, this choice of 
^max corresponds to a rapidity cut y"^^^* = 2 for £^cm = 7TeV or equivalently 7/"^"* = 2.4 for 
Ecia = lOTeV. This is motivated by the upper bound y^^^ = it Y , which follows from 



the factorization theorem Eq. (2.29) when we integrate tb < exp(— 2y^*). 



Figure 6-1 shows the integrated li, u and gluon bGa^m function ^-^^(tniax; x,iih) evaluated 
at the hard scale = Q = x7 TeV. For the fixed-order results at LO (lowest gray band) 
and NLO (light green band) the perturbative uncertainties are obtained by varying fiH 
by factors of two, since this is the scale at which the perturbation series for the matching 



coefficients in Eq. (3.37) is evaluated. At LO, the resulting variation is entirely due to 
the scale dependence of the PDF. This is clearly an underestimation, since the perturbative 
corrections in the NLO beam function contain large single and double logarithms of tmax/Q^j 
causing a large shift in the NLO central value with large uncertainties. For the gluon beam 
function the shift and uncertainties are particularly large, due to differences between Igl^ 
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and Iql^ such as an overall coefficient of Ca instead oi Cp- 

For the resummed results at NLL (dark blue band) and NNLL (medium orange band) 
the beam function OPE in Eq. (3.37) is evaluated at the beam scale /i^ ~ tmaxj and the 



beam function is then evolved to fin using its RGE, Eq. (3.28). In this way, the large 
logarithms of ~ t^a.x/ fJ-'n — are resummed. The perturbative uncertainties are 

now evaluated by varying the matching scale /i^, where the perturbation series is evaluated, 
while keeping [xu fixed. The uncertainty bands show the minimum and maximum variation 
in the interval ^ft^^jl < /^_b < 2-y/tmax (which due to the double-logarithmic series do not 
occur at the edges of the interval) with the central value given by the center of the bands. 
The NLL result is close to the NLO result, showing that the large logarithms make up the 
biggest part of the NLO corrections. Consequently, the corrections from NLL to NNLL 
are within the NLL uncertainties. Resummed perturbation theory is thus well-behaved and 
should be used for the beam function at the hard scale, rather than fixed-order perturbation 
theory. The large uncertainty bands for the resummed gluon beam function at small x are 
due to the strong double logarithmic running. The effect is much larger than for the quark 
beam functions because Cp gets replaced by Cyi, which for the running shows up in the 
exponent. The effect of the running is enhanced at small x because asif^s) becomes large. 
For X = 0.01 we are varying the matching scale fiB between 5GeV and 20GeV where Os 
is much larger than for the corresponding variation for the fixed order results, where 
varies between 35 GeV and 140 GeV. 

To study the perturbative corrections to the beam functions in more detail, we consider 
them at the scale /i^ ~ imax, where there are no large logarithms and we can use fixed-order 



perturbation theory. The u and d beam functions at LO and NLO are shown in Fig. 6-2 



the u and d beam functions in Fig. 6-3 and the gluon beam function in Fig. 6-4 The 
left panels show xBi{t^i^x,x, fis)- The right panels show the same results but as relative 
corrections with respect to the LO results. At LO, the only scale variation comes from the 
PDFs and the minimum and maximum variations are obtained for fiB = {\/^max/2, 2-y/tmax} 
with the central value at = \/t^nax- For the NLO results, the maximum variation in the 
range < hb < 2V*max is approximately attained for /i^ = {0■7^/t^,2.0y/t^} 

and the corresponding central value for ^b = 1.4-y/tmax- To be consistent we use the same 
central value ^b = 1.4\/tjnax for the NLO results in the threshold limit and without the 
mixing contribution. For the quark and anti-quark beam functions the NLO perturbative 
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Figure 6-2: The u (top row) and d (bottom row) beam functions at the beam scale n\ ~ tmax 
at LO and NLO, integrated up to tmax = (xe~^7TeV)^. The left column shows the functions 
times X. The right column shows the relative differences compared to the LO result. Also 
shown are the NLO beam functions in the threshold limit (dotted) and without the gluon 
contribution (dashed). The bands show the perturbative scale uncertainties as explained in 
the text. 



corrections are of 0(10%) and exhibit reasonable uncertainties. The NLO corrections for 
the gluon beam function are 0(30%), which is much bigger but still reasonable. 



The integration limits x < < 1 in the beam function OPE, Eq. (3.37), force z = x/^ 



1 in the limit x — >• 1. Hence, the threshold terms in Eq. (6.1 ) are expected to dominate over 



the non-threshold terms at large values of x. This can be seen in Figs. 6-2 and 6-3, where 
the threshold results shown by the dotted lines approach the full results towards large x 



values where the beam functions vanish. For the quark beam functions in Fig. 6-2, away 
from the endpoint, x < 0.5, the threshold corrections give a poor approximation to the 



full NLO corrections. For the antiquark beam functions in Fig. 6-3, the threshold result 
turns out to be relatively close to the full result even for small x. However, the reason for 
this is a relatively strong cancelation between the non-threshold terms in the quark and 
gluon contributions Xqq and Xqg at one loop. As shown by the result without the gluon 
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Figure 6-3: The u (top row) and d (bottom row) beam functions at the beam scale. The 
meaning of the curves is analogous to Fig. |6-2[ 
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Figure 6-4: The gluon beam function at the beam scale. The meaning of the curves is 



analogous to Fig. 6-2 



contribution (dashed lines) the non-threshold terms in the quark and gluon contributions 
each by themselves are of the same size or larger than the threshold contributions. This 
cancellation appears to be accidental, since it depends on both the relative size of the 
antiquark and gluon PDFs as well as the relative size of the non-threshold terms in Igq 
and Iqg. Note also that for the d beam function the threshold result approaches the no- 
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gluon result rather than the full result at large x. This is even more visible for the gluon 
beam function, where the threshold terms might seem a good approximation around x ~ 0.2. 
However, for large x the threshold terms approach the no-quark result which is very different 
from the full result. 

It has been argued UHl [67] that the steep fall-off of the PDFs causes a systematic 
enhancement of the partonic threshold region z — )• 1 even away from the hadronic threshold 



limit X — )• 1. This likely explains why the threshold terms in Figs. 6-2 and 6-3 start to 
dominate already close to the x values where the PDFs are close to zero, rather than 
strictly near x = 1 |42j . However, our results show that the same arguments do not apply 
in the relevant region of x where the PDFs and beam functions are substantially nonzero. 
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Chapter 7 



Isolated Factorization Theorem 



In this chapter we derive the isolated factorization theorem in Eq. (2.16), which was first 
presented in Ref. |166j . Our analysis is based on factorization in SCET, which rigorously 
and systematically separates hard, soft and collinear contributions [301 [351 [37] . We make use 
of a setup with SCETj and SCETn [36], carrying out the factorization in two stages at the 
scales Cf' and UafiB^h respectively. We have an SCETj analysis to factorize initial-state jets 
from soft radiation. The initial-state jets described by beam functions in SCETj are then 
matched onto initial-state PDFs with lower offshellness for the collinear particles in SCETn- 
In this chapter, we carry out the SCETj computation, while the matching onto SCETn was 
discussed in Ch. [3j Our analysis below uses similar tools as used in the derivation of the 
factorization theorem for hemisphere invariant masses for e~^e~ — >• dijets in Ref. [93j, but 
differs significantly due to the kinematics, and the fact that we have initial-state rather 
than final-state jets and a further matching onto SCETn. The soft dynamics of e"*"e~ — )• 
dijets was studied earlier in SCET in Refs. [33l|33j. We start with a brief overview of the 



necessary SCET ingredients in Sec. 7.1 and describe the relevant kinematics in Sec. 7.2 We 
derive the factorization theorem for isolated pp — )• XL in Sec. |7.3[ including arguments to 

we 



rule out contributions from so-called Glauber degrees of freedom. Finally in Sec. 7.4 
apply the factorization theorem to pp — )■ X£~^i~ and quote final results for the beam thrust 
cross section with one-loop corrections and logarithmic resummation. 
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7.1 SCET 



Soft-collinear effective theory is an effective field theory of QCD that describes the interac- 
tions of colhnear and soft particles [291 Ell 133 [37] . CoUinear particles are characterized by 
having large energy and small invariant mass. To separate the large and small momentum 
components, it is convenient to use light-cone coordinates. We define two light-cone vectors 

n'' = (l,n), n^ = (l,-n), (7.1) 

with = = 0, n ■ n = 2, and n is a unit three-vector. Any four-momentum p can then 
be decomposed as 

= n-p— + n-p— +p^^. (7.2) 

Choosing n close to the direction of a collinear particle, its momentum p scales as {n ■ p,n- 
PiPn±) ^ n ■ p 1, A), with A ^ 1 a small parameter. For example, for a jet of collinear 
particles in the n direction with total momentum px, n-px — 2Ex corresponds to the large 
energy of the jet, while n ■ px ~ p\/Ex ^ Ex, so A^ ~ p\/E\ <^ 1. 

To construct the fields of the effective theory, the momentum is written as 

pi^=pt^ + k'' = n-p—+p^^ + k'' (7.3) 

where n ■ p ^ Q and pn± ~ XQ are the large momentum components, where Q is the scale 
of the hard interaction, while k ~ X^Q is a small residual momentum. The effective theory 
expansion is in powers of the small parameter A. 

The SCET fields for n-collinear quarks and gluons, ^n,p(a;) and An,p{x), are labeled by 
the collinear direction n and their large momentum p. They are written in position space 
with respect to the residual momentum and in momentum space with respect to the large 
momentum components. Frequently, we will only keep the label n denoting the collinear 
direction, while the momentum labels are summed over and suppressed. Derivatives acting 
on the fields pick out the residual momentum dependence, id'^ ~ ~ X^Q. The large label 
momentum is obtained from the momentum operator Vn, e.g. Vn5,n,p = P^ (,n,p- If there 
are several fields, Vn returns the sum of the label momenta of all n-collinear fields. For 
convenience, we define Vn = n ■ Vn, which picks out the large minus component. 

Collinear operators are constructed out of products of fields and Wilson lines that are 
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invariant under collinear gauge transformations [31^ [37] . The smallest building blocks are 
collinearly gauge-invariant quark and gluon fields, defined as 



5{u-Vn)Wi{x)Ux) 

5{lO + Vn)Wi{x)iD'^.Wn{x) 



9 L 



where 



(7.4) 



(7.5) 



is the collinear covariant derivative and 



WJx) 



^ ^'y^vi^ n- An{x) 

perms 



(7.6) 



The label operators in Eqs. (7.4) and (7.6) only act inside the square brackets. Here, 



Wn{x) is a Wilson line of n-coUinear gluons in label momentum space. It sums up arbitrary 
emissions of n-collinear gluons from an n-collinear quark or gluon, which are 0{\) in the 
power counting. Since Wn{x) is localized with respect to the residual position x, we can 
treat Xn,i^{x) and Bn,w{x) as local quark and gluon fields. The label momentum oj is treated 



as a continuous variable, which is why we use a 5- function operator in Eq. (7.4). It is 



set equal to the sum of the minus label momenta of all fields that the 5 function acts on, 
including those in the Wilson lines, while the label momenta of the individual fields are 
summed over. 

In general, the effective theory can contain several collinear sectors, each containing 
collinear fields along a different collinear direction. To have a well-defined power expansion 
in this case, the different collinear directions rii have to be well separated |30j . 



Ui-rij ^> A 



for 



(7.7) 



which is simply the requirement that different collinear sectors are distinct and do not 
overlap. For pp — )• Xl^l~ , we need two collinear sectors, Ua and n;,, along the directions of 
the two beams. We use a bar to denote the conjugate lightlike vector, so nj • = 2. As the 
beams are back to back, we have Ua ~ fn,, so • ~ 2 and Eq. (7.7) is easily satisfied. 



Particles that exchange large momentum of 0{Q) between collinear particles moving in 
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different directions have to be off shell by an amount of 0(nj • rijQ'^). These modes can be 
integrated out of the theory at the hard scale Q by matching full QCD onto SCET, which 
yields the hard function. The effective theory below the scale Q then splits into several 
distinct collinear sectors, where particles in the same collinear sector can still interact with 
each other, while at leading order in the power counting particles from different collinear 
sectors can only interact by the exchange of soft particles. This means that before and after 
the hard interaction takes place, the jets described by the different collinear sectors evolve 
independently from each other with only soft but no hard interactions between them. 

The soft degrees of freedom, responsible for the radiation between collinear jets, are 
described in the effective theory by soft^ quark and gluon fields, qs{x) and As(x), which 
only have residual soft momentum dependence id^ ~ X^Q. They couple to the collinear 
sectors via the soft covariant derivative 

iD^ = id>^ + gA'^ (7.8) 

acting on the collinear fields. At leading order in A, n-collinear particles only couple to the 
n ■ As component of soft gluons, so the leading-order n-collinear Lagrangian only depends 
on n ■ Dg. For n-collinear quarks [3H [37] 




The leading-order n-collinear Lagrangian for gluons is given in Ref. |35j . 

The coupling of soft gluons to collinear particles can be removed at leading order by 
defining new collinear fields [35] 

X^nl(x) = Yi{x)Xn,M)^ (7.10) 

= Y^i^) Yn{x) = (x) 

where Yn{x) and ynix) are soft Wilson lines in the fundamental and adjoint representations, 

T^y:'{x) = Yn{x)T''Yi{x). (7.11) 

^In some situations it is necessary to distinguish two types of soft sectors, referred to as soft and ultrasoft 
in the SCET hterature. In this paper we only need what are usuaUy caUed uhrasoft particles, so we will 
simply refer to these as soft. 



y„(x) = Pexp 



15 



ds n-As{x + s n) 
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Figure 7-1: Definition of the different collinear momenta related to the incoming beams. 
The soft radiation is not shown. 



The symbol P in Eq. (7.11) denotes the path ordering of the color generators along the 



integration path. The integral limits in Eq. (7.11) with the reference point at — oo are the 



natural choice for incoming particles [TQj- The final results are always independent of the 
choice of reference point, and with the above choice the interpolating fields for the incoming 
proton states do not introduce additional Wilson lines |21j . 



After the field redefinition in Eq. (7.10), the leading-order SCET Lagrangian separates 



into the sum of independent nj-collinear and soft Lagrangians, 

/:SCET = 2]4?+^s + -- - , (7.12) 

rii 

with no interactions between any of the collinear and soft sectors. The ellipses denote terms 
that are subleading in the power counting. This decoupling is what will allow us to factorize 



the cross section into separate beam and soft functions. The field redefinition in Eq. (7.10) 



introduces soft Wilson lines in the operators, which because of Eq. (7.12) can be factored 
out of the matrix element and will make up the soft function. 



7.2 Kinematics 

Before deriving the factorization theorem, we discuss the relevant kinematics, as illustrated 



in Fig. 7-1 As already mentioned, we introduce a separate set of collinear fields for each of 
the beams, with the light-cone vectors and rif, aligned with the beam directions. To derive 
the factorization theorem we work in the center-of-mass frame of the hadronic collision, so 
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the momenta of the incoming protons are (neglecting the proton mass) 

Pa' = E,^'f, Pi: = E,^'^, (7.13) 

with fia = —fif,. In particular, nj, = fia and na-rih = 2. We will mostly keep the dependence 
on the two beam directions explicit, but one should keep in mind that Ua and are related. 

The collinear fields in the and n;, directions describe the interactions within each of 
the beams before and after the collision, and are also responsible for initiating the hard 
interaction. We define the momenta of the spacelike off-shell partons that go into the hard 
interaction as 

P'a = XaE^n. Y ' Y ' ' = ^^^^^ ^ " ^ " ' (7-14) 

where Xa and Xb are the light-cone momentum fractions at which the beam functions will 
be evaluated. The power-counting parameters for the collinear sectors are 



where the relevant momenta are those of the off-shell partons in Eq. (7.14), because these 
are the momenta carried by the n^- and n^-collinear fields. 

We write the momentum of the incoming partons that are taken out of the proton as 

Caii^cm Y + O(Aqcd) , ^^^cm ^ + O(Aqcd) , (7.16) 

which defines the light-cone momentum fractions ^a,b at which the PDFs are evaluated. The 
typical _L-momenta of partons in the proton are ©(Aqcd)) while the small plus components 
are ©(Aq^j-j/E'cm)- These momenta are much smaller than any soft or residual momenta in 
SCETi and are expanded, which precisely corresponds to the OPE for the beam functions 



in Eq. (2.31) when matching them onto SCETn. 



The momentum of the final-state remnant of the proton is thus given by 

= (1 - e,)i?e.n y , (7.17) 
while the remnant of the initial-state jet radiated into the final state by the beam function 
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has momentum 



(7.18) 



and similarly for the rih direction. The total na-collinear momentum in the final state is the 



sum of Eqs. (7.17) and (7.18), or equivalently, the difference between the proton momentum 



and Eq. (7.14), 



a± 



(7.19) 



In addition to the collinear momenta, we define kg as the total four-momentum of the 
soft radiation in the final state. Hence, the total hadronic momentum in the final state is 
given by 

P'x = (Pa' - P'a) + {Pt - Pt) + ' (7-20) 

and we can write total momentum conservation Pa + Pj^ = Px + as 



P^+p'^ = q' + k^, 



(7.21) 



where is the total leptonic momentum. 

The collinear and soft momenta, ba, b^, ks are not experimentally measurable quantities. 
Instead, the experiments can only measure hadronic quantities, such as the hemisphere 



momenta = Ua-Ba and Bt = rih-Bb introduced in Sec. 2.1 Splitting the total soft 



2-2 



momentum into its contributions from each hemisphere, kg = ka + k^ as shown in Fig. 
we then have 

Bii = b^a+< + kii, B^ = b^,+r^ + k^, (7.22 
and defining k^ = Ua ■ ka, k'^ = Ub ■ kf,, we get 



Bt = Ua-Ba = bt + kt, B+ = nyBb = b+ + k+ . 



(7.23) 



In particular, the remnant momenta r^^^ do not contribute to B'^j^. A physical argument 
for this was discussed in Sec. 12.11 

Next, we decompose the total leptonic momentum as 



n - na a 



(7.24) 
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where contains the two components of q^^ transverse to the beam direction. Taking the 
z-axis along the fta beam direction, we have 

5^ = 9°T9., g!^ = (0,gT,0), (7.25) 

where qx = {qx, Qy) is a two-vector in the transverse x-y-plane. The total leptonic invariant 
mass and rapidity are 

q"" =q+q- +ql = q+q- -q^, F = J In ^ , (7.26) 

z q 

with 

= e^^ ^Jq^ + q^, d^g = ^ dg+dg" d^T = ^ dg^ dF d^T ■ (7.27) 



As we will see in the next section, the derivation of the factorization theorem requires 
us to be insensitive to the transverse components qT such that we can freely integrate over 
them. Therefore, we have to expand the kinematics in the limit qx = 0. This expansion is 



justified because from Eq. (7.21) we have 



(1± 



-Pxi. 



XQ. 



(7.28) 



A parametrically large q'^ ^ Q would require a separate jet at large px ^ Q to balance the 
transverse momentum, which is not allowed in our setup. The kinematics of the hard matrix 
element in the factorization theorem is then given by the tree- level partonic kinematics, with 
the partonic momentum conservation 



+ 



(7.29) 



which implies 



X altera — Q 
XbEcm = 



r~2 Y 
re , 



/~2 -Y 

'q^e , 



Y 



~q = XaXbE^ 



2 g^ 



1 1 Xa 

- m — 

2 Xb 



(7.30) 



Equations (7.28) and (7.29) imply that parametrically the leptons are back to back in the 
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transverse plane. Since and q can differ substantially, the leptons do not need to be 
back to back in three dimensions. 



7.3 Derivation of Isolated Factorization Theorem 

We now proceed to derive the isolated factorization theorem for generic processes pp — )• XL, 
where the hadronic final state X has a restriction on the hemisphere momenta B^^^. The 
derivation is carried out using SCET without Glauber degrees of freedom. The proof that 
Glauber effects are not required is given at the end of this section. 



7.3.1 Cross Section in QCD 

We will generically refer to properties of L as "leptonic", even though L can contain any 
non-strongly interacting particles. We only consider processes where the hard interaction 
couples the strong and electroweak sectors through one two-particle QCD current. (This 
includes for example Drell-Yan or Higgs production through gluon fusion with the Higgs 
decaying non-hadronically, but does not include electroweak Higgs production via vector- 
boson fusion.) Then, at leading order in the electroweak interactions, we can factorize the 
full-theory matrix element into its leptonic and hadronic parts 

Mipp^ XL) = '^Lj{X\j\pp) . (7.31) 
J 

The sum runs over all relevant color-singlet two-particle QCD currents J, and Lj con- 
tains the corresponding electroweak matrix element, including the electroweak propagator 
coupling to J. For example, for Drell-Yan with L = (.^i" , the relevant currents are 

Ji}f = qfrqf, J'^f = q/ri^qf , (7.32) 



so in this case the sum over J in Eq. (7.31 ) includes the sums over the two Dirac structures. 



the vector index fi, and the quark flavor / = {u, d, . . .}. The corresponding Lyjr and 



are given below in Eq. (7.71). 



The cross section for some hadronic observable O in the center-of-mass frame of the 
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collision, averaged over proton spins, is 



^ ld<^L{27rf6U.q-PL) (7.33) 



dq^dYdO 2El^ J 2(27r 



i E T.MpP ^ ^^)|' '^[O - foiX)] {2n)^5\Pa + Pb-q-px). 

spins X 

Here, Pa^b are the incoming proton momenta, px and pi are the total hadronic and leptonic 
momenta, d^L denotes the leptonic phase space, and the phase-space integrations for the 
hadronic final states are included in the sum over X. The last 6 function is overall momen- 
tum conservation. The function fo{X) inside the second 5 function returns the value of 
the hadronic observable O for a given hadronic state X , so the 5 function picks out all final 
states that contribute to a certain value of O. The 5^{q—pL) under the leptonic phase-space 
integral defines the measured q as the total leptonic momentum. Expanding this 5 function 



for qT = 0, the leptonic part does not depend on qx at leading order, and using Eq. (7.31), 



we can rewrite Eq. (7.33) as 

dcj 



l^Y.Ljj,{q\Y)Wjr{q\Y,0). (7.34) 



The leptonic tensor is defined as 

LjAq\Y)=jd^LL\Lj.{2^f5\q-'^+q+^-pi)l, (7.35) 

where q^ = \fq^€^^ . The hadronic tensor contains the square of the hadronic matrix 
element 



WjAci\Y,0) = j;(pp|jt(o)|x>(X|j'(0)|pp> 

X {27r)^6\Pa + Pb-q-px) 6[0 - fo{X)] , (7.36) 



where as in Ch. [3] we keep the average over proton spins implicit in the matrix element. 
Since Wjj' is integrated over qx, it can only depend on and Y, as well as the hadronic 
observable O. 

We are interested in the hadronic observables = Ua-Ba and B^ = rib-Bb. The 
hemisphere hadronic momenta i?^^(X) can be obtained from the states \X) using the 
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hemisphere momentum operators pf^ ^ 

= B^^{X)\X) , = B^iX)\X) . (7.37) 

A field-theoretic definition of in terms of the energy-momentum tensor of the field 
theory was given in Ref. ^32j . The hadronic tensor for O = {B^, 

Wjj,{q\Y,Bt,B+) (7.38) 
1^ |d4xe-''-^(pp|jt(x)|X>(X|j'(0)|pp>J[i3+-n,.i3,(X)]5[i?+-n,.i?,(X)] 



2( 

dx+dx" 



X 

^_i(q+,-+g-,+ )/2 



J^[x-'^+x+'^)5{Bt-naPa)5{Bt-nb-pb)J\^)\pp) 



(47r)2 

In the first line we used momentum conservation to shift the position of j'^ , and in the second 



line we performed the integral over which sets xt to zero. We also used Eq. (7.37) to 
eliminate the explicit dependence on X, allowing us to carry out the sum over all states X. 
The restriction on the states X is now implicit through the operator 5 functions inside the 
matrix element. 

7.3.2 Matching QCD onto SCET 

In the next step, we match the QCD currents J onto SCET currents by integrating out 
fluctuations at the hard scale Q. At leading order in the power counting, the matching 
takes the form 

J{x) =Y, [d^i du;2 e-(^^+^^)- lY^Cf^^ih, h) 0,t (6i, 52; x) 

+ C%g{biM)Ogg^,{biM-,x)\ , (7.39) 

where a, /3 are spinor indices, /i, v are vector indices, and the sum over q runs over all 
quark flavors {u,d, . . .}. The Wilson coefficients and operators depend on the large label 
momenta 

6^ = ^ , 6^ = ..2 f . (7.40) 
They will eventually be set to either q~na/2 or q^n^j^ by momentum conservation, but at 



this point are unspecified, and the sums and integrals over m, n2 and wi, u;2 in Eq. (7.39) 
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run over all sets of distinct collinear directions and large label momenta. On the right-hand 



side of Eq. (7.39), the full x dependence of the current is separated into the x dependence ap- 
pearing in the overall phase factor with large label momenta and the residual x dependence 
of the SCET operators. 

The SCET operators Og^{x) and Ogg{x) are constructed out of the collinear fields in 



Eq. (7.4). At leading order in the power counting they contain one field for each collinear di- 
rection. Since the QCD currents are color singlets, the leading operators that can contribute 
are 

(7.41) 

where j and c are color indices in the fundamental and adjoint representations. We included 
appropriate minus signs on the labels, such that we always have loi^2 > for incoming 
particles. Here, x = Xq '^^ ^ quark field of flavor q, which for simplicity we keep implicit 
in our notation. Note that the entire spin and flavor structure of the current J is hidden 



in the label J on the matching coefficients in Eq. (7.39). The gluon operator is symmetric 
under interchanging both /i o and 61 62, so its matching coefficient must have the 
same symmetry, 

cZ,ib2M) = c[;;g{h,h) . (7.42) 

We define the conjugate quark operator and matching coefficient with the usual factors of 
7°, i.e., 

Ojf (61, 62) = £U{x) xZ-^M ' = b'clg,{K ^2)7T" • (7.43) 

The matching coefficients are obtained by computing the renormalized matrix elements 



(0|...|gg) and {Q\...\gg) on both sides of Eq. (7.39) and comparing the results. In pure 
dimensional regularization for UV and IR divergences all loop graphs in SCET are scaleless 
and vanish, which means the UV and IR divergences in the bare matrix elements precisely 



cancel each other. The renormalized matrix elements of the right-hand side of Eq. (7.39) are 
then given by their tree- level expressions plus pure 1/e IR divergences, which cancel against 
those of the full-theory matrix elements {Q\J\qq) and (0| Jj^^f) of the left-hand side. Hence, 
the matching coefficients in MS are given in terms of the IR-finite parts of the renormalized 
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full-theory matrix elements computed in pure dimensional regularization. 



7.3.3 Soft-Collinear Factorization 



The field redefinitions in Eq. (7.10) introduce soft Wilson lines into the operators in Eq. ( 7.41 ) , 



(7.44) 



The time ordering is required to ensure the proper ordering of the soft gluon fields inside 
the Wilson lines. It only affects the ordering of the field operators, while the ordering of 
the color generators is still determined by the (anti)path ordering of the Wilson lines. In 
the remainder, we use these redefined fields and drop the (0) superscript for convenience. 



Since the momentum operator is linear in the Lagrangian, Eq. (7.12) allows us to write 



the hemisphere momentum operators as the sum of independent operators acting in the 
separate collinear and soft sectors. 



Pa = Pa,na + Pa^n^ + Pa,s , Ph = Pb,na + Pfe.rifc + Ph,s ■ 



(7.45) 



The Ha {rib) collinear sector cannot contribute momentum in the Ub (ra^) hemisphere. Thus, 



Pa, 



Pb,na = 0, while pa,na = Pua and pb. 



P 



nt 



reduce to the total momentum 



operators for each of the collinear sectors. For the soft sector, the distinction between the 
two hemisphere operators is important. We can now write 

5{Bt - na-Pa) = j d6+ dkt 6{B+ - 6+ - 5{bt - n,-p„J 5(fc+ - n,-p,,,) , 

6(3+ - Ub-Pb) = j d6fe+ dA:+ 5{B+ - 6+ - k+) 6{b+ - Ub-pnJ S{k+ - nb-pb,s) ■ (7.46) 



Using Eq. (7.39) in the hadronic tensor in Eq. (7.38), the forward matrix element of the 



product of currents turns into the forward matrix element of the product of the operators 



in Eq. (7.44). Since the Lagrangian in Eq. (7.12) contains no interactions between the 



collinear and soft sectors after the field redefinition, we can use Eq. (7.46) to factorize the 



resulting matrix element into a product of independent na-collinear, n;,-collinear, and soft 
matrix elements. 
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We first look at the contribution from Oqq. Tlie x integral of the forward matrix element 



of Oqq becomes 



dx~'"d2: 



PriaPnt 

dx^dx 



^-i{q+x-+q-x+)/2 gi{6i+62)-a; 



Ojf (x) 6{B^ - na-Pa) 5{B^ - nypt) O,"/(0) 



PriaP'. 



rib 



(47r)2 

X /da;adw6e'('^'"^^+'^''^")/2 



X (0 



T[Ylix) Yn,{x)]'''6{k: - na-Va,s) d{k, 



n,-p,,,)r[y4(o)y„,(o)] 



Pnb) + ( 



(7.47) 



Here, \pna) ^-^d are the proton states with momenta P^^ = Ecmn^^hf^ as in Eq. (7.13). 



The two terms in brackets in Eq. (7.47) arise from the different ways of matching up the 



fields with the external proton states. The restriction to have positive labels uj requires 
the fields in Oqq to be matched with the incoming proton states and the fields in Oqq with 
the outgoing proton states. In principle, there are two more ways to match the fields and 
external states, yielding matrix elements with the structure {p\xx\p) and {p\xx\p)j which 
vanish due to quark flavor number conservation in QCD. For the same reason, in the full 
product (^^ Ogg)(^q; Oq'q') Only the flavor-diagonal term with q = q' survives. 



We abbreviate the collinear and soft matrix elements in the last three lines of Eq. (7.47) 
as M^^{x~), M^j(x+), Ms(x+,x^). The collinear matrix elements only depend on one 
light-cone coordinate because the label momenta uJa,b are defined to be continuous. We 
could have also started with discrete label momenta, ojafi, and then convert to continuous 
labels by absorbing the residual dependence as follows: 

e'^-^^l"^ McSx'^,x-) = Y^ jdk- e'('^"+''""^'^^/^M-^^^- (x") 

■''a 

'dwae'"'""''/'M^„(x-), (7.48) 
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and analogously for M^j^(x"'"). In the second step we used that by reparameterization 
invariance the Fourier-transformed matrix element can only depend on uja + = uja- 

As an aside, note that in the well-studied case where the collinear matrix elements are 
between vacuum states, giving rise to jet functions, the distinction between discrete and 
continuous labels is not as relevant. In that case, the SCET Feynman rules imply that 
the collinear matrix elements do not depend on the residual k" (and A;_|_) components, 
and therefore the label momenta can be treated in either way. In our case, momentum 
conservation with the external state forces the collinear matrix elements to depend on 

. Therefore, the only way to eliminate the residual k~ dependence is to absorb it into 
continuous uj labels. One can easily see this already at tree level. Replacing the proton 
states by quark states with momentum p = p + Pr-, we get 



■^ik~-^+>^^-')/^q{p)\^^{x+,x-)5^^^Xn{Q)\q{p)) 



(ix~^ dx 
(47r)2 

= uud^p- 6{k~ — Pr) 5{k^ — pt) ) 



(7.49) 



and the label and residual minus momenta are combined using 5^^^- 5{k~ —p~) = 5{uj—p~). 
Our continuous oj is physical and corresponds to the momentum fraction of the quark in 
the proton. 



Returning to our discussion, to perform the x integral in Eq. (7.47), we take the residual 
Fourier transforms of the matrix elements. 



M^(x",x+) 



dk^ 

dkt dk 



(2vr) 



2 jikt^++kIx+)/2 M,{kt, k~) . 



dk- 



(7.50) 



Just as x^ , the residual momenta k^ and kf- here are all defined with respect to the common 



n = Ua- The x integral in Eq. (7.47) now becomes 



/ 



da; "'"da 
~(4^ 



^ ^-(^'■^ ^'^^^^ 

X 6{iOa - q~ +k- +k~) d{ujb - ^+ + A;+ + kf) 
diiOa - q-) Siojb - g+) M^„(0) M^,(0) M,(0) . 



(7.51) 
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In the last step we expanded — — kf = q^[l + 0(A^)]. The remaining residual 
integrations are then simply the Fourier transforms of the matrix elements at x = 0. 



Note that without the integration over qx in the hadronic tensor Eq. (7.38), the currents 



would depend on x± , which would require us to include perpendicular components 6a,fe± in 
the label momenta, and the soft matrix element would depend on x±, too. (The residual 
k± dependence in the collinear matrix elements can again be absorbed into continuous 



ba,b±-) The corresponding x± integration in Eq. (7.51) would yield an additional 6 function 
^'^{ba± + bb± + qr — ks±). Integrating over qx effectively eliminates this 6 function, which 
would otherwise force us to introduce an explicit dependence on bafi± in the beam functions. 
If one considers the qx spectrum of the dileptons for q^ q'^i our analysis here provides a 
starting point but requires further study. One cannot just use px in place of with our 
arguments to impose an analogous restriction on the final state, because at O^a^) one can 
have two jets at high pj- that still have small total px- 



The ria-collinear matrix element now reduces to the quark beam functions defined in 



Eq. (3.10), 



Jl^'l^ §k'k . ^ — lia \ 



4 N, 



Air 



X (p 



Pria 



4 Nr 



(7.52) 



We abbreviated = • and in the last step we used e'^"^ /^Xri(0)e '^"^ 
Xn{y^n/2) and Pn\Pn) = 0. Similarly, for the antiquark beam function we have 



M^JO) = ^(a;,)^^,, x'^lix) S{b+ - n,-p,) 5{u:, - P„J ^n,' (0) 



rib 



4 A^. 



■UJbBq{uJbbt,UJb/Pb ) • 



(7.53) 



Since the collinear matrix elements are color diagonal, the soft matrix element reduces 
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to an overall color-singlet trace, which defines the qq incoming hemisphere soft function, 



Stuni.K,H) = ^ (0|T[y„t^(0) YnMW-naVa,s) 5{k+-nm,,)T[Yl{Q) Y^M] |0> . 

(7.54) 

The trace is over color and the factor of l/N^ is included by convention, such that at tree 
level we have S'^^^^^^{k'^ ,k^) = 6{k^) S{k^). The soft matrix element in the second term 



of Eq. (7.47) with a ^ h interchanged is equal to the above one due to charge conjugation 
invariance of QCD. Under charge conjugation, the Wilson lines transform as C~'^YnC = 
T\Yn'^^\- The explicit time ordering is required because the fields in Yn are time-ordered 
by default, and charge conjugation only changes the ordering of the color generators but 
not of the field operators. For us this is not relevant, because the ordering of the fields is 
determined by the overall (anti-)time ordering in the matrix element. Thus, for the soft 
matrix element with a f-)- 6 interchanged, we find 

ti{Q\T[YlYn^]5{kt - na-pa,s) 6{k+ - n,.p5,s) T[y„t^y„,] |0> (7.55) 
^ tr(0|T[y„^y„tT] _ _ nyp,^,)T[YlY^^] |0> = Sf^l^^ik^ , k+) , 

where the transpose refers to the color indices. In the last step we used tvlA^ B'^ D'^] = 
tr[BADC] and the fact that the fields in Yn^ and Yn^ are spacelike separated and thus 
commute. Under parity, we have P~^y„^P = and P^^n^ • Pa,sP = n-h • pb s- Therefore, 
CP invariance implies that S^^^^^ is symmetric in its arguments, 

SZmM^K) = ^J^{0\^YlYn^]6ikt -nh-pb,s)6{k+ -na-Pa,s)T[Yl^^^ 

= SZm-M,ki). (7.56) 



Having worked out the different terms in Eq. (7.47), we are ready to include the remain- 



ing pieces from Eqs. (7.38) and (7.39). The qq contribution to W becomes 



Wjj>,^{q^Y,B+,B', 



du}adu}bS{u;a - q )S{u;b-q~^) 



/ dwi duj2 duj[ duj2 C^j'^qibiM) C",^-(6'i, 62) 



ni ,n2 ,71.;^ ,71-2 
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4 4 iV, 
xSZ^,ikt,k^) + {a^b) 



(7.57) 



All label sums and integrations from Eq. (7.39) eliminate the label 5's from Eq. (7.47). In 
the second step we defined 



ba ~ XaEcm j 



W = XfjEc 



n 



Ecm Eq^ 



/ q'^ e 
~E 



b -u — 1,111 2 



Ecm Eqyh 



'g^ e 



(7.58) 



as in Eq. (7.30), and introduced the hard functions 



HjJ'qq{ba, h) 



1 1 



trs 



Hjj'gg{ba,bb) = Hjj/gg{bb,ba) 



^ Cjqqih, ba) y Cj'qq(bb, ba 



(7.59) 



Equation (7.57) is the final factorized result for the Oqq contribution to the hadronic tensor. 



Repeating the same steps for Ogg, we obtain for the forward matrix element 



/ 



dx'^'ix -i(q+x-+q-x+)/2 i{bi+b2)-x 
(4^)2 

X (pnaPn^ 



Olfix) 6{B+ - na-Pa,s) 5{B+ - nb-Pb,s) O^'(O) 
duja dujb^ii^a - q~) 5{uJb - q'^) 
X I dbt db+ dkt dk+ 6{B+ -bt- kt) 6{B+ - b+ - k+) 

X [^nma H^l - ^a) SniUi 5(^2 " Wfc) + (a O 6)] 
X Kin, H^l - ^a) Sn'^nt ^(^2 " Wfe) + (a O 6)] 



PriaPrit 



X (0 



T[yl (0) yn, (0)] 5{k+ - na-Pa,s) s{k+ - nb-pb,s) T [yl (0) yn^ (0)] 

X u:a e{Ua) (Pn„ |C±(0) '^(^a " ^a-p„J 5(a;„ - PnJ <l(0)|p„, 

X ub e{ujb) (pn, S::fx(0) 5(6+ - nfe.j5„J 5(u;fe - P„J ^^(O) 



(7.60) 
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where we already performed the integral over x. The four terms in the third line correspond 
to the four different ways to match up the gluon fields with the incoming proton states. 
The coUinear matrix elements reduce to the gluon beam function defined in Eq. (3.10), 



UJaBg{u}abt,(^a/Pa ) 



Pna 



2 iV2-l 



(7.61) 



Including the color traces from the beam functions, the soft matrix element defines the 
gluonic incoming hemisphere soft function, 



qgg (u+ 

1 



(7.62) 



iV2-l 



treoior { T [yi^ {Q)yn, (0)] <5(A:+ - na ■pa,s) S{k+ - n, •pt.s) T [yl {0)yn^ (0)] } 



where the normalization is again convention. Putting everything together, the gluon con- 
tribution to the hadronic tensor becomes 



Wjj^ggiq', Y, B+, B+) = Hjj,ggiba, h) / dA:+dA:+ q'Bg[xaE,^{B+ - k+),Xa] 



X Bg[xtE,^iB^ - k+),xt] Sfl^^ikt, k+) , (7.63) 



with the hard function 



HjJ'gg{ba, bb 



1 1 



7, {9±uu' 9Luu' + 9±iiu' 9Lua') ClZ(&a, ^b) C'i,''(ba, bb) . (7.64) 



iV2 _ 1 2 



.flW y± Vll' ) ^Jgg 



Here we have used the symmetry of the Wilson coefficients in Eq. ( 7.42 ) to simplify the four 



terms that arise from interchanging a -f-)- 5 in Eq. (7.60). 



To obtain the full result for the hadronic tensor all we have to do now is to add up the 
contributions from the different quark flavors and the gluon. 



Wjj.{q\Y, 5+, S+) = ^ Wjj,gg{q^ Y, 5+, S+) + Wjj,gg{q\Y, B+ , B+) . (7.65) 
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Inserting this back into Eq. (7.34), the final result for the factorized cross section becomes 



da 

dq'^dYdBtdBil 



H,,{q\Y) J dk+ dk^ q^B,[xaE,^{Bt - xj 
X B,[x,EUBt - kt), xt] Sg^^,{k+, k+) , 



(7.66) 



with XahEc 



-•"^e^^ as in Eqs. (|7.30|) and (|7.58|) and the hard function 



J, J' 

f Ha Ub 



(7.67) 



The sum in Eq. (7.66) runs over parton species ij = {gg,uu,uu,dd, . . .}, where Bi is the 
beam function for parton i in beam a and Bj for parton j in beam b. Equation ( 7.66 ) is the 



final factorization theorem for the isolated pp — )■ XL and pp — )• XL processes. In Sec. 7.4 



below we will apply it to the case of Drell-Yan, which will yield Eq. (2.16) 



The beam functions in Eq. (7.66) are universal and take into account collinear radiation 



for isolated processes with x away from one. Since the soft function only depends on the 
color representation, but not on the specific quark flavor, there are only two independent 



soft functions S?^^^-^ and Sf^^^^. In the sum over ij in Eq. (7.66), there are no mixed terms 
with ij corresponding to beam functions of two different quark flavors. Likewise, there are 



no mixed terms with quark and gluon beam functions. For example, a graph like Fig. 2-5(e) 
is part of the ij = qq term in the sum. Thus, cross terms between quark and gluon PDFs 
only appear via the contributions of different PDFs to a given beam function, as shown in 



Eq. (2.31). 



The only process dependence in Eq. (7.66) arises through the hard functions Hij{q'^, Y), 



and one can study any desired leptonic observables by inserting the appropriate projections 
in the leptonic phase-space integrations inside Ljji(q'^,Y). Since the hard function Hjjuj 
corresponds to the partonic matrix element {ij\J^O){0\J'\ij) and Ljji is given by the square 
of the relevant electroweak matrix elements L^jLji, Hij[q^,Y) can be determined from 
calculations of the partonic cross section ij — t- L. Furthermore, Hij{q'^,Y) is identical to 
the hard function in threshold factorization theorems and hence in many cases is known 
from existing computations. 
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Figure 7-2: RGE running including potential Glauber modes. 
7.3.4 Cancellation of Glauber Gluons 

In the above derivation we have implicitly assumed that contributions from Glauber gluons 
cancel in the final cross section, so that we do not need Glauber interactions in the effective 
theory. To complete the proof of factorization, we now argue that this is indeed the case. 
In principle, Glauber interactions add an additional term to the SCET Lagrangian 



5 Xna ) Xnj, ) 

A,) . (7.68) 

Glauber interactions in SCET have been considered in Refs. [90^1119] . but wc will not require 
an explicit construction of Co here. Our arguments will be based on the one hand, on the 
consistency with processes where it has been proven that Glauber interactions cancel, and 
on the other hand on systematic scale separation in the language of effective field theory. 
The scale separation is valid independently of whether it leads to a factorization into simple 
matrix elements, or whether it leads to a non-factorizable matrix element with complicated 
dynamics. 

The possible danger of the Glauber modes comes from the fact that they couple the two 
collinear sectors and nt, with momentum scaling Q(X'^ ,1, X) and (5(1, A^, A). With Cq, 
there will still be interactions between soft and collinear modes present in the Lagrangian 
even after the field redefinition, so we cannot a priori factorize the full matrix element into 
independent soft and collinear matrix elements. Therefore, we have to revisit each step in 
our derivation with Cq in mind. 

Our argument will be divided into three steps: (i) above the scale fiB, (ii) at the scale 
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fiB, and (iii) below the scale ^b- For (i) and (ii) we have to consider Glauber modes with 
momentum scaling Q{X'^ , X'^ , X), which we call Gi modes. Since they have virtuality ~ 
they are integrated out at this scale. Any residual effects of Glauber interactions below ^b 
could occur from modes with momentum scaling (5(A^, A^, A^), which we call G2. These 
modes are illustrated in Fig. |7-2[ 

Above the hard scale fin — Q, we have full QCD and no distinction between different 
modes is required, so in step (i) we are concerned with contributions of Gi in the region 
fj,H > fji'> fji-B- At the scale ^h-, we integrate out hard modes with virtualities or higher 
by matching the QCD currents onto SCET currents. For our process, the leading operators 



are given in Eq. (7.41), which contain only one field per collinear direction. For the theory 
in Eq. (7.68), all other possible operators are power suppressed. The matching onto these 
currents is valid at an operator level and can be performed with quark and gluon states. It 
is independent of the hadronic matrix element we are going to take later on. The key point 
is that the exact same matching calculation and resulting Wilson coefficients G occur for 
threshold Drell-Yan and e~^e~ — )• 2 jets. For these cases it is known [771 EH] that Gi modes 
do not affect the matching of the hard function H ^ |Cp at or the running of H in the 
region > fJ- > fJ-B shown in Figs. |2-6[b) or |2-6|^c). The hard function H gives a complete 



description of the physics down to the scale fiB whether or not the modes in the SCET 



matrix elements factorize further. In Fig. 7-2, this corresponds to taking the scale fi = fiB- 



Therefore, the Gi modes can give neither large ln{fiB/ fJ-ii) terms nor finite contributions 
above hb- 

In step (ii), we integrate out modes with virtualities Q^X^ at the scale ^b, which may 
involve matrix elements with Gi modes exchanged. This matching affects the na-collinear, 
nb-collinear, and Gi modes, whose momentum scaling below /j^b changes to Q{X'^,1, X^), 
(5(1,A^,A^), and (5(A^, A^, A^), respectively. Here we consider the theory right above fiB 
including Gi modes, leaving the discussion of the theory just below /ib and G2 modes to 
step (iii). Thus, we have to consider the matrix element of the composite operator 

[XnaXnb]{x'^,X')S{Bt - na-Pa)S{B^ - nb-pb)6{uJa - V Ua) K^b - V Uf,) [Xn^Xn,] (0) , 

(7.69) 



where we suppressed all spin and color indices for simplicity, and these collinear fields still 
couple to soft fields in their Lagrangians. Since //^ is a perturbative scale, we can carry 
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out the matching onto the theory below /i^ at the operator level and do not yet have to 
consider proton states. Since the Glauber gluons are spacelike, they cannot cross the final- 
state cut indicated by the 6 functions and only appear in virtual subdiagrams. We can 
therefore make a correspondence with the calculation in step (i) as follows. For any given 
final state with collinear and soft particles, the SCET computation for (ii) is identical to 
the SCET computation carried out for the matching in step (i) but using this particular 
choice of external states.^ Since that SCET computation cannot induce any dependence on 
Gi in step (i), there can also be no contributions from Gi for the forward matrix-element 
computation here. The result of the step (ii) matching is thus given by a Wilson coefficient 
times an operator of the form 

J dktdk+C{x+,x-,B+ - kt,B+ - A;+)x'„„(0)T[y„t^y„,]x;,(0)<5(A:+ -n,.p,) 

X 5{k+ - nyp,) 6{u;a - VnJ 5{lo, - PnJXn,(0) T[y4y„JxL(0) ' (7-70) 

where the primed collinear fields have scaling (5(A^,1,A^) and (5(1, A'^, A^), and the soft 
fields in the Y Wilson lines have scaling Q(A^, A^, A^). 



For step (iii) below /i^, we have to consider the {pp\ ■ ■ -{pp) matrix element of Eq. (7.70) 
and possible contributions from G2 Glauber gluons, which can now also connect to spectator 
lines in the proton (which are primed collinear modes). The G2 gluons may spoil the 
factorization of the two collinear sectors. To argue that this is not the case, we rely heavily 
on the original proof of the cancellation of Glauber gluons for inclusive Drell- Yan in Ref . [78] . 



By construction, for our observables the k^f^ variables in Eq. (7.70) are of 0{QX'^) and 
thus only get contributions from the soft gluons. Hence, we are fully inclusive in the 
Hilbert space of the primed collinear fields. Therefore, the G2 modes as well as possible 
"ultrasoft" (3(A^, A^, A^) gluons cancel in the sum over states, just as in the inclusive case. 
This discussion for the cancellation of G2 modes is identical to Ref. [22], where arguments 
were presented for the cancellation of G2 gluons up to the scale induced by the measurement 
on the final state, which in our case is fiB- 

Physically, one could imagine that Glauber modes kick the spectators in the proton 
remnant such that they can contribute to -B^^. The above arguments show that this is 



^In practice one would never make such a complicated choice, but if one does, it must give the same 
result as picking a minimal state for the matching. 
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not the case, so that our treatment of the proton and its remnant m the derivation of the 
factorization is correct. 

Note that the above arguments do not suffice to show that Glauber interactions cancel 
when there are additional hard central jets in the final state. 



7.4 Final Results for Drell-Yan 

In this section, we present the final results for the isolated Drell-Yan cross section. Our 
discussion is split into four parts: the leptonic tensor, the hard function, the soft function, 
and the final cross section for beam thrust. 



7.4.1 The Leptonic Tensor 



To give an explicit example, we now apply the final factorization result in Eq. (7.66) to 
the Drell-Yan process with L = The relevant QCD currents are the vector and 

axial-vector currents Jj^j with h = {V,A}, already given in Eq. (7.32). The corresponding 
leptonic contributions are 



47racm 



Vf_ 

1 — m% 



cp' 1 — m^/q- 



-u{p2)j'^{vi - aa5)v{pi) 



(7.71) 



where in this section pi = p£+ and p2 = Pi- are the lepton momenta, Qf is the quark charge 
(in units of |e|), and V£j, a^j are the standard vector and axial couplings of the leptons and 



quarks of fiavor / to the Z boson. We will include the width of the Z later in Eq. (7.81) 
The leptonic phase space integral is 



/ 



d^pi d^p2 
(2vr)2 



5(pi)5{pi)5Hq 



+ 



1 



2 32^2 j 1 + cosh Ay ' 



dAy dip 



(7.72) 



where (p is the azimuthal angle of the leptons in the transverse plane and Ay is the rapidity 
difference of the two leptons: 



1, nb-pi 

yi = In , Ay = yi - y2 ■ 

2 na-pi 



(7.73) 



Since we expanded qx = 0, the leptons are back to back in the transverse plane, which 



138 



implies that at the order we are working 



P2 Pi 



Pit 



-P2T ■ 



Y = -{yi+ y2) 



-2 _ -Q 

Pit — P2T 



2(1 + cosh Ay) ' 



(7.74) 



Thus, the leptonic kinematics is described by the four independent variables {q^, Y, Ay, (p}, 
with {Y,Ay} being equivalent to {yi,y2}- For simplicity, we assume that we do not dis- 
tinguish the two leptons, as one would for example by measuring their rapidities yi or 
transverse momenta piT- We can then integrate over < ip < 2tt and — oo < Ay < oo in 



Eq. (7.72), giving an overall factor of 47r. The leptonic tensor, Eq. (7.35), now becomes 



1 



dAy dip 



327r2 J 1 + cosh Ay 

spins 

3g2 



X] LhfiP^^P^)^h'f'iPUP2) 



(7.75) 



where 



Lvvff'{<f) 
LAVff'{q^) 



Q^Q^^ _ {QfVf +VfQf')vi ^ VfVf,{vj + a] 



1 — m\ /q^ 



(1 — m^/q^Y 



afaf>{vj + af] 
(1 — m|/g2)2 

1 — m| / g2 
LvAf'fio^) ■ 



Vf{vj + aj) 



(7.76) 



7.4.2 The Hard Function 



Using parity and charge conjugation invariance of QCD, the matching coefficients for the 
vector and axial-vector QCD currents can be written as 

C'Af,,iba,h) = CaMq') h'^l,r^ 

CTjba, h) = CAgiq^) (ba + ^^'\.baK ■ (7-77) 
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By Lorentz invariance (or reparameterization invariance of Uafi and fia^b |139] ). the scalar 
coefficients can only depend on ba ■ bb = XaXbE^^ = (f ■ In principle, parity and charge 
conjugation would also allow the Dirac structures (ba — bbY (^"^ and (ba — bb)^ (75)"^- How- 
ever, as the vector and axial-vector currents are chiral even and the matching from QCD 
conserves chirality for massless quarks, these cannot be generated. For the gluon operator. 



the symmetry of the Wilson coefficient [see Eq. (7.42)] requires it to be proportional to 



1^ = ba b^. Current conservation for the vector current requires QfiCvfqq ~ ^' '^^ich 
eliminates this term. Thus, as expected, the only contribution for the gluon operator is due 



to the axial anomaly, coming from the diagram in Fig. 2-5 (b) , Since we neglect the lepton 



masses, q^iL'^f = 0, and thus CAfgg does not survive the contraction of the leptonic and 
hadronic tensors for L = i'^i^ . Hence, the gluon beam functions do not contribute to Drell- 
Yan, and the gluon PDF only appears through its contribution to the quark beam functions. 



Inserting Eq. (7.77) into the general expression for the hard function in Eq. (7.59), we obtain 



Hhh'ff'gqiba,bb) 



2Nr 



C*hf,{q')Ch'Mq^) (for hh' = {VV, AA}) 



(for hh'= {VA,AV}). 

(7.78) 



At one loop, the vector and axial- vector coefficients are equal and diagonal in flavor and 
the SCET matching computation was performed in Refs. |33 | I138j . in agreement with the 
one-loop form factors 



Cvfgiq^)=CAfgiq^)=Sf,C(q^), 



C(q\fi) = l + 



2 . , as{fl)CF 



Air 



ln^( ' o ) +31n 



iO\ . /-(7^-iO\ TT 



fi^ J \ fi^ J 6 



21 



(7.79) 



The vector current coefficient at two loops was obtained in Refs. |41t I118| from the known 
two- loop quark form factor \100\ I130| 11441 I145j . Starting at three loops, it can have a 
contribution that is not diagonal in flavor, i.e., is not proportional to 6fq. The axial-vector 
coefficient can also receive additional diagonal and nondiagonal contributions starting at 
two loops from the axial anomaly |5H I122| I123j . The anomaly contributions cancel in the 
final result in the sum over / as long as one sums over massless quark doublets. Therefore, 
they will cancel when the hard matching scale is much larger than the top-quark mass, in 
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which case the top quark can be treated as massless. On the other hand, they have to be 
taken into account when the matching scale is below the top-quark mass, in which case the 
top quark is integrated out during the matching step and its mass cannot be neglected. 



Combining Eqs. (|7.78[) and (|7.79D with Eqs. (|7.75[) and (|7.76D, the coefficients Hij{q^,Y) 



in Eq. (7.67) become 



= C^oi/gg(g^ /i) , 



(7.80) 



where at one loop 



'2 ivl + al)ive+aj)-'^QqVqVi{^-ml/q^) 

(l-m|^2)2 + ^|r|/g4 



(7.81) 

\Ciq',f^)\\ 



with |C(g^,/x)p given by Eq. (7.79), and where we also included the nonzero width of the 
Z. The RGE for the hard function Hqq{q'^,fi) is 



^dif^,^2^ ^^^^^2^^)^^_(^2^^)^ 7/^(9',^) = 2reuspK(^)]ln^+7HK(/i)], (7.82) 



where Ecusp is the universal cusp anomalous dimension [127], and the one-loop non-cusp 



term is 7//[as(;u)] = — 3as(/i) Cf/t^ |138j . The solution of Eq. (7.82) has the standard form 



2 Tj ill ^) 

Hqq{q\fi)=Hqq{q^fio)UH{q^fio,f^), UH{q^^lo,^^)=e''»^^'^'>'U^y" , (7.83) 



where JCH(/io,/u) and riH{fio,fj,) are 



4 Ecusp (ttg 



dal 



7T +7//(as 



/■"s(m) das 

mif^O, ^) = 2 / -r-A- Ecusp(as) • 



(7.84) 



Together, Eqs. (7.83) and (7.84) sum the large logarithms occurring in isolated Drell-Yan 



between the scales fin and fiB- Electroweak corrections to the hard function Hqq{q^ , fj) can 
be included using the results of Refs. [7TH73] . 
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7.4.3 The qq Soft Function 



The incoming hemisphere soft fmiction contains incoming Wilson hnes stretching from — oo 
to along Ha and Uf,. Under time reversal, each incoming Wilson line transforms into a 
corresponding outgoing Wilson line stretching from to oo along the opposite direction, 



T-iy„,T = Pexp 



/•oo 

-ig / dsrib-Asisnb) 
Jo 



, (7.85) 



where P denotes anti-path ordering. Since T itself does not affect the original ordering of 
the field operators, time ordering turns into anti-time ordering and vice versa. In addition 
Tua • Pa,s = Ub • Pb^s- Therefore, time-reversal invariance implies 

SfLnikt,kt) = ^^tr{0\T[YlYn^]6{kt -nb-h,s)S{k^ -na-Pa,s)T^^^ 

= l-tr{0\T[YlYr,,]5{kt-na-Pa,s)S{k+-nb-pb,s)T[^^^^ (7.86) 

In the second step, the complex conjugation has no effect since the matrix element is 
real, and we used parity to switch Tib^a back to n^ b- For comparison, the hemisphere soft 
function with outgoing Wilson appearing in the double-differential hemisphere invariant- 
mass distribution in e+e" 2 jets [93l [Ml HIS UMl IISZ] is 

SlLM, K) = ]^ (0|T[i;t^y„,] 5{kt - na-pa,s) 5{kt - nb-pb,s)T[YlYn^] |o) . (7.87) 



This is almost the same as Eq. (7.86), the only difference being the opposite time ordering. 
Thus, S'ihemi and S'hemi are equal at one loop, where the time ordering is still irrelevant. 
Beyond one loop, S'ihemi and S'hemi may in general be different. However, since the beam 
and jet functions have the same anomalous dimension, the combined anomalous dimension 
of the hard and beam functions in isolated Drell-Yan agrees with that of the hard and jet 
functions for the e"'"e~ hemisphere invariant-mass distribution. The consistency of the RGE 
in both cases then requires that S'ihemi and Shemi have the same anomalous dimension to 
all orders in perturbation theory. In addition, the purely virtual contributions, obtained by 
inserting the vacuum state, are the same in both cases, 

S^^ikt.H) = j^Jikt)5{kt)iT\{^\T[YlY^,^^^^^^ (7.88) 
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Hence, 5'ihemi and ^hemi can only differ by finite real-emission corrections at each order in 
perturbation theory. 

Using the one-loop results for S^l^^ from Refs. [9ll I157j , we have 



as{^i) Cf 



47r 



e{k+/fi)ln{k+/fi) 



+ 



(7.89) 



The plus distribution is defined in Eq. (B.l ). The one- loop soft function for beam thrust in 



Eq. ([2331) then becomes 5B(A:+,/i) = Sik+) + 2 5"°°P(fc+). 



7.4.4 Final Cross Section for Beam Thrust 



The differential cross section for beam thrust in Eq. (2.29) including the RGE running is 



da 



dg2 dydrs 



(7.90) 



tg + tb 



(JQY^Hij{q^,^XH)UH{q^■,^J.H,^J's) j dta dtbQ Sb{Qtb 

X j dt'^Bi{ta - t'^,Xa,fJ-B) U sit'^, HB , IJ'S) j dt'f^Bj{tb - tb,Xb,HB) UBitb,fJ'B,f^s) ■ 



We now consider its fixed-order Og expansion. To our knowledge da/dq^dYdTB has not been 
considered in perturbation theory in full QCD even at one loop. To obtain an expression for 
da/dq'^dYdTB at NLO in Ug and leading order in the power counting, we drop the evolution 
factors Uff and Ub and expand all functions to NLO at a common scale fj,. From the above 
NLO results for the hard and soft functions and the NLO results for the beam functions 
from Ch. [3j we find 



da 



dg2 dYdTB 



'^2 ivf + aj)ivj + aj)-2QiViVi{l-ml/q^) 
(l-m|/g2)2 + ^|r|/g4 



dCa d^b ^ fXa Xb 2 



, q , rB,n) fi/aiia, fJ-) fj/biCb, IJ) ■ 



(7.91) 



Here, fi/aiCa, fJ') and fj/biCb, fJ-) are the PDFs for parton i in proton a and parton j in 
(anti-)proton b. At tree level, the nonzero coefficients are 



Clf^ {Za,Zb,q^,TB,fJ-) = CfJ"" {Za,Zb,q^,TB,n) = S{tb)5{1- Za)5{l- Zb) 



(7.92) 
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At one loop, we obtain 



asin) Cf 

e{TB) 



5{l-Za)e{zb)< 



{ 


-2 




3 










+ "2 





S{tb){4 



IT 



+ 



TB 



+ S{tb) In ■ 



(1 - zt. 



1- Zb 

+ {Za O Zb 



+ S{tb) 



[l-Zb)Hl-Zb) 
I - Zb 



S{1 - Zb) 

(1 + z!) 



/ 1 + ^2 

C5g°°''(2a, Zb, q^,TB, /i) = Cg5°°P(Za, Zb, g^ Tg, ^) , 



6ii-za)eizb)e{i-zb) 



0{tb) 

TB 



+ 6{tb) In 



[^6 + (1 - ^6) 



^6 



2-K 
+ S{tb) 

C|g°°^(^a, 2^6, g^ Tb, fl) = Cgg°°P{Za, Zb, g^ Tfi, /i) , 

Ci^°°P(z„ Zb, q\TB, fi) = Cf^'^'iza, Zb, q\ TB,fi) = Cil°°P(z,, Z^, q\ TB,f^) . (7.93) 

The single logarithms of are multiplied by the QCD splitting kernels and are re- 

summed by the PDFs. Thus, in fixed-order perturbation theory the PDFs should be evalu- 
ated at the hard scale fi = Q, such that there are no large logarithms when integrating over 
< ^ 1- However, if the integration is restricted to tb < t^"* ^ 1, the plus distributions 
in Tb produce large logarithms In^ rj^^* and In r^^* , which make a fixed-order expansion un- 
reliable. These are precisely the logarithms that are resummed by the combined RGE of 



hard, jet, and soft functions in Eq. (7.90) 
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Chapter 8 

Beam Thrust Cross Section for 
Drell-Yan 



In this chapter we present and discuss results for the Drell-Yan beam thrust cross section, 



to appear in Ref. |169j . Beam thrust was defined in Eq. (2.25) and may be written as 



TB = ^Y.\P>'T\e-\'"~'^\ , (8.1) 

k 

where and Y are the dilepton invariant mass and rapidity. The sum runs over all 
(pseudo) particles in the final state except the two signal leptons, where \pkT\ and rjk are the 
measured transverse momenta and rapidities with respect to the beam axis (for simplicity 
we took all particles to be massless). Beam thrust is the analog of thrust for e~^e~ — )• jets 
with the thrust axis fixed to the proton beam axis. For <C 1, the hadronic final state 
consists of two back-to-back jets along the beam axis due to the initial state radiation. 
The measurement of tb can provide a test of our understanding of the initial state at the 
Tevatron and LHC. 

Experimentally, beam thrust is one of the simplest hadronic observables at a hadron 
collider. It requires no jet algorithms, is boost invariant along the beam axis, and can be 
directly compared to theory predictions without utilizing parton showering or hadronization 
from Monte Carlos. 

Beam thrust is also theoretically clean. It is infrared safe, and simple enough to include 
in theoretical calculations. In Sec. 17.31 we derived an all-orders factorization theorem for 
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the Drell-Yan beam thrust cross section for smaU tb 



da 
dQdYdTB 



Y^Hij{Q'^,flH)UHiQ'^,f^H,fJ') j dtadt'^Bi{ta-t'^,Xa,fJ-B)UBit'a,fJ'B,fJ') 



where Xa = {Q / Ecm)^ and = (Q/i^cm)^"^ are the partonic momentum fractions trans- 
ferred to the leptons, and the sum runs over quark flavors ij = {nn, mm, dd, . . .}. The hard 
function Hij{Q^ , fin) contains hard virtual radiation at the hard scale hh — Q (and also 
includes the leptonic decay). The soft function SB{k,fis) encodes the effect of soft virtual 
and real radiation on the measurement of tb and is sensitive to the soft scale fis — tbQ- 
The b QdbUi functions Bii^taj Xa,fiB) and Bj{th,Xb, fiB) describe extracting the colliding par- 
tons out of the proton and the formation of the initial state jets. As we showed in Sec. |3.3[ 
they can be matched onto PDFs Bi = ^^.lik ® fk at the beam scale fj-B — \/ta,b — \'^Q- 
The cross section for small ^ 1 contains large double logarithms a" In™ tb with 



m < 2n. These are summed to all orders in Eq. (8.2) by evaluating the hard, beam, and 
soft functions at their natural scale where they contain no large logarithms and using the 
RG evolution kernels Uh-, U^b ■• ^'^^ to run them to a common scale ^. The combination 
of the evolution kernels is independent and sums the logarithms of tb- At NNLL we 
need the next-to-leading order (NLO) expressions for Hij^ lik, and Sb, as well as the NNLL 
expressions for Uh, U^b ^ and Us, given in Sec. 7.4 and App. [Ej The required convolutions 



in Eq. are carried out analytically following Ref. JT35\. We use the MSTW2008 [Tl3] 
parton distribution functions at NLO for as{mz) = 0.117 and with two- loop, five- flavor 
running for as(/i). Details on our order counting, as well as a technical point related to the 
running of as(/x), can be found at the beginning of Ch. |6j 

To estimate the perturbative uncertainties we use the minimum and maximum variation 
under the following separate scale variations 

a) fJ-H = rQ, HB = r^fr^Q , y,s = ttbQ , 

b) i^H = Q, f^B = T'^^'^^^^I^^Q , ^5 = tbQ , 

c) ^XH = Q, [iB = V^Q , fis = t-^'^^^^I^tbQ , (8.3) 
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Figure 8-1: Drell-Yan cross section da/dQ at NNLL with a cut on beam thrust tb < 
exp(— 2y^*) with y^* = 1,2. The sohd hues show the cross section at the LHC with 
Ecia = 7TeV and the dashed hues the cross section at the Tevatron. 



with r = {1/2, 2}, and r = 1 corresponds to the central value. Only considering the simul- 
taneous variation in a) produces unnaturally small scale uncertainties due to cancelations 
between the running of the hard, beam, and soft functions. The exponent of r for cases b) 
and c) is chosen such that for tb = e.~* the scales ^b or /x^ vary by factors of two, with 
smaller variations for increasing tb and no variation for — )■ 1. In this limit, there should 
only be a single scale = /^_B = l^s^ ^-ncl thus the only scale variation should be case a). 
For the cross section integrated up to tb < t^*^^, the scales are chosen by replacing tb with 



r]j^^=^ in Eq. ([8^). 

As already mentioned, small tb describes events with energetic initial-state radiation 
in the forward direction but no hard jets at central rapidities. More explicitly, a cut on 
Tb < exp(— 2y^*) vetos energetic radiation with total energy > Q in the rapidity region 
I?/ — y| < y^*, while allowing energetic radiation for |y — y| > y^* with a smooth transition 



between the two regions. In Fig. 8-1 we show the Drell-Yan cross section da/dQ at NNLL for 
ygit = 1 and 2 at the LHC with Ecm = 7TeV and the Tevatron. Increasing the restriction 
on the hadronic final state from y^* = 1 to 2 reduces the cross section by a factor of few 
at larger Q which becomes a factor of 10 at small Q. The Z resonance is visible around 
Q = 90 GeV. 

Figure 8-2 shows the same cross section for = 1 normalized relative to the NLL 
result at the LHC with Ecm = 7 TeV (left panel) and the Tevatron (right panel) . The bands 



147 




Q [GeV] Q [GeV] 



Figure 8-2: The Drell-Yan cross section da/dQ at NNLL with a cut on beam thrust 
TB < ex.p{-2yf^) with yf^ = 1 at the LHC with = 7TeV (left panel) and the 
Tevatron (right panel). All results are normalized relative to the NLL result. Shown are 
the full NNLL result (solid), NNLL without the gluon contribution (dashed), and NNLL in 
the threshold limit (dotted). The bands show the scale uncertainties at NLL and NNLL as 
explained in the text. 



show the perturbative uncertainties at NLL (medium blue band) and NNLL (dark orange 
band). From NLL to NNLL they are reduced by about a factor of two. As illustration, the 
dashed line shows the NNLL result without the contribution of the gluon PDF to the quark 
beam function, Iqg in Eq. (3.37). The dotted line shows the NNLL result in the threshold 
limit X —7- 1, where in addition to neglecting the gluon contribution we only keep the 



leading terms as x — )• 1 in the quark contribution Iqq, shown in Eq. (6.1 ). Hence, the quark 
non-threshold contributions shift the result from the dotted to the dashed line, and the 
(non-threshold) gluon contributions shift it from the dashed to the solid line. At the LHC, 
there is a noticable cancelation between the non-threshold quark and gluon contributions, 
which is less prominent at the Tevatron. The reason is that the gluon PDF has a much larger 
effect on the antiquark beam function than the quark beam function. The non-threshold 
quark contributions are as large or larger than the threshold quark contributions, so there is 
no indication that the threshold terms are numerically dominant. The only exception is for 
the Tevatron, where because of the lower center-of-mass energy the threshold corrections 
dominate for large Q ~ 1 TeV, corresponding to x ~ 0.5. In this region the absolute cross 



section in Fig. 8-1 however is very small. 



The left panel of Fig. 8-3 shows the beam thrust cross section da/dQdYdrB at the LHC 
for small tb and at fixed 1" = and Q = 100 GeV. The bands show the perturbative 
uncertainties and results are shown at NNLL (dark orange) , NLL (medium blue) , and fixed 



148 




> 
o 

T3 

-0 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




: 1 
^^^^^ — 


7TeV 


^ ' 


Q = 

: Y = 


100 GeV ^ 




















MNNLL J 






(no ff) E 






(^^1) = 






HNLL 






NLO ^ 


, , 1 







-1.5 



,,cut 



-0.5 



Figure 8-3: The beam thrust cross section for fixed Q = 100 GeV and y = at the LHC 
with Ecra = 7TeV at NLO, NLL, and NNLL. The bands indicate the perturbative scale 
uncertainties as explained in the text. Left: The differential cross section at small . Right: 
The cross section integrated up to tb < exp(— 2y^*) as function of y^*. Also shown are 
the NNLL results in the threshold limit x — )■ 1 (dotted) and without the gluon contribution 
(dashed) . 



order NLO (light green). The NLO result is obtained from Eq. (8.2) by using a common 



fixed scale = /"_b = /^s = with central value fj, = Q and varying ^ = {2Q, Q/2}. This 



reproduces the leading singular terms for — )• shown in Eq. (7.91 ). The full NLO result, 



which is not yet known, contains additional nonsingular terms that are relevant for ~ 1, 



but for ^ 1 constitute power corrections and are neglected in Eq. (8.2). The summation 
has a large effect at small tb and effectively regularizes the IR singularity in the fixed-order 
result. The curves are not plotted for tb < 0.01, because at this point the soft scale drops 
below 1 GeV. Near this cutoff, the soft function becomes nonperturbative and so we expect 
large corrections to our purely perturbative results. 



In the right panel of Fig. 8-3 we plot the corresponding cross section integrated up to 
Tb < exp(— 2y^*) as function of y^*. For nonzero (small t^^^) the summation is again 
important, while for y^* — )• {t^^^ — )• 1) no summation is necessary. The full NNLL result 
provides a smooth interpolation between these two regimes, approaching the singular NLO 
result for yg"* — )- and the NLL result for yg^* > 1-5. 

In conclusion, beam thrust in Drell-Yan provides an experimentally and theoretically 
clean measure of initial-state radiation in qq — )• e^e~ , similar to how usual thrust measures 
final-state radiation in e~^e~ — t- qq. The experimental measurement of beam thrust will con- 
tribute very valuable information to our understanding of initial-state radiation at hadron 
colliders. It should be a priority in the early LHC runs when backgrounds from pile-up 
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effects are still small. We have presented theoretical results for the Drell-Yan beam thrust 
cross section at NNLL, showing that the effects of initial-state radiation in quark initiated 
processes are important but also under good theoretical control. We will see an example of 
a gluon initiated process in the next chapter. 
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Chapter 9 

Beam Thrust Cross Section for 
Higgs Production 

In this chapter we apply the factorization theorem of Ch. [7] to Higgs production through 



gluon fusion. In Sec. 9.1 we discuss the experimental relevance of a central jet veto for 
Higgs searches in the H — t- W~^W~ — t- and — t- 77 channel and present the 

corresponding factorization theorems. Plots for the cross section are shown and discussed 



in Sec. 9.2 These results will appear in Ref. 



9.1 Introduction 

The discovery of the Higgs boson is a major goal of the Large Hadron Collider (LHC) and 
current analyses at the Tevatron. The H — )• W~^W~ — )• H.'^vH.^u channel has strong discovery 
potential and plays a very important role for early searches that are statistics limited. It 
is the dominant channel in the current Tevatron exclusion limit [5J. The presence of the 
final-state neutrinos does not allow the reconstruction of the Higgs invariant mass, and 
hence for this channel sideband methods cannot be used to determine the backgrounds 
directly from data. At the LHC, tt — t- W~^W~bb events constitute by far the largest 
background, dominating the signal by a factor of more than 30 : 1 after applying lepton 
selection criteria [31 H] (see Ref. [IH] for a sensitivity update at 7 TeV) . Requiring a minimum 
missing energy is not effective against this background since it also contains two neutrinos. 
To eliminate the huge background from top-quark decays one imposes a stringent jet veto to 
define a 0-jet sample for the search. For example, ATLAS rejects all events that contain any 
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hard jet with px > 20GeV and pseudorapidity |??| < 4.8, which reduces the background 
by a factor of more than 300 [3]. After the jet veto, the main background is from the direct 
production pp — )• W~^W~ , which at this point stiU dominates the signal by a factor of 4 : 1, 
and the final background discrimination is achieved by a fit to several kinematic variables. 

Theoretically, the inclusive Higgs production cross section has been studied extensively 
in the literature and is known to NNLO [H ESI EZl [ml El ESal USSl [l63] (for reviews 
and additional references see e.g. Refs. [511 EE])- However, Higgs production in a 0-jet 
sample differs substantially from inclusive Higgs production. In particular, the jet veto 
induces large double logarithms that are not present in the inclusive cross section, and also 
induces dependence on the choice of jet algorithm used to define the veto. The current 
method employed by experiments to study the effect of the jet veto and the accompanying 
large logarithms is to use parton-shower Monte Carlos, such as MC@NLO |981 |99] and 
Pythia [160]. This allows one to take into account the dependence of the 0-jet sample 
on the choice of jet algorithm, but for the large logarithms it limits the accuracy to the 
leading- logarithmic summation inherent to parton showers. 

Theoretical studies of the jet veto at fixed NNLO |171 [65l and including additional 
kinematic selection cuts |11H131 1106j are available (see also Ref . |49j ) . The comparison |1H 
dSllMj of the fixed NLO results with MC@NLO, which combines NLO with a LL parton- 
shower resummation, shows indeed discrepancies that indicate the importance of resumming 
the phase-space logarithms caused by the jet veto. The disagreements are reduced but 
not eliminated when the fixed NNLO prediction is compared with those from MC@NLO, 
Herwig |801 ISlj. and Pythia (reweighted to the same inclusive NNLO cross section). There 
are also notable differences between the Herwig and Pythia parton-level results, which is 
another indication that subleading phase-space logarithms are important. 

Theoretically, one can also study the Higgs production as a function of the Higgs trans- 
verse momentum, pff , both in fixed-order perturbation theory for large p^ |18 | l86 iri051ll54| 
and with a resummation of logarithms of pff at small p^ \25\ [50l [551 [561 [^ 11311 1142j . The 
latter is motivated by the fact that the jet veto automatically forces p^ to be small, see e.g. 
Refs. im [131 183j . However, a restriction to small ptf by itself does not provide a jet veto 
as it still allows for configurations with back-to-back hard final-state jets, so a study of the 
small-p^ spectrum can only provide a qualitative template for the effect of the jet veto. 

Here, we use beam thrust to impose the jet veto in pp — )• HX and pp — >• HX. This 
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allows us to directly predict a 0-jet Higgs production cross section, including higher-order 
logarithmic resummation and without relying on parton showers or hadronization models 
from Monte Carlo. The corresponding cross section for the pp — )• WW + Oj background can 
be computed in a similar manner, which is left for future work. 

While H —7- WW provides the most obvious motivation for studying the effect of jet 
vetos, we will also consider the case oi H ^ 77. Here, sideband methods are available 
to experimentally control the overwhelming QCD background (which is up to six orders 
of magnitude larger than the signal at the LHC), but it may still be interesting to study 
sidebands in combination with an explicit jet veto using an inclusive variable such as beam 
thrust. By suppressing the backgrounds from pp — )• jj and pp — )• J7 this could reduce the 
sensitivity of the measurement to the details of the photon isolation cuts. For our purposes, 
we use this channel to study the effect of the jet veto for the case where the Higgs invariant 
mass distribution and rapidity can be measured. With the ability to measure these variables 
the optimal beam thrust variable for — t- 77 differs slightly from that for H — t- WW, as 
we discuss below. 

It should be emphasized that Tevatron Higgs searches, which exclude the range mn = 
162 — 166 GeV at 95% confidence level [SHZ], already analyse their data using a jet algorithm 
and Monte Carlo to implement a jet veto on jets withp^^ > 15 GeV and \r]j\ < 2.4 — 2.5. The 
motivation is again to eliminate backgrounds with additional jets. The CDF analysis [6] 
explicitly uses separate 0-jet, 1-jet, and > 2-jet signal samples, where the latter are included 
to maximize the discovery reach due to the presence of ZH, WH, and vector boson fusion 
production channels. The DO search previously used a (0 + l)-jet sample, where in the 
latest update [7] the number of jets is used as input to a neural network. Even though 
the different jet samples are combined in the final exclusion limit, they are likely to have 
different efficiencies and systematics, with most of the sensitivity coming from the 0-jet 
sample. Hence, it is important to have a theoretical handle on the 0-jet production cross 
section also at the Tevatron. Furthermore, it will be important to use the available Tevatron 
data to test our method in preparation for the LHC searches. 

At the LHC, the Higgs plus 2-jet channel with H — )• WW is also considered. In this 
thesis we restrict ourselves to studying an inclusive 0-jet veto using beam thrust. An 
analogous inclusive variable, A^-jettiness, that works in the presence of A^-signal jets will be 



discussed in Ch. [I0]and can be used to study the exclusive H + Ij and H + 2j cross section. 
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To implement a jet veto in H ^ WW where missing energy plays an important role, 
the appropriate version of beam thrust is defined in the hadronic center-of-mass frame 



^cm 

Tb 



'Tcra 
I B 

VflH 



-T-cm 
'B 



fcT e 



-1%! 



(9.1) 



For the case studied here where the mass of the Higgs is unknown, the dimension-one 
variable T§™ is more convenient than the dimensionless r^™. The sum over k runs over 
all particles in the final state, with pkT and r]k the measured transverse momentum and 
rapidity of each particle, excluding the signal leptons from the W decays. For simplicity we 
take all particles to be massless. 

The production cross section from gluon fusion, gg — )• H, for Tb ^ fnn (i.e. <^ 1) 



is given by the factorization theorem [see Eq. (2.22)] 
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where 



rriH -Y 



V2Gf m]j 



(9.3) 



727r(iV2 - l)i?2^ ' 

and the limits on the Y integration are ln(m///£'cm) ^Y < — In^rriH / Ecm)- 

For 77—7-77 the rapidity Y of the Higgs is measurable and provides an estimate of the 
boost of the partonic hard collision relative to the hadronic center-of-mass frame. In this 
case it makes sense to account for the boost of the hard collision when imposing the jet 
veto. Therefore we use the usual beam thrust 
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In this case, the factorization theorem for the production cross section is 



da 



dYdTB 
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(9.5) 



We will focus on the Higgs production cross section. The leptonic decay of the Higgs, 
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which is of course important in practical applications, does not alter the structure of the 
factorization theorem. It can be included straightforwardly as was done in Ref. |166j for 
the simpler case of pp — )• Z/7 — )• i'^l^ . In the narrow width approximation one can simply 
multiply the production cross section by the relevant differential branching fraction. 



9.2 Plots 

We will now show results for the 0-jet Higgs production cross section from gluon fusion. 



We restrict ourselves to the cross section for Tb in Eq. (9.5), because there is no distinction 
between Tb and T§™ at NLL and the difference at NNLL is very small for the cross section 
integrated over Y (the distinction is important for large Y but this only gives a small con- 
tribution to the total cross section). We evaluate the hard, beam and soft functions at their 
natural scale, where they contain no large logarithms, and run them to a common scale 
using their RGE, which effectively sums the large logarithms. All the necessary perturbative 
results are collected in App. |Ej and we evaluate the convolutions of plus distributions using 
the identities from App. B of Ref. [I35]. We use the MSTW2008 [EE] parton distributions 
at NLO for as{mz) = 0.12 and with two-loop, five-flavor running for as(/u). 

The perturbative uncertainty bands in the plots are estimated by taking the envelope 
of the following three separate scale variations 

a) HH = rruH , l^B = ry^TsmH , IJ-s = rTs , 

b) ^H = mH, fiB = r-^'-^^l'^^^I^^TBmu, P^s = Tb. 

c) iXH = mH. [iB = ^/fmm, ixs = r-(i°'^s/™«)/^rB , (9.6) 

with r = {1/2, 2}, and r = 1 corresponds to the central value. This procedure for estimating 



the uncertainties is identical to Drell-Yan and is discussed further around Eq. (8.3). For 
the cross section integrated up to Tb < T§"*, the scales are chosen by replacing Tb with 



T§"* in Eq. (9.6) 



In Fig. |9-1| we show the Higgs production cross section from gluon fusion at the LHC 
with Ecra = 7TeV (left panel) and the Tevatron (right panel), where we imposed the jet 
veto by requiring Tb < 20 GeV. The large shift from NLL to NNLL is due to the large NLO 
corrections rather than NNLL resummation. This is clear from the NLL' result (dashed), 
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Figure 9-1: The Higgs production cross section from gluon fusion at the LHC with Ecm = 
7 TeV (left panel) and the Tevatron (right panel) with a cut on beam thrust Tb < 20 GeV. 
Shown are the LL, NLL and NNLL results with perturbative uncertainties as well as the 
NLL' results, which combine NLO fixed order corrections with NLL resummation. 



which combines the fixed order NLO corrections with NLL resummation and accounts for 
almost the entire shift. These large perturbative corrections are well-known from fixed order 
calculations. Although it is not very visible, the relative uncertainties decrease slightly when 
going from NLL to NNLL. The larger relative uncertainties compared to Drell-Yan are not 
surprising since Ca/Cf = 9/4. 



Figure 9-2 shows the beam thrust dependence of the production cross section for a Higgs 
with mass niH = 165 GeV at the LHC with E^ra = 7 TeV (left column) and the Tevatron 
(right column). Li these plots we include the singular NLO result, which corresponds to 



evaluating the factorization theorem in Eq. (9.5) at = Q without any running, i.e. without 



any resummation. The top row of Fig. 9-2 shows the cross section as function of Tb- We 
observe that the resummation of logarithms regularizes the IR singularity for Tb — )• in the 
fixed-order NLO result, as was the case for Drell-Yan. Compared to the Drell-Yan result 



in the left panel of Fig. 8-3, the peak occurs at a larger value of Tb and the tail is larger 
as well. This implies that there is more initial state radiation, gluons radiate more than 
quarks. We did not plot the curves for Tb < 1 GeV, because fJ-s = Tb and we therefore 
expect large corrections to our perturbative soft function in that region. 



The bottom row in Fig. 9-2 shows the corresponding cross sections integrated up to 
Tb < 7^"* as a function of T§"*. For T^*^* — s- the resummation is most important, while 
for large Tg"* the perturbative corrections are most important. We see that the NNLL result 
interpolates between these regions, approaching the sing ular NLO result as T§"* becomes 
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Figure 9-2: The Higgs production cross section though gluon fusion for mn = 165 GeV at 
the LHC with Ecm = 7TeV (left column) and the Tevatron (right column). Top row: The 
cross section differential in Tb- Bottom row: The cross section integrated up to Tb < Tg^* as 
function of T§"*. Shown are the LL, NLL and NNLL results with perturbative uncertainties 
and the NLL' and singular NLO result. 

large and the NLL and LL result for Tn'^^ — s- 0. Compared to the results for Drell-Yan 



shown in the right panel of Fig. 8-3 the size of the corrections as well as the effect of 
resummation is much larger. 
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Chapter 10 

AT-Jettiness: An Inclusive Event 
Shape to Veto Jets 

In this chapter we will extend our work to events with N signal jets and consider the 
veto on additional unwanted jets. We introduce a global event shape 'W-jettiness" r^r to 
impose this veto, which for = reduces to the beam thrust. We discuss the experimental 
relevance of jet vetos and the benefits of using r^v in theory calculations. In particular, r^r 
allows us to sum large logarithms due to phase space restrictions and leads to a factorization 
formula with inclusive jet and beam functions. The work presented in this chapter was first 
reported in Ref. |167j . 

10.1 Introduction 

At the LHC or Tevatron, hard interactions involving Higgs or new-physics particles are 
identified by looking for signals with a characteristic number of energetic jets, leptons, or 
photons The backgrounds come from Standard Model processes producing the same 

signature of hard objects possibly with additional jets. An example are top quarks decaying 
into W plus 6-jet, which is a major background for H — t- WW [5]. When reconstructing 
masses and decay chains of new-physics particles additional jets can cause large combinato- 
rial backgrounds. Standard Model processes can also fake a signal when a jet is misidentified 
as lepton or photon, a typical example being H — t- 77. 

Thus, a veto on additional undesired jets is an effective and sometimes necessary method 
to clean up the events and discriminate signal and the various backgrounds. More generally, 
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one would like to measure an "exclusive" A^-jet cross section, pp — )■ XL{Nj), to produce 
N signal jets j where the remaining X contains no hard (central) jets. Here, > and L 
denotes the hard leptons or photons required as part of the signal. 

We introduce an inclusive event shape 'W-jettiness" , denoted r^v and defined below in 



Eq. (10.1). For an event with at least N energetic jets, tn provides an inclusive measure of 
how A^-jet-like the event looks. In the limit tn ^ the event contains exactly N infinitely 
narrow jets. For r^r ~ 1 the event has hard radiation between the signal jets. Requiring 
ttv ^ 1 constrains the radiation outside the signal and beam jets, providing an inclusive 
way to veto additional central jets. It yields an inclusive definition of an exclusive A^-jet 
cross section with a smooth transition between the case of no jet veto, rjv ~ 1, and the 
extremely exclusive case rjy — )• 0. 

Vetoing additional jets imposes a phase-space restriction on the underlying inclusive 
A^-jet cross section to produce A^ or more jets with the same L. Irrespectively of its precise 
definition, the jet veto introduces a jet resolution scale fij that characterizes this restriction, 
i.e. the distinction between A^ and A^ + 1 jets. Hence, the exclusive A^-jet cross section 
contains phase-space logarithms a" In"^ {fi'^j / fijj) , where m < 2n and ij-h is the scale of 
the hard interaction. For r^r, fj.'j/fi'jj ~ r^v ^ 1. Generically there is always a hierarchy 
/^J ^ A*// 5 which becomes larger the stronger the restrictions are. These large logarithms 
must be summed to obtain reliable predictions. 

Jet vetoes are typically implemented by using a jet algorithm to find all jets in the 
event and vetoing events with too many energetic jets. Jet algorithms are good tools to 
identify the signal jets. However, they are not necessarily well-suited to veto unwanted 
jets, because the corresponding phase-space restrictions are complicated and depend in 
detail on the algorithm. This makes it difficult to incorporate the jet veto into explicit 
theoretical calculations and in particular inhibits a systematic summation of the resulting 
large logarithms. In this case, usually the only way to predict the corresponding exclusive 
A^-jet cross section is to rely on parton shower Monte Carlos to sum the leading logarithmic 
(LL) series. 

In contrast, vetoing jets by cutting on an inclusive variable like tn has several advan- 
tages. First, we can go beyond LL order, because the logarithms from the phase-space 
restriction, a^lrJ^T^, are simple enough to allow their systematic summation to higher 
orders. Moreover, the theory predictions with factorization can be directly compared with 
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(c) pp — leptons plus jets. 

Figure 10-1: Different situations for the application of A^-jettiness. 



experiment without having to utilize Monte Carlos for parton showering or hadronization. 
Experimentally, tn reduces the dependence on jet algorithms and might help improve the 
background rejection. 



10.2 Definition 

A'^-jettiness is defined as 



TN = ^ 2Zmin{Q'a -Pk, Qb-Pk, Qi ■ Pk, ■ ■ ■ , QN ■ Pk} ■ (10.1) 
^ k 

As we discuss below, this definition of r^v yields a factorization formula with inclusive jet 
and beam functions and allows the summation of logarithms to next-to-next-to-leading 



logarithmic (NNLL) order. The sum over k in Eq. (10.1) runs over the momenta pk of all 
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measured (pseudo-)particles in the final state excluding the signal leptons or photons in L. 
(Any other leptons or photons, e.g. from hadronic decays, are included in the sum.) For 
simplicity we take all pk to be massless. The qa, qt, and qi, q^ are a fixed set of massless 
reference momenta for the two beams and the A*" signal jets, 

<6 = 5^a,b^cm<,6, < = (1,0, 0,1), < = (1,0,0,-1), 

q^j = Ej{l,hj) , J = {l,...,iV}. (10.2) 

The Ej and hj correspond to the energies and directions of the signal jets. Their 
choice is discussed below. The beam reference momenta qa and q^ are the large momentum 
components of the colliding partons along the beam axis (taken to be the z axis). They are 
defined by 

Xa-Ecm = nb • (^1 H VqN + g), (10.3) 

and analogously for Xh with a b. Here, q is the total momentum of the non-hadronic 



signal L. In Eq. (10.1), Q = XaX},E^^ is the hard interaction scale, and the distance of a 
particle with momentum pk from the jets or beams is measured by qm ' Pk- If contains 
missing energy, so q and Xa^b are not known, one can use a modified distance measure as 



we discuss below Eq. (10.11). 



The minimum for each k in Eq. (10.1) associates the particle with the closest beam 
or jet, appropriately dividing the hadronic initial-state radiation (ISR) and final-state ra- 
diation (FSR). Soft particles and energetic particles near any jet or beam only give small 
contributions to the sum. For 2 — )• A^ scattering of massless partons, tn = 0. Energetic 
particles far away from all jets and beams give large contributions. Hence, for ttv *C 1 the 
final state has A^ jets, two forward beam jets, and only soft radiation between them. In this 
limit Xa^b are the momentum fractions of the annihilated partons, and Y = \n{xa/xi,)/2 is 
the boost of the partonic center-of-mass frame. 



10.2.1 AT = 2 for e+e" Jets. 



In e~^e collisions there is no hadronic ISR, so we drop the qa,b ' Pk entries in Eq. (10.1). 
Now is the total invariant mass of the leptons and Y = 0. In the two-jet limit, the jet 
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directions are close to the thrust axis t, so we can choose 



q^=^Q{l,-i) 



(10.4) 



as reference momenta, and Eq. (10.1) becomes 



rl*^ = ^ ^ Ek min{ 1 - cos 9k, 1 + cos ^a;} 



(10.5) 



where 9k is the angle between pk and t. The minimum divides all particles into the two 



hemispheres perpendicular to t as shown in Fig. lO-l(a) For rl"^ <^ 1, the total invariant 



mass in each hemisphere is much smaller than Q, so the final state contains two narrow 
jets. In this limit, = 1 — T, where T is thrust, and a factorization theorem exists for 
dfj/drl'^, which can be used to sum logarithms of |431 [Ml US]. 



10.2.2 AT = for Drell-Yan 



Next, consider isolated Drell-Yan, pp — t- Xi~^i~ with no hard central jets, shown in Fig. 10- 
l(b)[ We now have ISR from the incoming partons, but no FSR from jets. From Eq. (10.3) 
we have 

XaEcm = e^^J q^ + q^, XbE^m = e~^J q^ + q^ , (10.6) 



where q^ and qx are the dilepton invariant mass and transverse momentum, and Y equals 
the dilepton rapidity. Now, Q"^ = q^ + q^ and Eq. (10.1) becomes 



To 



Q 



PkT\ mmie 



-Y+Vk 



}• 



(10.7) 



Thus, the qa and qi, dependence in Eq. ( 10.1 ) explicitly accounts for the boost of the partonic 



center-of-mass frame. For Y = 0, the minimum in Eq. (10.7) divides all particles into two 



hemispheres perpendicular to the beam axis (analogous to the thrust axis in e^e ). For 
Y ^ 0, the hemispheres are boosted with their boundary now at Y, and the beam jet in the 



direction of the boost is narrower than the other, as depicted in Fig. 10-1 (b) Contributions 
to tq from large rapidities are exponentially suppressed by |pfcT|e~'''*' ~ 2£'fce~^l'"^l, so 
particles beyond the detector's rapidity reach give negligible contributions. 

Beam thrust is given by tb = \/l + q^/q'^ tq. It is obtained by choosing XafiEcm = 
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g2g±y gj^gg q2 g^j^^ Y are measured rather than the longitudinal components • g in 



Eq. ( 10.6 ). For tq ^ 1 the hadronic final state can only contain soft radiation plus energetic 
radiation in the forward directions. In this limit the leptons become back-to-back in the 
transverse plane, qx ^ Q, so tb = tq. In Ch. [7| we derived a factorization theorem for 
dfj/dr^ at small r^, which we used to sum logarithms of tb to NNLL. 



10.2.3 General Case 



For pp — )• XL(Nj) we have both ISR and FSR. We select candidate signal events by 
measuring L and running a jet algorithm to find the N signal jets and their momenta pj. 
The conditions on the jets and L that define the signal are encoded in the cross section 
by a measurement function F/v({pj},L). Generically, -F/v will enforce that there are at 
least energetic jets that are sufficiently separated from each other and the beams. We 
now use the measured jet energies and directions to define the jet reference momenta qj in 



Eq. (10.2), 



Ej=p°j, nj=Pj/\pj\, 



(10.^ 



while qa and qh are given by Eqs. (10.2) and (7.30) 



Taking the minimum in Eq. (10.1) combines the previous cases in Eqs. (10.5) and (10.7). 



It divides all particles into jet and beam regions that are unique for a given set of reference 



momenta and whose union covers all of phase space, as illustrated in Fig. 10-1 (c) The 
boundary between any two neighboring regions is part of a cone and is such that the sum 
of the total invariant masses in the regions is minimized (or in case of a beam region the 
virtuality of the incoming colliding parton). 

For events with small r^v all jet algorithms should agree how energetic radiation is split 
up between the jets and beams, and only differ in their treatment of softer particles. Thus, 
they all give the same hj and Ej up to power corrections, while the split up of the soft 
radiation is determined by tn itself. Hence, the dependence of tat on the jet algorithm is 



formally power suppressed, t^^'^ = t'^'^"^ -\- 0{t'^), as seen in Eq. (10.14) below. 

To measure tat, we still rely on having a suitable jet algorithm to find the N signal jets 
but not more so than if we were not measuring rjy. Imagine the jet size in the algorithm 
is chosen too small such that the algorithm divides what should be a single signal jet into 
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several narrow jets ^. In this case, the jet algorithm yields a poorly reconstructed signal 
irrespective of measuring tm- 

Since the jet veto is now provided by r^v, this situation can be avoided because we do 
not have to rely on the jet algorithm to identify additional jets and so can use an algorithm 
that can be forced to always yield at most N jets. This is in fact the most natural thing to 
do when one is looking for N jets. Therefore, using rjv as jet veto could also help improve 
the signal reconstruction. 

10.2.4 Generalizations 

We can generalize tn to 

Tn = '^'aim{da{pk),db{pk),di{pk),. . .,dN{pk)} , (10.9) 



where dm{pk) can be any infrared-safe distance measure. In Eq. ( 10.1 ), dm{pk) = '^Qm.-Pk/Q^ 
with 

"^Qa-Pk = Ip/ctI Qe^"''* , 

2gj • Pk = \PkT\ WjtI (2 cosh Arjjk - 2 cos Acpjk) . (10.10) 

Here, Aijjk and A(j)jk are the rapidity and azimuthal distances between qj and pk- If these 
are small, the factor in brackets reduces to the familiar = (Ar/)^ + (Ac/))'^. 

Raising a given measure to a power, dmiPk)°^"^ changes the relative weight between the 
region's center and periphery in the sum over k, while the powers Om can be used to change 
the division between beams and jets. Different measures that are boost-invariant along the 



beam axis can be obtained by raising each factor in Eq. (10.10) to a different power. A 
geometric measure, which is independent of I^jtIi is 

daiPk) = ^ e^-^" , dbipk) = ^ 6-^+"^ , 

dj{pk) = ^ (2coshA?7jfc - 2cosA0jfc) . (10.11) 



'^This can be tested by comparing the total energy in each region defined by rjv with the energy from 
the jet algorithm. If these are very different, but at the same time rjv is small, then there are additional 
energetic particles near the signal jets that the algorithm should have included. 
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It evenly divides the area rather than invariant mass between neighboring regions, such that 
more energetic jets also get more invariant mass. 



If L contains missing energy, Xa^h in Eq. (10.3) and thus Q and Y are not known. For Q, 
one can use any hard scale, like the \qjT\ of the hardest jet, since it only serves as an overall 



normalization. In the beam measures da^b{Pk) in Eq. ( 10.11 ) we can simply set y = 0, which 
defines them in the hadronic center-of-mass frame and effectively averages over all possible 
boosts Y . 

A^-jettiness does not split events into A^, + 1, + 2, etc. jets like a traditional jet 
algorithm. But we can consider using to define an "exclusive A^-jet algorithm" as follows: 
First, we use a geometric measure and find the directions hj and boost Y that minimize r^v. 
This is analogous to determining the thrust axis in e^e~ — )• jets. It might actually allow 
one to get an estimate of Y even in the case of missing energy by exploiting the asymmetry 
in the beam jets. Second, we determine the jet energies by summing over the particles in 
each jet region. (To reduce the sensitivity to the underlying event and pile-up, one can 
weigh the sum over energies by the distance from hj.) 



10.3 Factorization formula 



For simplicity, we now use r^r again as defined in Eq. (10.1). For tat <C 1, QCD ISR and 
FSR can be described in soft-collinear effective theory [29H3H [351 EZ] at leading power by 
A^ + 2 independent sectors for coUinear particles close to each with m = {a, 6, J} and a 
separate sector for soft particles. By power counting, J-collinear particles are closest to qj, 
so the sum over k for the J-collinear sector is 

vtmi[2qm ■ Vk] = ^ 2qj-pk = sj, (10.12) 

fcScoUj fcScoUj 

where (up to power corrections) sj is the total invariant mass in the J-collinear sector. 
Similarly, the sum over the beam collinear sectors yields the total (transverse) virtuality of 
the colliding partons, ta and t^. Therefore, 

tnQ'^ = ta + tb + '^sj + ^ mm{2qm-pk}- (10.13) 

J fcGsoft 
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The sum in the last term is now restricted to the soft sector. Combining Eq. (10.13) with 



the analyses in Ref. [32] and Ch. [t] yields the factorization formula for A^-jettiness ^ 

^= jdxadxb jd\d^L{q) jd^N{{qj])FN{{qm],L) (2TT)H^{qa + qb-Y,qj-q 
ti Hij^f,{{qm], L, /i) jdtaBi{ta,Xa,n) j dth Bj{tb, Xb, fJ-)Y{ j (isj Jf,j{sj, 



X 



X Sr\r^ - ^"^''Q2^''^ {g>n},/iJ . (10.14) 

Hij^K{{qm},L) contains the underlying hard process i{qa)j{qb) L{q)Ki{qi) ■ ■ ■ KN{qN), 
where i,j,Kj denote the parton types, and the sum over ij,K. is over all relevant partonic 
channels. It is a matrix in color space given by the IR-finite parts (in pure dim. reg.) of 
the squared partonic matrix elements in each channel. The A^-body phase space for the 
massless momenta qj is denoted d^ni{qj}), and that for L by d$L(g). 

The distributions in sj and ta^b are described by inclusive jet and beam functions, 



Jkj{sj) and Bij{ta^b, Xa,b)- The last term in Eq. (10.13) is the contribution to tn from soft 
particles in the "underlying event". It is described by the soft function 5'^~^'^(t^^*, {gm}), 
which depends on the jet's angles hi ■ hm and energy fractions Ei/Em- Like H, it is a color 



matrix, and the trace in Eq. (10.14) is over tr(HS). 

In Eq. (10.14), all functions are evaluated at the same renormalization scale fi. Large 
logarithms of r^v in dcr/drjv are summed by first computing H{^h), J{fJ'j), B{^b)i S{fis) 
at the scales fin ~ Q, ^J'J ~ fJ-B ~ V^nvQ, ^5 — tnQ, where the functions contain no large 
logarithms, and then evolving them to the scale /i. This evolution is known analytically |94| 
I135t I151| and the required anomalous dimensions are already known to NNLL |92[ I112| 
11271 11301 11461 1147j . because we have inclusive jet and beam functions. NNLL also requires 
the 0(0^) corrections for each function, which are known for J and B. The 0{as) hard 
function is determined by the one-loop QCD matrix elements. For tj^ » ^qcd/Q, S{ijs) 
can be computed perturbatively and will be given in a future publication. 



^Here, Fn enforces distinct coUinear sectors with 1 — ■ fim ^ tjv and Em/Q ^ tat. We assume Fn only 
depends on the large components qj of the jet momenta, pj — qj [1 + C'(rjv)] , and that L only couples to the 
QCD subprocess via a hard interaction. We also assume that Glauber gluons do not spoil this factorization. 
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Chapter 11 



Conclusions and Outlook 



Typical Higgs and new physics searches at the LHC or Tcvatron require a certain number of 
hard jets, leptons and photons. Jet algorithms are good for identifying jets, but, as we have 
seen, a veto on additional unwanted jets is most easily incorporated in theory calculations 
by using an event shape. This approach allowed us to sum the corresponding large phase 
space logarithms and make predictions without relying on Monte Carlo programs for parton 
showers or hadronization models. 

A lot of our work has focussed on the 0-jet case, for which we implemented the jet veto 
by restricting the beam thrust variable -C 1. A central jet veto is needed to reduce the 
large background in H ^ WW — )■ from top quarks decaying into a W plus 6-jet. We 
derived a factorization theorem that allowed us to sum the large phase space logarithms 
a" In"* tb- The central jet veto also probes the initial-state radiation, which requires beam 
function to describe the initial state. Beam functions describe extracting a parton out of the 
proton as well as the formation of an initial-state jet. They can be related to PDFs through 
a perturbative matching calculation, which we carried out at one loop. Our calculations 
allow us to study the effect of initial-state radiation in perturbation theory and the next 
step is to compare our analytical results with Monte Carlo programs. Another interesting 
extension of our work would be to include the measurement of the recoil of the initial-state 
jets, which requires px dependent beam functions. 

We concluded this thesis by considering the extension to N signal jets, for which we 
introduced the event shape A/"-jettiness to veto unwanted additional jets. Many of the 
features of the 0-jet case carry over, including the resummation of large logarithms and the 
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presence of beam functions. In the 0-jet case there was no jet algorithm dependence, in 
this case a jet algorithm is used to identify the jets but the dependence on the choice of jet 
algorithm is power suppressed. We still need to show the cancellation of Glauber gluons. 
Once the one-loop soft function has been calculated we can use our factorization formula 
to determine cross sections for processes with final state jets at NNLL. This will allow us 
to study the jet energies and directions, which are of interest to new physics searches. To 
study other properties of the jets, such as their invariant masses, requires an extension of 
our factorization formula. 
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Appendix A 



Notation Cheat Sheet 



A.l Lightcone Coordinates 

We use lightcone vectors and with = = and n-n = 2 to decompose four- vectors 

p^ = p — + p^ — + p'^ , p^=n-p, p~=n-p. (A.l) 

We have two coUinear directions ria = (1, 0, 0, 1) and = (1,0, 0,-1) corresponding to the 
directions of the incoming protons, and take = and = n^. Lightcone coordinates 
of a momentum with subscript a or 6 will be with respect to or Ub, e.g. = Ua-Ba, 
B^ = ni)-Bb. We will sometimes write px = P±- 

For an energetic particle collinear to the direction we separate its momentum p^ into 
a (discrete) large p^ and (continuous) small momentum pr in SCET, p^ = p'^ + Pr- The 
large momentum p^ = (0,p~,p^) has components p~ = n-p and ~ ^P~ , where A is the 
power counting parameter. The continuous small momentum has scaling pr = }?p~ . The 
large momentum appears as a label on the field and is picked out by the label operator 
v. The small momentum is picked out by p and is often written in the conjugate position 
space coordinate y. We absorb the residual p^ and p^j_ into the label p, to give continuous 
labels. This leaves us with just one residual coordinate y~ . More details on SCET can be 



found in Sees. 1.4, 3.1 and 7.1 
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A. 2 Symbols 

Below we list the most common symbols and their definition. Whenever there is a subscript 
a there will be a corresponding symbol with a subscript 6, corresponding to the other pro- 
ton, parton or hemisphere. The subscripts i and j run over flavors (quark or gluon). 





Symbol Definition 


General 


£ polarization vector 

L non-hadronic final state for Drell-Yan) 

ria lightcone vectors in the direction of the incoming protons 

Pn n-coUinear proton 

Unip) n-coUinear spinor 

X total hadronic final state 


Momenta 


ba collinear radiation from the incoming parton a 

Ba hadronic momentum in hemisphere a, Ba = ba + rH + ka 

ka soft momentum in the hemisphere a 

Pa momentum of proton a, Pa = -E'cm'T'a/2 

Pa momenta or mconniig parion a 

p'^2 momenta of the leptons for Drell-Yan 

Ta momentum of the proton remnant for proton a 

qf^ total non-hadronic momentum (g'^ = Pi + P2 Drell-Yan) 


Kinematics 


oja^b large momentum component of incoming parton, coa ~— Pa 
Q inv. mass of non-hadronic particles, hard scale, Q = 
t transverse virtuality of incoming parton, ta = oJaba 

r = qVEl^ 
tb beam thrust 

Xa momentum fraction of parton, Xa = Ua/Pa 
^a momentum fraction of parton in the PDF 

y'^^ rapidity "cut" corresponding to beam thrust, tb < exp(— 2y^*) 

Y total rapidity of non-hadronic particles 

z parton momentum fraction of in partonic calculation 


Scales 


fiB beam scale, /ns ~ ^/tbQ for beam thrust 
fiH hard scale, f^H Q 

fi\ low scale at which the PDFs are defined 
jis soft scale, jJ-s ^ tbQ for beam thrust 
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Symbol 


Definition 


Factor- 


B^ 


beam function 


ization 


Bi 


beam function integrated up to t < tmax 




fi 


parton distribution function 




Hij 


hard function (same for threshold and isolated processes) 






Wilson coefficient for initial state jet, B = X® f 




Ji 


(final state) jet function 




A 


power counting parameter, A ~ /Q ~ tb 




Oi 


beam function operator, Bi = {P~\Oi\P~) 




Qi 


PDF operator, fi = {P~\Qi\P~) 






soft function (differs between threshold and isolated) 






Born cross section 


Renorm- 




shows up in RGE solution for cusp piece 


alization 


7(...,^) 


anomalous dimension 




7("s) 


non-cusp anomalous dimension 




F 

cusp 


cusp anomalous dimension 




Kr 


exponent in RGE solution for cusp piece 






exponent in RGE solution for non-cusp piece 




Pjk{z) 


splitting functions in PDF evolution 




U 


evolution function 




z 


renormalization factor 


SCET 




n-colinear gluon field 




As 


soft gluon field 




I3n± 


gauge invariant collinear gluon field 




Xn 


gauge invariant n-coUinear quark field, Xn{y) = '^n{y)in{y) 




A 


power counting parameter, A ~ B^ jQ ~ tb 






residual momentum operator 






label momentum operator, V^^n,f>iy) = P^in,p{y) 




V 


minus label momentum operator, V = V~ 






n-collinear Wilson line 






n- collinear quark field 




y~ 


residual minus position coordinate 




Yn 


soft Wilson line along the direction (fundamental repr.) 




yn 


soft Wilson line along the direction (adjoint repr.) 
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Appendix B 



Plus Distributions and 
Discontinuities 



The standard plus distribution for some function g{x) can be defined as 

[e{x)g{x)] ^ = lim [9{x - /3) G{x)] with G{x) = dx' g{x') , (B.l) 
satisfying the boundary condition Jq dx [9(x)g{x)]-\. = 0. Two special cases we need are 



Cn{x) = 
a^{x) = 



9{x) In" X 

X 

e{x) 



+ 

lim 



lim 



X n + 1 



■n 



(B.2) 



In addition, we need the identity 



d{x) 



X 



1+e 



6{x)+Co{x)-eCi{x) + 0{e^), 



(B.3) 



the Fourier transform 



and the two limits 



C^{x) = - j^e'-y ln[i(i/-iO)e^-], 



(B.4) 



lim 

/3->-o 



/3^0 X^ 



TT 



6 



= 5{x) . 



(B.5) 
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Away from x = these relations are straightforward, while the behavior at x = is obtained 
by taking the integral of both sides. General relations for the rescaling and convolutions of 
Cn{x) and C^[x) can be found in App. B of Ref. [135J. 
The discontinuity of a function g{x) is defined as 

Disca; g{x) = lim \g{x + i/3) - g{x - i/3)l . (B.6) 

/3->0 

If we are only interested in the discontinuity in some interval in x, we simply multiply 



the right-hand side with the appropriate 9 functions, as in Eq. (3.54). If g{x) is real then 
D\scxg{x) = 2\lmg{x + iO). Two useful identities are 

— Disc^^ = ^Disc^(-a;)"-^ = (B.7) 

2-K ""x^+i n! ^ ^' 27r ^ ' ^ ' n ^ ' ^ ^ 

To derive the last identity, note that 

(-X - iO)"-" = exp[(n - e) ln(-x - iO)] = |x|"-" exp[-i7r(n - e)6'(x)] , (B.8) 

so taking the imaginary part gives Im(— x — iO)"""^ = (— 1)" sin(7re) 9{x) x""*^. To calculate 
the discontinuities of the various graphs in Sec. |5.3| we need the relations, 



- 9{z)^ DisCi>o da A'^-^ = 9{z) ^ 9(1 - z)z'+%l - z)"^ , (B.9) 

- 9iz)^ DisCi>o da (1 - a)A-i- = 9{z) ^ ^(1 - z)z^+^il - zf-^ , 
_ e{z)^ DisCi>o da (1 - a) tA-^- = 9{z) % 9{l - z)z'+%l - z)-^ , 



where A = t{l — a/z). 
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Appendix C 

Renormalization Structure of the 
Beam Function 



In this appendix we derive the general structure of the beam function RGE in Eq. (3.25) 
to ah orders in perturbation theory, which was presented in Ref. |168j . The two essen- 
tial ingredients will be the known all-order renormalization properties of lightlike Wilson 
lines \57\ I124^ \126\ I127j and the factorization theorem for the isolated pp — )• XL cross 
section, where X is the hadronic and L the non-hadronic final state. In Ch. [7] we proved 
that to all orders in perturbation theory and leading order in the power counting this cross 
section factorizes as 

daMyd^+di?+ = ^ ^' I ^^'^ '^''^ ^^'^^ 

b ij 

X q^Bi[uJaiB+ - kt),Xa,^]B,[u:,{B+ - A;+),Xb,^] . (C.l) 

The sum over ij runs over parton species ij = {gg, uu, uu, dd, dd, . . .}. The soft function does 
not depend on the quark flavor, and its superscript only refers to the color representation. 
The three ingredients in Eq. ( |C.l ) are the renormalized hard, beam, and soft functions, 



Hij{q'^ ,Y, fi), Bi{t,x,fi), 5**1^0011 (^o" ' ^fc"' /^)- Their dependence on the renormalization scale 



/i must cancel in Eq. (C.l ), because the cross section must be /j independent. The structure 
of the RGE for the hard and soft functions thus uniquely determines the allowed structure 
of the beam function RGE. 

The hard function is a contraction between the relevant leptonic matrix element squared 
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and the square of the Wilson coefficients of the color-singlet qq and gg local SCET currents 
(see Sec. 7.3.2) 



where a and /3 are spin indices. In each collinear sector, total label momentum and fermion 
number for each quark flavor are conserved. Thus, the currents cannot mix with each other 
and are multiplicatively renormalized. Furthermore, RPI-III invariance implies that the 
RGE for the currents can only depend on q^ = uJa^h- The renormalization of these SCET 
currents also does not depend on their spin structure, so the RGE for the hard function 
must have the same structure as for the currents. Therefore, to all orders in perturbation 
theory we have (with no sum on ij) 

fi^H,,{q'',Y,fi)=j'^{q'',fi)Hij{q'',Y,fx). (C.3) 

Next, the incoming hemisphere soft function, Sl^^^^{k^ , , fi) , is given by the vacuum 
matrix element of incoming soft lightlike Wilson lines along the Ua and Uf, directions. In 
position space, 

slU,{ya,yt,f^) = J dktdk+ e-'^'iy^+'^y^y^ siu,{kt,k+,f,) (c.4) 

has two cusps, one at spacetime position and one at y = y~na/2 + y^nb/2. The renor- 
malization properties of lightlike Wilson lines with cusps |57L 11241 1126^ I127j then imply that 
to all orders in perturbation theory. 



€ iVa , yb > ^) = 2rL,p {as) \- In Uy^^^e^A - In i^-^i^e'- 



where E^^gp is the cusp anomalous dimension for quarks/antiquarks or gluons, and 7^ [as(/u)] 
and E*^gp[as(/x)] depend only indirectly on /i via ««(/«). Dimensional analysis and RPI-III 
invariance imply that the single logarithm multiplying 2E*ygp scales like ln(y~y^/i^). (The 
additional dimensionless factors are chosen for convenience. Any change in them can be 
absorbed into 7^(0^).) The correct overall sign and iO prescription for the logarithms can 
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be deduced from the explicitly known one-loop result see Eq. (7.89) and \94: \ I157j . 



Taking the Fourier transform of the cross section in Eq. (C.l) with respect to and 



-B^ and differentiating the result with respect to ^ yields 



= II 



Va 



d/i 



(C.6) 



The factorization theorem for the cross section neither depends on the choice of L, which 
affects the form of Hij for different ij, nor the type of the colliding hadrons. This implies 
that each term in the sum over ij must vanish separately. (For example, choosing Drell-Yan, 
L = £'^i~ , there is no contribution from ij = gg, so the quark and gluon contributions are 
separately zero. Then, by assigning arbitrary electroweak quark charges, the contribution 
from each quark flavor must vanish separately. Finally, the ij = qq and ij = qq contributions 
for a single quark flavor q must vanish separately by choosing various different incoming 
hadrons.) Therefore, the RGE for the product of the two beam functions is 



d^ 



Bi( ^!^,Xa,fj)Bj(^,Xb,fJ-) =0. 



'\2uJa' 



(C.7) 



which shows that the beam functions in position space renormalize multiplicatively and 
independently of Xa,b- The RGE for each individual beam function can only depend on the 
RPI-III invariant y~ /2uj and obviously cannot depend on the variables of the other beam 
function. Hence, we find that to all orders in perturbation theory 



(C.8) 



which is the result we set out to prove in this Appendix. Using Eq. (C.8) together with 



Eq. (C.7), the anomalous dimensions must satisfy the consistency condition 



(C.9) 



Given the form of 7^ in Eq. (C.5), it follows that the anomalous dimensions are given to 
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all orders by 



(C.IO) 



Taking the Fourier transform using Eq. (B.4), the momentum-space anomalous dimensions 
become 



■^-^ cusp 



t 



1^' 



(C.ll) 



The same all-order structure of the soft anomalous dimension as in Eq. (C.ll) was 
obtained in Ref. [93j for the hemisphere soft function with outgoing Wilson lines in e~^e~ — t- 2 
jets using analogous consistency conditions. In fact, the hard SCET currents here and there 



are the same and in Sec. 3.2 we proved that the anomalous dimensions for the beam and 



jet function are the same, 7^ = jj. Hence, the hemisphere soft functions with incoming 
and outgoing Wilson lines have in fact identical anomalous dimensions to all orders. 
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Appendix D 



Quark Beam Function Matching in 
Pure Dimensional Regularization 



Here we repeat the NLO matching calculation from Ch. |4] using dimensional regularization 
for both the UV and IR, which was reported in Ref. |168j . Since we only change the IR 
regulator, the final results for the matching coefficients Iij{t,z,fi) should not be affected. 

In pure dimensional regularization all the loop diagrams contributing to the bare matrix 
elements of Qq vanish, since by dimensional analysis there is no Lorentz invariant quantity 



they can depend on. Hence, including the counter terms in Eq. (4.17) to subtract the UV 
divergences, the renormalized matrix elements consist of pure IR divergences with opposite 
signs to the UV divergences, 

{gn\Qq{oj,f^)\gnf^ = -^"^^!^9iz)Pgg{z). (D.l) 

e ZTT 

This shows explicitly that the conventional MS definition of the PDFs in QCD, which also 



yields Eq. (D.l), is indeed identical to the SCET definition used in our OPE for the beam 
function. 



Considering the beam function matrix elements, the bare results for Figs. 4-2 (c) and 



4-2(d) now vanish, because their loop integrals are again scaleless. For the remaining dia- 



grams we can reuse the intermediate results from Sec. 4.2 before carrying out the Feynman 



parameter integrals and taking the discontinuity. Setting t' = the denominator in the 
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Feynman parameter integrals in Eq. (4.22) becomes (1 — a)A — aB = t{\ — a/z). In this 
case it easier to carry out the integral after taking the discontinuity. The discontinuity we 
need is 



27r 



-Disc 



i>0 



l--\t 

z 



^i-_sinvre0(t)^^a_^^^a_^^-i-e^ (D.2) 



where we used Eq. (B.7). Since we require z > 0, the first 6 function becomes 9{a — z), and 
so we have 



-e{z)—Bisct>oh{A,B,e) 
Ztt 



vre 



d{l- z)z^+'{l- z) 



(D.3) 



4-2 (a) , using Eqs. ( 


4.23 


) and ( 


D.3 



{qn\9{uj)O^g^''^{t,L0)\qr, 
as{ij)CF 9{z) 



27r z 
as{n)CF 



r(l + e){e^'^^?f{l - ef [--^Disct>o h{A, B, e) 



27r 



9{z)e{l - z){l - z)T{l + e){e<^^y{l - e) 



sin vre 6{t) / z \^ 



vre 



Vi 



where in the last step we used Eq. (B.3) to expand in e. For Fig. 4-2(b) , we start from the 



third line in Eq. (4.29) and using Eqs. (D.3) and (B.3) we get 



as{fJ,)CF Biz) 



r(l + e)(e^^/i' 



7e ,,2\e 



Disct>o/i(A5,e) 



vr 1 — z^''^ '"'L27r 



l+e 



as{p.)CF 



-^{z)T{l + e){e<^^l'r 



TT 



9{z) 



vre (1 - 

- 5(1 -z) + £0(1-2)2 



+ 4^)5(1 -z) + 6(1) \Ci{l - z)z - £0(1 - z)zhiz- ^5{l - z) 

11"^ Xfi^J L 12 



(D.5) 



Adding up Eqs. (D.4) and (D.5), the bare quark matrix element in pure dimensional regu- 



182 



larization becomes 



{qn\eioj)0^^'''^%t,u)\qnf^ 
_ as{ii)CF 



2tt 



2 1 / t 

2-^0 I ~2 



+m 



IT 



£i{l-z){l + z')- —6{l-z) + 9{l-z)(l 



1 + z^ 
1 - z 



Inz 



. (D.6) 



We can now proceed in two ways to obtain the matching coefficient Iqq{t, z, ji). 



First, we can subtract 5{t) times Eq. (D.l) from Eq. (D.6) to obtain the bare matching 



coefficient. This simply removes the [1/ e)5{t)Pqq{z) in the first line of Eq. ( D.6| ). Assuming 
that the IR divergences between the PDF and beam function cancel (and including the 
vanishing zero-bin) the remaining poles in the fo'st line are of UV origin and determine the 



necessary MS counter term, reproducing our previous result for Zg(t,/i) in Eq. (4.40). 



Alternatively, we can use our general result that the beam function has the same renor- 
malization as the jet function. In this case, we subtract the one-loop counter term for 
Qhnre -j^ ( 4.40| ) (which is already known from the jet function's renormalization) from 



Eq. (D.6) to obtain the renormalized quark matrix element, which equals Eq. (D.6) without 



the [•.•](5(1 — z) term in the first line. The remaining 1/e pole must then be of IR origin, so 
we again have an explicit check that the IR divergences in the beam function match those 



of the PDF in Eq. (D.l). Either way, the finite terms in the last two lines of Eq. (D.6) 
determine the renormalized matching coefficient Iqq{t, z, fj,), which agrees with our previous 



result in Eq. (4.47) 



For the gluon matrix element. Fig. 4-2(f) again does not contribute. For Fig. 4-2(e)| 



starting from the third line of Eq. (4.37), we find 



{gn\e{u:)0^g^'^'{t,u)\gn 



(D.7) 



^ r(l + e)(e^V^)^(— - 2z 

as{fi)TF 



2tt 



Bisct>ol2iA,B,e) 



2tt 
as{fJ.)TF 
2^r 



z)e{l - z)r{l + e){e"'^iJ^y{l -2z + 2z^ - e) 



sinvre 9{t) ( z \^ 



1 



P,,(z)(-^+ln ^ 



7re(l - e) \\- z 
1 - z 



1+0(1-2) 
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The same discussion as for the quark matrix element above can be repeated for the gluon ma- 



trix element. The {\ / e)6{t)Pqg{z) term matches the IR divergence in the PDF in Eq. (D.l ). 
Since there are no further poles, no UV renormalization is required and the quark and gluon 



operators do not mix. The finite terms in Eq. (D.7) then determine the matching coefficient 



Iqg{t, z, fi), reproducing our previous result in Eq. (4.47). 
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Appendix E 



Perturbative Results 



In this appendix we collect perturbative results relevant for the Drell-Yan beam thrust cross 



section in Eq. (8.2) and the cross section for isolated Higgs production through gluon fusion 



in Eq. (9.5). These were included as appendices in Ref. [38 1 1168j . 



E.l Fixed-Order Results for Drell-Yan 



We summarize the results of Sec. 7.4 The one-loop hard function for Drell-Yan is given by 
the square of Wilson coefficients in SCET for which the relevant one-loop matching from 
QCD onto SCET was carried out in Refs. [33l I138j 



1 



{vl + al){vj + af) - 2QqVqVi{l - m\/q^) 



(l-m|/g2)2 + ^|r|/g4 



27r r 



In^ ^ + 3 In ^ 



+ 



7^ 
6 



(E.l) 



We included the prefactor from the leptonic matrix element, where Qq is the quark charge 
in units of |e|, vt,q and a^^q are the standard vector and axial couplings of the leptons and 
quarks, and mz and Tz are the mass and width of the Z boson. The beam thrust soft 
function can be extracted form the one-loop incoming hemisphere soft function [941 1157| 
and is given by 



Sf{k+,,s) = 5{k+) + 



27r 



+ ^ S{k-^ 
b 



(E.2) 
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Our one-loop results for the matching coefficients in the beam function OPE in Eq. (3.37) 
are for the quark beam function given in Eq. (4.47). 



E.2 Fixed-Order Results for Higgs Production 

For Higgs production through gluon fusion we will integrate out the top quark and hard 
off-shell modes in one step. This avoids the often used expansion mn *C rrit, but does not 
allow us to sum logarithms of rnj^/mf. In the narrow width approximation the virtuality 
of the Higgs is simply = m^, but we will keep q'^ general in the below discussion. 
In pure dimensional regularization the matching coefficient CggH{'nit,(f') is given by the 
infrared-finite part of the on-shell ggH form factor |10 | IllH I115| [T52] . At NLO this yields 



-In- 



2 9" 77r 



(E.3) 



where p = m\ / m| and 
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For /9 — )• (corresponding to rrit — s- oo) we have F(o)(0) = 1 and F(i)(0) = 5Ca - 3C7f. 

The gluon beam thrust soft function has Wilson lines in the adjoint rather than funda- 
mental representation. At one loop we can simply replace Cj? by Ca (Casimir scaling) in 



the quark result in Eq. (E.2) 



S%\k+,p.) = 5{k+) + 



asjp) Ca 
2tt 



Kk+/p) \n{k+/p) 



The one-loop coefficients for matching the gluon beam function onto PDFs are given in 
Eq. ([5^. 
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E.3 Renormalization Group Evolution 



The RGE and anomalous dimension for the hard function are [see Eqs. (7.82) and (C.IO)] 



7]^(^^ /^) = 2r*u3p(a^) hi ^ + 1h{ 



a. 



(E.5) 



The expansion coefficients of r*^gp(as) and 7^(q!s) are given below in Eqs. (E.15) and 



(E.18). Note that ] cusp — Tcusp so r*^gp — Tcusp. The RGE in Eq. ( |7.82[ ) has the standard 
solution 



where the functions K^{fj,o,fi), r/p(;Uo,/^) and are given below in Eq. (E.12) 



The beam function RGE is [see Eqs. (3.25) and (3.26)] 



d f 

x,fj,) = J dt' 7^(t - t', fi) Bi{t\ X, , 

-fUt,f^) = -2r^usp(«s) ^^o(^) +7M«s)<5(t) 



(E.7) 



and its solution is [261 [Ml ttSS [m] [see Eq. (3.28)] 



Bi{t,x,n) = J dt' Bi{t - t',x,^o) C/]5(t',/io,/u) , 



-2r/f(/io,^). (E. 



The beam thrust soft function is given in terms of •S'jV ■ by 



{k+,^) = I dktdk+ St^ikt, k+,f,) Sik+ - kt - k+) . 
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(E.9) 



Its RGE is easily obtained by integrating Eqs. (C.5) and (C.ll), 



(E.IO) 

-7i^(a,)-27lj(a,), 



whose solution is completely analogous to Eq. (E.8) 



^0 ^/^o 



-K's (("0,^) = -4:Kl^{fio, ^J.) + K^^j{fMo, fj,) , r]'^{no,n) = 4r/r(^o,/^) • (E.ll) 
The functions Kl^{fio,fi), r7p(/Uo,/u), K^{fiQ,fi) in the above RGE solutions are defined as 
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Expanding the beta function and anomalous dimensions in powers of as, 
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their explicit expressions at NNLL are, 
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Here r = as(^)/as(/io) and we have suppressed the superscript i on K^, r/p and FJ^. 



Up to three loops, the coefficients of the beta function |133t 1172] and cusp anomalous 
dimension [1271 fl^ in MS are 
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The MS non-cusp anomalous dimension for the hard function H'^^ can be obtained |4H 
I118j from the IR divergences of the on-shell massless quark form factor which is known to 
three loops |148j . Similarly, the anomalous dimension for H^^ can be extracted 00] from 
the gluon form factor, which has also been calculated to three loops |149j . 
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Denoting 7^ the coefficient of ttie 6{l — z) in the quark PDF anomalous dimension, Eq. ( 3.17| ) 
(which gives the non-cusp part of the anomalous dimension in the threshold limit z — )• 1), 
the factorization theorem for DIS at threshold implies that 7H+7J+7J = 0, which was used 
in Ref. [H] to obtain 7j at three loops from the known three-loop result for 7J |147j . This 



argument can repeated for a gluons. As we showed in Sec. 3.2, the anomalous dimension for 
the beam function equals that of the jet function, 7^ = 7j, so the three-loop result for 7^ 



together with Eq. (E.18) yield the non-cusp three- loop anomalous dimension for the beam 
functions 
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At NNLL, we only need the one- and two-loop coefficients of 7^ and 7]^. The three-loop 
coefficients, 7]^ 2 and ^b2i given here for completeness. They are required for the 
resummation at N^LL, where one would also need the four-loop beta function and cusp 
anomalous dimension, the latter of which is has not been calculated so far. In addition, the 
full N^LL would also require the two-loop fixed-order corrections, which are known for the 
hard function, but not yet for the beam and soft functions. 
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